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The Effect of Changes in Altitude on the 
Controlled Behavior of a Gas-Turbine Engine 


SEYMOUR C. HIMMEL* ano RICHARD P. KREBS?t 
Lewis Flight Propulsion Laboratory, N.A.C.A. 


SUMMARY and propeller torque is a function of blade angle and 
For the various phases of flight operation of an aircraft pow engine speed 
ered by a gas-turbine engine, a considerable range of power out- 
put is required of the engine. Not only are different power levels Q, = Q,(8, N) (2) 
required, but changes from one power level to another must be 
accomplished in a short period of time. It is the function of the The only energy storage units are the engine and pro- 
control to conduct the engine quickly and safely from one oper- peller inertias. The engine operates dynamically in a 
ating level to another without encountering excessive values of quasi-static state. 
speed, torque, and temperature. The design of a control system 
for a gas-turbine engine is further complicated by the fact that 
the reaction of the engine to a given change in one of the inde- 
pendent variables changes with the initial control setting and 
with altitude. As a result, the control is expected to regulate 


what amounts toa variety of power plants. (I, + I,) d AN 4 | (oe) + (*) | 
wre B 


This paper presents a theory of the change in engine response dl oN oN 


The relation between engine speed and fuel flow with 
constant propeller blade angle, based on the foregoing 
assumptions, can be expressed by the following: 


analysts, together with some experimental data corroborating 
these analyses. The data show that the response of the engine OW; 
at any altitude may be computed from measurements made at 
sea level. This equation relates the time derivative of speed 
The sea-level and altitude transient behavior of an engine change, and the speed change itself, with the change in 
operating under a proposed control are studied with the aid of an : , : ‘ 
fuel flow, and J, is the moment of inertia of the engine 


(3) 


with altitude which has been derived independently by several 
od, 7 
AW, 
N 


electronic analog computer to illustrate the effect of the changes : ‘ ; 
alone and J, is the moment of inertia of the propeller 


pears necessary to use compensation to vary the gain and time alone. The term multiplying the change in speed is 
variant parameters of the control loops with altitude if satisfac- made up of the sum of two partial derivatives—the 
wy pertormance is to be obtained changes of torque with speed at constant fuel flow and 
at constant blade angle. The term modifying the 
change in fuel flow is the partial of torque with respect 
F* IM THE ANALYSIS OF THE EQUATION describing the to fuel flow at constant speed. From this differential 

dynamic behavior of a turbine-propeller engine, we equation the so-called ‘‘response’’ of engine speed to 
will develop the variation of engine dynamic perform- fuel flow at constant blade angle can be written in oper- 
ance with changes in altitude. The analysis is based ational notation as given in Eq. (4). 


in engine dynamics with altitude. From these studies, it ap- 


ANALYSIS 


on the following assumptions: For small excursions P ow 
from a point of equilibrium, engine torque is a function ( AN ) _ (00./OWs)w 
8 ( 


AW, re) ), {Oo iE 
“™ [) +) ] 
ON We ON B 


l 
Presented at the Flight Propulsion Session, Nineteenth Annual +1 
e DP 
[ + p| 


Meeting, I.A.S., New York, January 29-February 1, 1951. 

* Aeronautical Research Scientist, Controls Branch, Engine (OQ, ON) w, + (OQ, ON), 
Research Division 

t Aeronautical Research Scientist, Engine Research Division 


£33 


of fuel flow and engine speed 


Q. = Q.(W;, N) (1) 
(4) 


where the AN and AW, are, respectively, the change 
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RESPONSE TO A STEP FUNCTION FREQUENCY RESPONSE 
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Fic. 1. Theoretical response of speed to fuel flow. 


in speed and the change in fuel flow. The symbol p is 
a differential operator. ; 

The right-hand side of this equation can be considered 
as the product of two factors—a gain or sensitivity 


factor 


(00./OW))~ 
[(0Q./ON)w, + (0Q,/ON)g] 


and a time variant factor 


l 
(I. + I,) 
[(0Q./ON)w, + (0Q,/ON)g] 


The gain factor tells how large a change in speed results 
from a given change in fuel flow; the time dependent 
factor tells in what manner and how fast this change is 
effected. This equation gives the response of engine 
speed to fuel flow and is independent of the nature of 
the change in fuel flow. 

The dynamic response of a system can be determined 
in a number of ways. One method is to apply a step 
change to the independent variable of the response (in 
this example fuel flow) and to study the time history 
of the dependent variable (in this case speed). Were 
this to be done, the response of speed to a step change 
in fuel flow would be as shown in the left-hand portion 
of Fig. 1. In this figure we have plotted only the time 
variant or ‘‘how fast’’ portion of the response. This is 
equivalent to dividing the response by the gain factor 
or “normalizing” the response. It is seen that the re- 
sponse is zero at zero time and approaches its final 
value of unity exponentially. A second method of 
determining dynamic response is to apply a constant 
amplitude sinusoidal oscillation over a series of fre- 
quencies to the independent variable (fuel flow) and to 
measure the speed oscillations. The normalized re- 
sponse as determined from such an experimental pro- 
cedure is shown on the right-hand side of Fig. 1. When 
the frequency of the applied sinusoidal fuel flow is low, 
the speed follows the change in fuel flow yielding a 
normalized response close to unity, but as the fre- 
quency increases, the speed change is reduced in ampli- 
tude, for a constant magnitude of fuel flow change, and 
the response decays. 

A system that yields responses such as those shown 
is quantitatively described by a single entity known 
as the time constant. This time constant is the time 


1+ p 


for the speed to reach 63.3 per cent of its final value in 
the case of a step change in fuel flow or, in the case of 
the sinusoidal input, the reciprocal of the frequency at 
which the amplitude of the response falls to 70.7 per 
cent of its zero frequency value, the so-called break 
frequency. In the particular system under consider- 
ation, the time constant is given by 


7 I.+1, 
"° QQ/0N)w, + (0Q)/0N), 


()) 


The foregoing time constant is descriptive of the re- 
sponse of the combination of engine and propeller and 
is therefore referred to as the time constant of the com- 
bination, 7,. 

Thus, the dynamics of the system can be determined 
by subjecting the engine to either a step change in fuel 
flow or a series of sinusoidal varying fuel flows. Other 
responses within the engine, such as speed to torque 
or temperature to propeller blade angle, can be obtained 
by proper choice of the dependent and independent 
variables. Some of the responses—for example, that 
of speed to torque at constant blade angle 
lent to cutting the shaft between the engine and pro- 
these components 


are equiva- 


peller and considering either of 
separately according to laws of rigid body dynamics. 

From the response of speed to torque at constant 
fuel flow and at constant blade angle, two new time 
constants are defined—namely, the time constant of the 
engine alone and that of the propeller alone. Any 
other response, whether between a dependent and an 
independent variable or between two dependent vari- 
ables, is defined by a gain factor and a combination of 
the three time constants. Therefore, knowledge of the 
three time constants is essential to the complete descrip- 
tion of the power-plant dynamics and to the design of a 
control system. 

When the flight altitude is increased, the relation be- 
tween increments of the engine variables is changed. 
For example, the change of torque for a given change 
of engine speed at constant blade angle is decreased 
because of the decreased air density at high altitude. 
The inertia of the components, however, is independ- 
ent of altitude. Hence, a change in component dy- 
namics is to be expected with changes in flight alti- 
tude. These changes in dynamics are twofold. There 
is a change in the speed of response and a change in the 
sensitivity of response. Quantitatively, the change 
in the speed of response can be evaluated by a con- 
sideration of the change in time constant. The time 
constant involves the relation between torque and speed, 
as is shown in Eq. (5), and, hence, the time constant 
will vary as the V 0/6, where @ is the ratio of ambient 
static temperature to sea-level temperature and 46 
is the ratio of ambient static pressure to sea-level static 
pressure. In other words, if the time constant is multi- 
plied by a generalization factor 6/+/@, the product is 
invariant with altitude. This product is referred to as 
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OF ALTITUDE 


EFFECT 


generalized time constant. This variation of time con- 
stant with altitude is perfectly general and is appli- 
cable to all time constants involved in the responses. 

Consideration of the quantities involved in any given 
response will make it possible to determine the general- 
ization factor for a particular response sensitivity. No 
single generalization factor is applicable to all response 
sensitivities, but each can be considered individually. 
For example, the gain factor for the response of speed 
to fuel flow has the dimensions of speed over fuel flow. 
Because speed is generalized by the +/6 and fuel flow 
is generalized by 6+/ 0, the generalization factor for the 
gain term in the speed-fuel flow response is 6. Because 
it is possible to generalize both the sensitivity and the 
time constant individually, it should be possible to 
generalize the complete dynamic response. 

Thus far, no physical significance has been attached 
to the partial derivatives appearing in the response 
equations. These partial derivatives have the value of 
the slopes of their respective steady-state performance 
curves. Thus, it should be possible to determine dy- 
namic response from the steady-state performance 
curves of the engine and the moments of inertia of its 
components. 

The foregoing is an analytical approach to the de- 
termination of the dynamics of the engine. The direct 
method for determining the dynamic response and its 
variation with flight altitude is to subject the engine to 
controlled changes in the independent variables and 
to record the time history of the dependent variables. 
Such an investigation has been conducted in the Alti- 
tude Wind Tunnel of the N.A.C.A. Lewis Flight Pro- 


pulsion Laboratory. 


EXPERIMENTAL PROGRAM 


The Lewis Altitude Wind Tunnel is of the return 
type with a 20-ft. diameter test section and is so de- 
signed as to permit the simulation of altitude condi- 
A turbine-propeller engine was mounted on an 
For 


tions. 
airfoil that spanned the test section of the tunnel. 
the dynamic response investigations the frequency- 
response technique was used, because, with the applica- 
tion of a sinusoidal oscillation to either fuel flow or pro- 
peller blade angle, all of the other engine variables exe- 
cute sinusoidal fluctuations permitting a simple deter- 
mination of the response of any one variable to another. 
Application of a sinusoidal change to either fuel flow 
or blade angle is experimentally feasible. Therefore, 
a separate fuel system was designed to impress sinus- 
oidal fuel flows of variable frequency and amplitude 
on the engine. 

Dynamic response tests were run at pressure altitudes 
of 10,000, 20,000, and 30,000 ft. at a fixed tunnel air 
temperature. For the dynamic response runs, the 
engine variables were measured by highly sensitive fast- 
acting instruments and were recorded on a galvanomet- 
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Fic. 2. Experimental response of speed to fuel flow. 


ric oscillograph or strip chart recorders. These in- 
struments had sufficiently fast response so that no 
correction of the data for instrument dynamics was 


necessary except for the temperature data. 


RESULTS 


Some of the results of the dynamic response investi- 
gation are shown in the succeeding figures, wherein the 
amplitude response is plotted against the frequency of 
the applied sinusoidal fuel flow. In these figures the 
amplitude responses have been normalized and the fre- 
quency nondimensionalized. Normalization consisted 
of dividing the amplitude of the response by the gain 
or sensitivity factor for the response. The nondimen- 
sional frequency consists of the product of the fre- 
quency of the applied oscillations, w, and the sea-level 
time constant of the engine-propeller combination, 
TcSL+ 

The amplitude response of engine speed to sinus- 
oidally varying fuel flows with the propeller blade 
angle held constant is shown in Fig. 2. Responses are 
shown for three different altitudes. The general shape 
of the response curve is similar to that of a first-order 
dynamic system of the type shown in Fig. 1. It can 
be seen that the break frequencies for the three curves 
decrease with increasing altitude, indicating a slower 
response as the flight altitude increases. The time 
constant of the engine-propeller combination at 10,000 
ft. is approximately one and four-tenths times the value 
at sea level, while for 30,000 ft. the time constant is 
approximately three and one-tenth times that at sea 
level. 

In a similar manner, the speed response of the pro- 
peller alone and its associated time constant can be 
determined by a consideration of the speed-shaft 
torque relationship at constant propeller blade angle. 
The nature of this response is similar to that of the 
response of the engine-propeller combination—that is, 
a first-order lag system—and the response is affected 
in the same manner by changes in altitude. 


The dynamics of the engine alone and its associated 
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Fic. 4. Generalized experimental responses of temperature and 
speed to fuel flow. 


time constant may be determined by a mathematical 
synthesis of the two aforementioned responses or by 
the application of the frequency-response technique 
using a fixed fuel flow and variable blade angle and 
considering the speed-shaft torque relationship under 
these conditions. The engine alone is also found to 
have a response similar to that of the engine-propeller 
combination. 

The foregoing responses have all dealt with speed as 
the dependent variable. Another dependent variable 
of considerable interest from the standpoint of engine 
longevity and control application is the turbine-outlet 
temperature. A nondimensional amplitude response 
of turbine-outlet temperature to sinusoidally varying 
fuel flows at fixed blade angle is shown in Fig. 3. These 
data have had the thermocouple lag removed. A fre- 
quency response of this nature can be described by a 
system whose response equation is shown in the left- 
hand inset. Such a system arises from the combination 
of a lead network with a lag network, where the time 
constant of the lag is that of the engine-propeller com- 
bination and the time constant of the lead is some multi- 


TICAL SCIENCES JULY, 19391 
ple a times this time constant. For the case shown, the 
value of a is greater than unity and is equal to the ampli 
tude of the normalized frequency response at infinite 
frequency. The response of such a system to a step 
input is shown in the right-hand inset, and it can be 
seen that the dependent variable overshoots its final 
value. It can be shown that the ratio of the peak value 


to the final value is the constant a in the response equa 
tion. 

the speed-fuel flow response, the 
toward the 


As in the case of 


progression of the temperature response 
decreasing frequencies with increasing altitude indi 
altitude. The high 


the response is approximately the 


cates increased response times at 


frequency value of 


same for all altitudes. This indicates that the per 


centage overshoot for a step change in 


fuel flow is independent of altitude. 


temperature 


The data of the preceding responses for all three alti 
tudes have been replotted in Fig. 4 against a nondi 
mensional abscissa generalized for altitude in the man 


ner discussed in the ‘‘Analysis’’ section. Shown in 
Fig. 4 are both the temperature and speed responses ti 
fuel flow. It is seen that the application of the general 
ization factor has rendered the responses invariant with 
altitude 


sponse at any 


as predicted by theory. The dynamic re 
altitude can be predicted from the dy 
namic response at sea level. 

It is of interest to compare the values of the time 
constants involved in the preceding responses as deter 
mined from the experimental dynamic 


vestigations with those that may be computed from the 


response 1n 


moments of inertia and the slopes of the steady-state 
engine performance curves. The comparison is shown 
in Fig. 5. Here, the data on the left-hand side for the 
time constants of both the engine-propeller combination 
and the propeller alone and the overshoot ratio for the 
temperature-fuel flow response were determined directly 
from the frequency response investigations. The time 
constant for the engine alone was determined by syn 
thesis. It is interesting to note the relative values of 


the three time constants. The time constant of the 
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EFPFECT OF ALTITUDE 
combination falls between that of the propeller alone 
and that of the engine alone. The propeller time con 
stant is approximately one-third that of the combi- 
nation, and the engine time constant is about two and 
one-half times that of the combination. 

From the appropriate slopes of the generalized steady 
state performance curves and the moments of inertia of 
time constants on 


An examination 


the components, the generalized 
the right-hand side were computed. 
of the two sets of values show that they agree within 
10 per cent of one another. The values for the com 
bination are extremely close, while the greatest dis 
crepancy appears in the determination of the general 
ized time constant for the propeller alone from the data 
taken at 10,000 ft. As a result of this good agreement 
over the range of altitudes investigated, it is concluded 
that the steady-state performance data afford a reli 
able means of determining the dynamic performance 


of the engine. 


CONTROLLED BEHAVIOR 


The significance of the changes of engine dynamics 
with changes in altitude can be illustrated by consider- 
ing engine behavior when acted upon by a linear con 
trol system. <A linear control system was chosen so 
that effects of changes in engine dynamics would not 
be masked by control nonlinearities. In order fo study 
these effects, we have computed the behavior of a com 
plete controlled system at two different altitudes using 
an electronic analog computer. For this purpose we 
have chosen a control wherein the change in fuel flow is 
regulated in proportion to the sum of the speed error 
and its time integral and the change in propeller blade 
angle is regulated in proportion to the sum of the shaft 
torque error and its time integral. Such a control is 
referred to as a double-loop proportional-plus-integral 
control. The complete control system is shown dia 
grammatically in Fig. 6. 

The engine is represented diagrammatically by the 
large block. It has inputs of fuel flow and propeller 
blade angle and outputs of speed, torque, and temper 
ature. The speed is fed back and compared to the 
speed setting, and the difference or speed error energizes 
that portion of the control which regulates the fuel 
flow. Included in the operational representation of 
the control is the time delay of the fuel system, repre- 
sented by a first-order lag of time constant r,,;. Ina 
similar manner, the torque error controls the blade 
angle. A time lag in the propeller has been included 
in this loop, represented by a first-order lag of time 
constant rT, ». 

Such a controlled system possesses what is known 
as interaction that is, because each of the dependent 
variables is influenced by both of the independent vari- 
ables, a change in the setting of one of the control loops 
will cause the independent variable in the other loop 


to change. Interaction can be studied by representing 
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lic. 6. Schematic diagram of engine-control combination 


the engine by a combination of transfer functions 
representing the response of a dependent variable to 
one of the independent variables, as was done in the 
diagram. Included in the diagram is a representation 
of the turbine-outlet temperature as a function of the 
two independent variables fuel flow and propeller 
blade angle 

In the interest of simplicity of engine operation, we 
have made the control “single lever’’ operated by sched 
uling the speed and torque setting to ensure steady 
state operation of the engine along a predetermined 
operating line and have incorporated some form of 
barometric compensation in the scheduling to keep the 
engine on the operating line as the altitude is changed. 

The engine gains or sensitivities have been chosen as 
those determined from the experimental investigation 
of the turbine-propeller engine. The time constant of 
the engine-propeller combination has been arbitrarily 
selected to be | sec. at sea level. The other time con 
stants were chosen in the same proportion to the com 
bination time constant as those determined in the ex 
perimental investigation. 

The control has four adjustable parameters. 
are: the two control gains, A; and Ke, and two integral 


The magnitudes of these 


They 


time constants, 7; and To. 
parameters were chosen to give a good response at sea 
that is, fast response in the speed and torque and 
A sea 


level 
tolerable temperature and fuel flow overshoots. 
level response that was considered to be desirable is 
shown in Fig. 7. On this and succeeding figures the 
variables shown are speed, torque, temperature, fuel 
flow, and blade angle. The abscissa for all figures rep 
resents a 20-sec. time interval, and the final value of 
all variables represents the steady-state value corre 
sponding to the final value of the power setting. At 
sea level the change in power setting consisted of a 
100-r.p.m. change in speed and a 215-Ib.-ft. change in 
It can be seen that the speed and torque 
The first speed- 


torque. 
quickly come up to their final values. 
pip is reached in 0.5 sec., whereas the torque hits its 
peak value in 0.2 sec. and is practically constant after 
There is a temperature overshoot of 86 per 
This temperature over- 


11/4 sec. 
cent which lasts for 0.7 sec. 
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shoot amounts to approximately 45°R. above the final 
Such a fast speed response requires a 68 per 
This overshoot is not con- 


value. 
cent overshoot in fuel flow. 
sidered of sufficient magnitude and duration to precipi- 
tate burner blowout. The change of blade angle is 
small, only about 1.7°, and the fastest response de- 
manded by the control is approximately 5° per sec., a 
figure easily met by today’s propellers. 

Let us now take our engine and control up to 40,000 
ft. (Fig. 8) and observe the response to an equivalent 
change in control lever setting. The change of speed 
is still 100 r.p.m., but, because of the barometric com- 
pensation in the scheduling device, the change in 
torque is reduced to 40 Ib. ft. to keep the engine on the 
steady-state operating line. There is still a fast re- 
that is, the speed crosses its final value 
There are, however, 


sponse of speed 
in 0.3 sec. with a slight overshoot. 
intolerable overshoots in torque, temperature, and fuel 
flow. The severe torque overshoot endangers the re- 
duction gearing. The temperature overshoot stresses 
the turbine blading, and the fuel flow overshoot might 
very well blow out the combustors. Such excesses 
are not to be tolerated because they represent condi- 
tions severely endangering the engine. This poor re- 
sponse is the result of mismatch between the engine and 
assumed control, because the dynamics of the engine 
have changed with the change in altitude while the con- 
trol has retained its sea-level characteristics. 

The foregoing results have shown the need for some 


SCIENCES-—JULY, 1951 
compensation of the control to match the change in 


It has 
been previously stated that the change in the engine 


engine performance with a change in altitude. 


response consists of a change in engine sensitivity and a 
change in time of response. Likewise, the control has 
a sensitivity parameter and a time response parameter, 
It is the change in engine sensitivity which is primarily 
This 


change in engine sensitivities has long been recognized 


responsible for the large overshoots in Fig. 8. 


by control designers, and this recognition has resulted 
in compensation being built into the gains of the con- 
troller. 

The purpose of the control-gain compensation is to 
keep the loop gain of the system constant at all alti- 
This loop gain is the product of the engine and 
When the 
gains of both loops are kept constant by means of a 


tudes. 
the control sensitivities within the loop. 


barometric compensation on A, and Ko, the response of 
the controlled engine at a flight altitude of 40,000 ft. is 
as shown in Fig. 9. 

The gain compensation has succeeded in greatly re- 
ducing the torque, temperature, and fuel flow over- 
shoots. The fuel flow and temperature overshoots 
have now regained values comparable to those at sea 
level, but the compensation has rendered the responses 
The damping in both 
Two 


of all the variables oscillatory. 
loops is of approximately the same magnitude. 
cycles of oscillation are observed in all variables, and 
about 9 sec. are required before the variables retain 
their final values. 

The application of barometric compensation to the 
gains of the control produced an oscillatory system. 
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EFFECT OF ALTITUDE 
This response is the result of applying compensation to 
only one of the two sets of control parameters. The 
application of barometric compensation to the integral 
time constants, as well as the gains of the control, pro- 
duces a much more to be desired engine behavior at this 
altitude, as is shown in Fig. 10. All oscillations in 
speed and torque have been removed, and the tem- 
perature and fuel flow overshoots are even smaller 
than those encountered at sea level. The speed of re- 
sponse has been considerably reduced below that of 
the controlled response at sea level. The torque 
reaches SO per cent of its final value in 0.4 sec. but re- 
quires almost the full 20 sec. to attain its final value. 
The speed reaches 94 per cent of its final value in 5.5 
sec. Therefore, with a completely barometrically 
compensated control, we have an oscillation-free sys- 
tem in which the speed of response is three times as 
fast and the temperature overshoot one-third as large 
as for the uncontrolled engine at this altitude. 

It may be expected that the results shown in Fig. 10 
should be similar to those shown for sea level (Fig. 7) 
with the exception that the time response should be 
approximately V 0/6 or four and six-tenths times as 
However, compensation has not rendered the 
One reason for the 


long. 
two systems exactly proportional. 
lack of similarity is that barometric compensation does 
not compensate for the change in temperature with 
altitude. Secondly, reference to Fig. 6 shows that, 
when the gains of the individual loops are maintained 
constant, the contributions of the interaction do not 
retain their same proportional magnitude in the system. 

Let us look at the effect of compensating A, baro- 
metrically. A, acts on those portions of the engine 
whose sensitivity is designated by Aw and Kew. Ay 
is the sensitivity of the speed-fuel flow response, and 
Kew is the sensitivity of the torque-fuel flow response. 
While A, has been decreased by the factor 6 to com- 
pensate for the change in the factor Ay, which in- 
creased by the factor 1/6 due to the change in alti- 
tude, the factor Kgyw increased by only 1 V0. There- 
fore, barometric compensation of A, does not properly 
compensate for the change in Agy. Thus, even though 
barometric compensation does not preserve absolute 
similarity of response at all altitudes, the resultant en- 
gine behavior is satisfactory. 


SUMMARY 


We have arrived at a number of facts concerning the 
dynamic behavior of the engine. Specifically, we have 
shown that the dynamic behavior of the engine varies 
both as to sensitivity and time of response with changes 
in altitude. The variation of response with flight alti- 
tude is predictable, in that it was possible to generalize 
the experimental responses in a manner similar to that 
used to generalize the steady-state performance of the 
engine. Further, it has been shown that it is possible 
to analytically determine the dynamic responses of the 
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engine, within a reasonable margin of error, from the 
steady-state performance curves of the engine and the 
moments of inertia of its components. 

Applying this information, we have demonstrated 
the effect of the change in flight altitude on the con- 
trolled behavior of the engine, indicating the serious 
mismatch of engine and control which can arise from 
the changed engine dynamics at altitude, and the need 
for compensating the control to accommodate the 
changes. 

Gain compensation of the control, while removing 
some of the undesirable features of the uncompensated 
control action, gave rise to another objectionable fea- 
ture-—an oscillatory system. This situation was allevi- 
ated by barometric compensation of the integral time 
constants to account for the changed engine response 
times. This compensation removed the oscillations 
from the system at the expense of the speed or response. 
The slower response is, however, more rapid than that 
of the uncontrolled engine and is acceptable. 

The preceding discussion has been limited to the tur 
bine-propeller engine. However, the basic results and 
the method of approach indicated are equally applicable 
to any gas-turbine engine. 


APPENDIX (A) 


SYMBOLS 


temperature overshoot ratio 


I = polar moment of inertia (referred to engine speed ) 

Ku = gain of speed-fuel flow response at constant blade 
angle 

Kg = gain of speed-blade angle response at constant fuel 
flow 

Ke.g = gain of torque-blade angle response at constant fuel 
flow 

Ke.w = gain of torque-fuel flow response at constant blade 
angle 

Kr.w = gain of temperature-fuel flow response at constant 
blade angle 

Kr.g = gain of temperature-blade angle response at constant 


fuel flow 
N = engine speed 


p = differential operator d/d1 

Q = torque (referred to engine speed ) 

t = time 

= temperature 

W, = fuel flow 

B = blade angle 

A = incremental change from steady-state 

ry = pressure correction factor = ambient static pres- 
sure/N.A.C.A. standard sea-level pressure 

0 = temperature correction factor = ambient static 
temperature/N.A.C.A. standard sea-level tem- 
perature 

T = time constant 

w = frequency 

Subscripts 

c = engine-propeller combination 

e = engine 

| a = lag 

m = measured 

p = propeller 
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s = setting 
OL = sea level 
1 = fuel control loop 
2 = propeller control loop 


APPENDIX (B)—DERIVATION OF DyNAMIC-RESPONSE 
EQUATIONS AND TIME CONSTANTS 


The dynamic responses and the time constants are 
developed from the following assumptions: 

The engine torque is assumed to be a function of fuel 
flow and speed; propeller torque is assumed to be a 
function of propeller blade angle and speed. 


QO. = O.(W,, N) 
Q, = Q,(8, N) 


Similarly, any engine variable may be expressed as a 
function of two other variables. 

More explicitly, Eqs. (B1) and (B2), when linearized 
for incremental excursions about a steady-state point, 


(Bl) 
(B2) 


can be written as 


AQ. = (00,/0W)y AW, — (00,/0N)w,AN (B3) 
AQ, = (2Q,/08)v 4B + (0Q,/ON)g4N (B4) 


where the subscripts on the partial derivatives indicate 
the constant factor and where the signs are determined 
by the physical nature of the partial derivative. 

The difference between engine torque and propeller 
torque accelerates the engine; therefore, referring all 
terms to engine speed by multiplication by the ap- 
propriate function of the gear ratio, 


0. - Q = U, + I.) (dN/dt) 


If only incremental changes in the torques are consid- 
ered, Eq. (B5) may be written as 


AQ, — AQ, = UI, + I.) (dAN/dt) 


(BD5) 


(B6) 


Substitution of Eqs. (B1) and (B2) in Eq. (B6) yields 


20. _ (222) 7 202) 
= asi aa), ee , + 


2) | : dAN 
AN= (, tZ) 
be wy add dt 


(B7) 


At constant propeller blade angle, substituting p for 
d/dt in Eq. (B7) and rearranging terms yields the re- 
sponse of speed to fuel flow 

(“*) 2 (0Q./O0W,)x 

AW // [(0Q,/ON), + (020./ON)w | 


l 
( ) = Ky ( : ) (BS) 
l + TP l + TP 


where the time constant, 7,, is given by 


. 1 -+- f, = 
[(0Q,/0N) 4 + (OQ, ON) wl 
This time constant is descriptive of the dynamics of the 
engine-propeller combination and is therefore referred 
to as the time constant of the combination. 


(B9) 


Te = 
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In a similar manner, at constant fuel flow, the re- 
sponse of speed to blade angle is 


(=) = — (00, 08)x 
AB/ wy [(0Q,/ON)g + (02./0N)w,] 


f/ 4 l 
= K 3 ) 
(, 4 J ” (, + —) ata 


It is to be noted that the time constants in Eqs. (B8) 
and (B10) are the same, the difference in the responses 
being the gain or sensitivity term. 

The dynamics of the engine alone and of the pro- 
peller alone may be determined by introducing a new 
variable, shaft or measured torque, Q,,, and consider- 
ing the response of speed to the variable. This is equiv- 
alent to cutting the shaft between the engine and pro- 
peller and considering either of these components sep- 
arately according to the laws of rigid body dynamics. 

The difference between engine torque and _ shaft 
torque accelerates the engine; thus, 


Qe — Om = I-(dN/dt) (B11) 


Considering only incremental changes in engine and 
shaft torque, Eq. (B11) may be written 


AQ, — AQmn = I.(dAN/dt) (B12) 


Substitution of Eq. (B3) in Eq. (B12) vields 


oQ, ) (2) ee 
AQn = mw, =~) ant 
e ed soos ON] w, dt 
(B13) 


At constant fuel flow, substituting p for d/dt and re- 
arranging terms yields the response of speed to meas- 
ured torque. 


(=) = ! ; ) (BI4 
AQn 7 (OQ, ON) w, ( —- T oP 


Wy 
where the time constant involved is given by 


tT. = I,/(0Q0./0N)w, (B15) 


Eqs. (B14) and (B15) describe the dynamics of the 
engine alone. 

The difference between shaft torque and propeller 
torque accelerates the propeller. Considering incre- 
mental changes, this difference can be written as 


AQn — AQ, = I,(d AN/dt) (B16) 
Substituting Eq. (B4) yields 
AO, = (ser) Ag + (20) er (42%) 
can Op / ON /. e dt 
(B17) 


At constant blade angle, this expression can be reduced 
and rearranged to yield the response of speed to meas- 
ured torque: 


( AN l | ) 
) e , ( (B18) 
AQm 8 (002, ON )g l + TpP 
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SFRPECT OF ALTITUDE 


where the time constant, r,, is given by 


ty = 1,/(0Q,/ON)g (B19) 


Eqs. (B18) and (B19) describe the dynamics of the pro- 
peller alone. 

A variable of importance in the control of an engine 
is some characteristic temperature downstream of the 
combustion chamber. The of that tem- 
perature to the independent variables of the engine 


responses 


may be derived as follows. 
Temperature is assumed to be a function of fuel flow 
and speed. 
T = T(W,, N) (B20) 
Considering only incremental changes of the variables, 
Eq. (B20) can be written as 
AT = (OT/OW,)y AW, — (OT/ON)w,AN (B21) 
At constant propeller blade angle, the relation between 
speed and fuel flow is given by Eq. (BS). Substitution 
of Eq. (BS) in Eq. (B21) yields the response of tem- 
perature to fuel flow at constant blade angle. 


~) . ‘ba 
ey , Os 
| eS (22) 
ON] w AOW,/ x 


(B22) 
| 2G) . (22) | (1 + rp) 
Fe B ON / wy 
Eq. (B22) can be written in the form 
AT : 1 + ar.p 
( -) = Krw ( ore) (B23) 
AW, B ] + T-P 
where 
(-) (= ) 
= or ON/ wy OW s/n 
Ar.w = . se (B24) 
OW;/ x | (se) + 2“) | 
ON B oN Wy 
and 
a = (O7/OW,)y/Kr. x (B25) 


At constant fuel flow, the relation between speed and 
blade angle is given by Eq. (B10). Substitution of Eq. 
(B10) in Eq. (B21) yields the response of temperature 


eed (2) 
a] - ON Jw A 08 /y ‘ 
ably, [/2Q, =) 
(se) + (35), 
| 
| | = Kr. 2 ( ) (B26) 
(1 + tp) 1 + rp 


to blade angle. 
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For the analysis of the controlled system, it is neces- 
sary to know the response of the measured or shaft 
torque to the independent variables. This may be 
done in the following manner: 

The response of measured torque to fuel flow at con- 
stant blade angle may be obtained by eliminating the 
AN between Eqs. (BS) and (B18) in this manner: 
(B27) 


[AQn/AWy]p = [AQn/AN],[4N/ AW,]q 


Substituting the responses of Eqs. (BS) and (BIS) in 
Eq. (B27) yields 


( 20r) ( 00. ) 
Fa . ON/,\O0W7/ wn 
AW, Je | (Se) ; () | 
, v y 
ON B ON Wy 
(1 ) l 
| = aA = Ke v( A ) (B28) 
(1 + tp) 1+ tep 
Similarly, the response of measured torque to pro- 


peller blade angle may be obtained by eliminating the 
AN between Eqs. (B10) and (B14); thus, 
[ AQm ABlw, —= [ AQm AN lw, AN/ AB lw, ( B29) 


Substitution of the appropriate responses yields 
(22) (2) 
| A= sa ON wr op N x 
wh 12) + @).] 
ON/ ON/ wy 
l l eP 
( + 4 = Kos ( z = (B30) 
1 + rp 


1 + r.p 
The foregoing responses and time constants may all 
be presented in a generalized or ‘“‘corrected’’ form. 
For example, the response of speed to fuel flow may be 


written in generalized form as 


| a \ *| 7 
AW, ON 6 B 


( OV,/6 ) 
OW, /50/ 6) nro 1 
o0,/6 ) ( dQ, /6 ) ‘| ( 5 
l 'e ? 
(Ss V6 Meg ON/V/0 Wy/5% | dle V6 ! ) 
(B31) 


° ° e a . 
where the generalized time constant, 7,(6/Y/@), is given 
by 


j i. a 
-_ = r (B32) 
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On Turbulent Fluid Motion 


MAX M. MUNK 
Naval Ordnance Laboratory 


ABSTRACT 


Turbulent fluid motion is a daily and common experience, and 
still it comes as a surprise to the inquisitive student of hydrody- 
namics. It caine to him as a surprise even in the early period of 
boundary-layer study at a time when the possibility of turbulence 
had already been well recognized. By now we are resigned to the 
frequent occurrence of turbulence but have still to suppress the 
feeling of slight mental frustration or incongruence whenever we 
give serious thought to it. In spite of half a century of intense 
mathematical effort, the frequent instability of laminar flow is 
not clearly understood, and turbulence has not been laid bare 
as a simple consequence of the kinematics and dynamics involved. 
It has so far been denied to our inquisitive minds to see turbu- 
lence clearly, distinctly, and intuitively as a direct consequence 
of the underlying laws of motion. 

There are presented in this paper some broad thoughts toward 
clarification of the turbulence problem. There are thoughts 
that may perhaps suggest to one or the other reader a more pre- 
cise and perhaps more heuristic inquiry. These thoughts do not 
fully satisfy. They carry the inquiry forward without bringing 
it to full conclusion. Still they set the mind a little bit at ease 
and blunt the edge of consternation and frustration. 

The arguments will be chiefly directed to the straight channel 
flows of the incompressible fluid, which are the plane Couette 
flow and the Poiseuille flow. These are the simplest cases, but 
the territory of inquiry is much broader and the thoughts pre- 
sented are indeed fairly general. Also, in the following, it will 
not be necessary to give much attention to the three-dimensional 
aspect of the flow; the present issue is basically two-dimensional. 
Actual turbulence is, of course, always three-dimensional. 

The following is directed to readers reasonably prepared for 
It seems therefore unnecessary to define com 
nor does it seem necessary to describe the 


the discussion. 
monly used terms; 
phenomenon in question. 


ENERGY CONSIDERATIONS 


_ TWO-PLANE CHANNEL FLOWS just referred to 
satisfy a minimum condition with respect to their 
work dissipation. The laminar channel flows do not 
involve the active manifestation of mass forces. The 
viscosity forces are in equilibrium with the pressure 
distribution. It has been shown by Lamb in //ydro- 
dynamics that, in consequence of this, the laminar un- 
disturbed channel flows have less work dissipation than 
any other flow consistent with the boundary condi- 
tions. 

A straight plane laminar flow not having the Couette 
or Poiseuille velocity profile is not in equilibrium. 
Aside from this, it does not have absolute minimum 
dissipation; this would contradict the proposition just 
It has, however, a relative minimum of dis- 

As long as the boundary conditions are 


made. 

sipation. 
complied with and the time average of the flow is main- 
tained (whereby the “‘basic’’ flow or velocity profile is 
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maintained), the dissipation is increased by the super- 
position of any “‘disturbance’’ flow periodic in the flow 
direction such that its space average is zero. This fol- 
lows directly from the fact that the rate of dissipation 
per unit volume is proportional to the sum of the 
squares of certain space derivatives of the velocity 
distribution and because, in this case, the total rate of 
dissipation is proportional to the space integral of the 
vorticity square. 

The kinetic energy of the channel flow is clearly not 
a minimum, but the flow of the same flux having mini- 
mum kinetic energy has constant velocity. 

The work supplied to the flow can conveniently be 
figured out by using a reference system such that the 
total channel flux is zero. The channel walls move 
then in their own direction, and the work supplied per 
unit time is equal to the product of the wall velocity 
times the surface friction force. The latter is linear 
in terms of the vorticity at the It follows that 
the’ total energy supply is not changed by a purely 
it remains unchanged 


wall. 


space periodic disturbance, but 
as long as the flow itself does not change. 

From energy considerations it follows, therefore, that 
no disturbance can grow on an unchangeable laminar 
basic flow initially at equilibrium. That is to say, any 
instability of the steady laminar channel flow must in- 
volve more than a sole increase of the magnitude of the 
disturbance. It must involve a change of the basic 
flow itself or, in other words, a change of the basic 
velocity profile. Any increase of a space periodic dis- 
turbance is to be interpreted as being associated with 
a proper velocity profile change. This is already 
necessary during the first millionth of a second, other- 
wise the disturbance could not increase for want of 
The basic flow must change to one 
immaterial that the 


energy supply. 
supplying more energy. It is 
work deficiency involved is at first small, being of second 
order. The lack of even this small amount of energy 
would prove fatal and would invalidate the result. 

The linearized theory does not consider the energy 
balance or the velocity profile change. This hampers 
the clarification furthermore, 
doubts as to the significance of the results. As to the 
clarification, the linearized theory leaves a gap just 
where the greatest fullness is required. We look first 
and almost exclusively to the energy balance to settle 
If an initial 


process and, raises 


the question of stability or instability. 
tapping of an available energy source has to be care- 
fully avoided lest the energy stream opened cannot be 
stopped again, we recognize instability without too 
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ON TURBULENT 
much curiosity concerning the kinematics of the prob- 
lem. The kinematic aspect is not likely to satisfy us 
and becomes problematical if it is arrived at without 
regard to, or even in defiance of, the energy balance. 

We can easily see how the basic flow must change 
in case of instability for increased energy supply. 
There must be vorticity increase at the wall and, 
hence, a vorticity decrease at some region spaced from 
the wall, since the total vorticity is invariant. We are 
at this moment discussing the space average—that is, 
a definite basic velocity profile into which the initial 
undisturbed laminar velocity profile has been changed. 

Concentrating the vorticity near the wall tends to 
even out the velocity profile and is thus not inconsistent 
with a diminished kinetic energy of the flow. Thus, 
a certain amount of energy that may provide for in- 
creased dissipation becomes available. This effect is, 
however, not significant, and it is not important for 
explaining the occurrence of turbulence. This kinetic 
energy is finite and is available one time only, whereas 
continued turbulence requires a continuous supply of 
excess work to be continuously dissipated by excess 
dissipation. The kinetic energy is a secondary source 
of energy and cannot explain or support continuous 


excess dissipation. 


VoRTICITY DISCRIMINATION 


We have arrived at a basic flow, not necessarily in 
equilibrium by itself, associated with a whirl or dis- 
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The total sum of the dis- 


These whirls or disturbance 


turbance distribution. 

turbance vortices is zero. 
vortices deteriorate by viscosity, and new ones spring 
up from the basic flow and from the instability of the 
disturbance flow itself, so that with developed turbu- 
lence the disturbances remain statistically constant. 
Likewise, through viscosity effects, the basic flow tends 
to approach the laminar basic flow. A mechanism is 
therefore needed to keep the turbulent basic flow differ- 
ent from the laminar basic flow. The disturbance it- 
self must have such an effect; there is nothing else that 
may. The cardinal point of inquiry, or at least an 
important one, is, accordingly, a search for a mechanism 
that brings that about. On either side of the channel 
there is an equal number of disturbance vortices of 
both directions. The vortices corotational with the 
vorticity of the basic flow last longer near the walls and 
are shorter lived near the channel center than the coun- 
terrotational vortices, whereby the laminar basic flow 
is changed in profile and whereby said changed basic 
profile is conserved. We want to understand how and 
why that happens, for, if we do, we see a reason for 


turbulence. 


It isimmediately apparent that the linear terms of the 
equations of motion cannot have such a discriminative 
effect. The clarification has to lean on the quadratic 
terms. In view of the fact that the solution of the 
linearized equation took a lifetime, it seems to be an 
almost hopeless undertaking to try to solve the com- 
plete or correct equation. We shall, therefore, content 
ourselves with qualitative discussions of some kind as 
the only alternative at this time to giving up the prob- 
lem altogether. 


It will be important to find out whether any vortex 
discrimination operates diffusively or with local con- 
centration. If the discrimination leads merely to a 
change of the velocity profile without maintaining or 
bringing about a pronounced disturbance motion, it 
would not support turbulence. The increased energy 
supply would then be promptly dissipated in an orderly 
fashion by way of the changed velocity profile reverting 
to the initial laminar velocity profile. If, however, con- 
centration of certain vorticity whirls takes place, a case 
is made for instability of the laminar flow and for the 
creation or maintenance of the turbulent flow. For 
whatever the flow may be, the sum total of its vorticity 
is invariant. characterized by loval 
vorticity concentrations as contradistinguished from 
Local vorticity con 


Turbulence is 


a mere velocity profile change. 
centration is itself turbulence. 


THE STREAMWISE VORTEX RoW 


The quadratic character of the equations oi motion 
precludes superposition of individual solutions. The 
different aspects of the ensuing motion are not repre 
sented by particular solutions, and yet we must spe 
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cialize. We shall, therefore, study a special flow that 
we believe to be typical. 

Any disturbance vorticity concentrated within a 
small area will be associated with some opposite dis- 
turbance vorticity at another point. The channel 
flows we are considering are basically shear flows. 
They will stretch out any compact region of disturb- 
ance vorticity in a streamwise direction. The combina- 
tion of the basic flow and disturbance vorticity periodic 
in streamwise direction and approximately concen- 
trated along one undisturbed stream line appears ac- 
cordingly as a typical combination flow. Most of the 
disturbance is approximately nonrotational and can, 
therefore, be represented by a potential. The disturb- 
ance flow can be approximated by two different poten- 
tial flows—one on either side of the vortex line—of the 
type cosh y sin x. In addition, a boundary-layer 
flow will occur at each wall to eliminate slipping. The 
vorticity at the inner vorticity line per unit length is 
equal to the velocity jump of the two flows at their 
common boundary. It is supposed that the inner 
vorticity has already spread out by reason of viscosity 
and that it keeps on doing so similar to a boundary 
layer. 

So far, there is symmetry with respect to the two 
rotational directions. We try to concentrate our at- 
tention on possible vortex discrimination. The vis- 
cosity effect is basically a linear diffusion effect. Within 
an unlimited region, the square of the radius of gyra- 
tion of anything diffusing increases at a constant rate 
proportional, in the present case, solely to the kinematic 
viscosity having the proper physical dimension—area 
divided by time. If the region available for spreading 
out is bounded, naturally the square of the radius of 
gyration cannot continue to increase, and the rate of 
increase must be smaller. These are all details; the 
essential effect is that increase of said square takes place 
until positive and negative vortices merge, whereby 
they annihilate each other. Above all, diffusion is 
linear and does not discriminate between the two vor- 
ticity directions. 

The square of the radius of gyration of the disturb- 
ance vortex concentration would not remain equal 
even in absence of viscosity. Convection by the basic 
flow and self-convection by the disturbance takes place 
and changes said square. Here is a possible key to 
turbulence. The self-convection of our special dis- 
turbance flow exposes the disturbance vortices to con- 
vection by the basic flow whereby they spread out or 
move together (‘“cumulate’’), depending on whether 
the vortices are corotational or counterrotational with 
the basic flow vorticity. Here, we clearly have vor- 
tex discrimination, and this is a sufficient reason for 
subsequent change of energy supply. The inner vor- 
tex line is bent by self-convection into a wave, and 
the nodal points of the wave coincide with the points of 
maximum vorticity. Let the vortex line initially coin- 
cide with the stream line of zero velocity of a Couette 


AERONAUTICAL 


SCIENCES—JULY, 1951 

flow. After the waving has taken place, the vortices 
of corotational direction become flowwise compressed 
or telescoped (or cumulated), and the vortices of coun- 
terrotational direction become spread out in flow direc- 
tion by the convection of the basic flow, whereby the 
points of zero disturbance vorticity become pairwise 
and alternately more spaced and less spaced. Thereby, 
the diffusion effect of viscosity is aided by, or added 
to, with the counterrotational vortices, but it is coun- 
teracted with the corotational vortices. All disturb- 
ance vortices travel toward the channel walls, but the 
counterrotational vortices become spread out thinner 
than the corotational vortices. The counterrotational 
vortices will slow down and acquire relatively more 
area near the vortex line, and the corotational vortices 
will speed up and penetrate farther away from the 
vortex line and will acquire relatively more area near 
the walls. 

This takes place everywhere and, again, secondarily 
locally in smaller scale as an incident of a sufficiently 
large disturbance itself. The result can only be sta- 
tistically understood. The vortex discrimination does 
not cancel out statistically, being essentially a quadratic 
effect and having everywhere the same sign. It leads 
ultimately to increase of work supplied by the wall fric- 
tion and thus favors instability of the laminar flow or 
statistical equilibrium of the turbulent flow, provided 
the cumulation is large enough to exceed the effect of 
viscosity diffusion. It is a priori clear that the Reyn- 
olds Number of the flow determines whether cumula- 
Viscosity is dominating 
This calls for 


tion is relatively large or not. 
at extremely small Reynolds Numbers. 
stability of the laminar flow. At extremely large Reyn- 
olds Numbers the pressure modifications caused by 
viscosity are comparatively small. A much smaller 
scale of whirling is necessary to make the influence of 
viscosity on the pressure distribution comparable to 
the mass effect influence on the pressure distribution. 
Conditions are unfavorable for larger whirls. The 
smaller ones have more tendency to cancel out statis- 
tically because secondary cumulation becomes more 
dominant. This is again more favorable to stability 
of the laminar flow. Viscosity effects and mass effects 
in balanced combination are most favorable to cumula- 
tion at the right scale. At extremely large Reynolds 
Numbers there is, finally, not enough viscosity effect 
left with the Poiseuille flow for breaking up the line of 
fluid particles of maximum vorticity along the wall. 
There must be a limit to cumulation by vortex dis 
crimination. After a certain intensity of turbulence 
has been established, a statistical equilibrium is reached 
Cumulation 
Vorticity 


and turbulence no increases. 


reaches saturation because it is not linear. 


longer 


in shear flow gives rise to less cumulation if the vor- 
ticity is concentrated along a narrow channel strip 
than if it is evenly spread out through the entire chan- 
nel. Most of the turbulent channel flow has only 
small shear, and the flow region of large shear becomes 
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ON TURBULENT 
narrow. The increasing cumulation required becomes, 
finally, equal to the decreasing cumulation taking 


place. 


DIFFUSIVE SPREADING OUT 


The following vector analytical equations are written 
with the use of the positive and negative signs as orig 
inally proposed by Hamilton. The equation of dif- 


fusion of a scalar quantity F appears then as 
Fy, = wW-VF (1) 


Vv F denotes the antigradient—that is, a vector equal 

and opposite to the vector conventionally designated 

as the gradient. V -V F designates the Laplacian with 

the conventional sign, so that an ordinary three-dimen- 

sional maximum has a negative Laplacian. Eq. (1) 

indicates, therefore, a smoothing out or leveling out. 
The space integral; 


S ds 


r-rF (2) 
represents the negative moment of inertia of the dis 


tribution F. We assume that the space integral 


JS dS F (4) 
is finite and that the space in which diffusion 
occurs is uninterrupted and unlimited. We _ ex- 
amine the time change of the negative moment 


of inertia of F with respect to the origin of the refer 
Applying Green’s theorem two 
times in that is, applying two times 
Gauss’ theorem by way of integration by parts 


d . ‘ 
— aSrer F =» dor-rVF—- 
al . 
2 do-rF+ 6 / ds F) (+) 


a 


ence system chosen. 
succession 
gives 


The space integral (3) is constant provided the surface 
integral {“do-V F approaches zero with increasing space 
inclosed by the surface 0. Hence, if all surface integrals 
occurring vanish as the smallest distance of the surface 
from the origin approaches infinity, there results 


d/d7 (fdS r-rF S ds F) = Ov ()) 


That is to say, the square of the radius of gyration of 
the quantity F increases uniformly at the rate 6v. In 
the two-dimensional space, the equivalent rate of in- 
crease is only 4», and in the one-dimensional space, it is 
only 2p. 

If the space available for the diffusion is finite, there 
is naturally an upper limit to the square of the radius 
of gyration. The rate Gv is then the upper bound of the 
rate of increase. 


QUANTITATIVE ESTIMATE 


~ 


We proceed to a rough quantitative estimate of the 
intensitv of the cumulation effect to compare it with the 


FLUID MOTION 445 
diffusion. We assume a vortex row periodic in the 
stream direction X situated parallel to the wall of a 


Couette flow at a distance therefrom equal to 2.5 


length units. The channel width may be 10 length 


units. The disturbance flow potentials ®; and #, may 
be 
®, = sin x cosh (5 — y) 
, sinh 2.5 (6) 
@, = sin x cosh (5 + y) ——— 
sinh 7.5 


The cross velocity component O® Oy of the two poten- 
tial flows agrees at the vortex line and is equal to 


I” = sin x sinh 2.5 (7) 
The streamwise components do not agree at the vortex 
line, so that the combined flow is associated with the 
vorticity per unit length of the vortex line of intensity 
equal to the difference of that velocity component. 


cosh 7.5. y 
+ - sinh 2.; ) ze 


Q = cos x cosh 2.5 (1 : 
sinh 7.5 


2 cos x cosh 2.5 (8) 


Furthermore, the shear at the channel walls will give 
rise to suitable additional boundary layers, whereby 
slipping is eliminated. The vorticity [Eq. (8)] will also 
spread out in cross-stream direction so that they will 
soon occupy an interface layer of finite thickness simi 
lar to a boundary layer. 

The general features of such vortex row have already 
above. We direct our attention to 
The distance from the initial 


been described 
the cumulation only. 
vortex line attained by the vortices is proportional to 
the time, and the streamwise distance X traveled is 
proportional both to the cross distance that determines 
the streamwise time. In all, 
therefore, it is proportional to the square of the time. 
Cumulation is not a linear effect, but it is a quadratic 


velocity and to the 


effect relative to time. 

We proceed to compute the increase of the square of 
the radius of gyration of the corotational vortices 
caused by cumulation during a certain time interval 
7. The integral 


{a + 2x 
/ xsinxdx = 2x (9) 


a designating any arbitrary constant, will occur in the 
computation. 

The vortex of strength [Eq. (8)]| traveling with the 
velocity [Eq. (7)] has been displaced through the 
distance 


/ Vdt 


and has thereby acquired a streamwise velocity 


- [ova 


. 
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The streamwise displacement X is therefore equal to 
X= SS (1/2)w sin x sinh 5 dt dx 


The square of the radius of gyration has been decreased 
by 


°x/2 t t 
J w sin x cos x dx f i sinh 2.5 cosh 2.5 dt? 
see 0 0 


(11) 


The first integral in Eq. (11) has the value mw. The 
length of the row occupied by each group of corotational 
vortices is x. The chosen width, 10, of the channel is 
reasonably commensurate with the length mw. The 
difference of the two wall velocities is 10w, w designating 
the angular velocity of the cross fluid lines. The Reyn- 
olds Number of the Couette flow using the width and 
said velocity difference is, accordingly, 100w. The 
maximum velocity V given by Eq. (7) is then sinh 
2.5 = 6. An acceptable average velocity may, per- 
haps, be 4. An average travel distance for the vortices 
may be one-fifth of the width—that is, 2. This gives a 
time interval of travel equal to '/2, during which cumu- 
The total corotational vorticity consid- 
All this leads to a decrease of 


(10) 


lation acts. 
ered is 2 cosh 2.5 = 12. 
the square of the radius of gyration equal to 


mw(1/2)? 36/(1/2) = 187w (12) 


The viscous diffusive increase of the same square is 
perhaps 


dty = (3/2)y (13) 


Setting Eqs. (12) and (13) equal to each other leads to 
a Reynolds Number indicative of the critical one where 
stability has its limit. This gives 18mw/(3 2)y = 1, 
w/v ~ 1. Hence, R = 100. 

Considering the side portion of a turbulent Poiseuille 
flow as similar to a Couette flow, the Reynolds Number 
defined by the width of the Poiseuille flow is perhaps 
5 to 10 times as large as the number previously found 
that is, R = 500 — 1,000. 

The result of this rough estimate depends greatly on 
the choice of the cumulation distance and of the cumu- 
lation time. All that can be said is that a reasonable 
choice of both is conceivably consistent with the cor- 
rect Reynolds Number. 

The Poiseuille flow has variable vorticity, and the 
disturbance flow would disarrange that vorticity. This 
disarranging plays an important part in the linearized 
It may, by itself, contribute to instability, 
This dis- 


theory. 
but only if it is supported by cumulation. 
arranging would constitute a detail of the destabilizing 
effect of cumulation. 
THE LINEARIZED THEORY 
The ultimate cumulation-diffusion ratio 
independent of the initial strength of the disturbance. 
The linear kinematic effects are indeed dominant during 
the initial time elements, but even then they are vitally 


appears 
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dependent on the velocity profile changes. In the fina] 
result, diffusion and cumulation are side-by-side with 
equal order of magnitude, their relative importance 
not being dependent on the initial magnitude of the 
disturbance. This raises serious doubts about the sig- 
nificance of the results yielded by the pure linearized 
theory disregarding cumulation. The Reynolds Num. 
ber computed as critical is not necessarily critical. 
Linearized instability does not appear to be demon- 
strated to be either a necessary or a sufficient condition 
for actual instability. 

The linearization seems to have been resorted to by a 
routine striking out of all quadratic terms. These terms 
may, however, be of significance and may not even be 
small. Since the absence of a solution of the exact 
equation makes it impractical to examine the justifi- 
cation of linearization by mathematical procedure, such 
examination may be performed by the intuitive faculty. 
Linearization is open to great doubt when resorted to 
without any examination. 

The importance of cumulation is even greater in con- 
nection with the question of statistical equilibrium of 
the velocity profile of the turbulent flow. Linearized 
theory has never been applied to this question, although 
it may be if completely and broadly applicable. There 
is little doubt that this statistical equilibrium is closely 
allied with cumulation and vortex discrimination. 


CONCLUSION 


It is of significance that the instability of the laminar 
flow has been traced back to shear and not to vorticity. 
The distance between properly selected pairs of fluid 
A fluid body rotating like a 
Such motion is stable 


particles must diminish. 
solid has rotation but no shear. 
and cannot give rise to turbulence because it cannot 
have dissipation. The couple that would supply the 
energy to be dissipated would not be associated with an 
increase of the moment of momentum, in contradiction 
to Newton’s laws of dynamics. The mechanism de 
scribed accordingly not only makes us expect turbu 
lence of shear flows under suitable conditions but also 
meets the test of not predicting turbulence of the pure 
rotational flow. 

We have reached the end of the inquiry. We started 
out without seeing a mechanism that produces turbu 
lence, and now we see such a mechanism reasonably 
clearly. We see the disturbance 
vortices slowed down in their race with the corotational 


counterrotational 


vortices on their way toward the walls, because positive 


cumulation is reserved for the corotational vortices. 
This vortex discrimination leads to increased wall ro 
tation and wall friction and, hence, to increased energy 
supply. Nothing indicates that the excess energy su; 
plied is dissipated so uniformly as to exclude the form: 
tion of new disturbances. Thus, a foundation for con 
tinual turbulence is laid. The increasing velocity pro 
file steepness near the wall finally leads to saturation vi 


the vorticity and to statistical equilibrium. 
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Simplified Design Comparisons of Axial 
Compressors 


FRANK L. WATTENDORFY 


Air Engineering Development Division, US.A.F. 


ABSTRACT 


This paper reports a continuation of studies initiated in 1944 
regarding potential improvements of turbojet propulsion systems 
for supersonic flight as suggested originally by Dr. Theodore von 
Karman, Chairman of the U.S.A.F. Scientific Advisory Board 
and Chief Scientific Advisor to the Air Engineering Development 
Division 

Simplified design comparisons of axial compressors with par 
ticular emphasis on the attainment of a high rate of mass flow 
per unit frontal area are discussed. This feature is important 
not only in the application of turbojets to the propulsion of 
supersonic aircraft but also in the study of future potential re 
quirements for jet engine component testing and for application 
to the design of more compact compressor systems for wind tun 
nels such as those at the Arnold Engineering Development 
Center 

A comparison is made of some aspects of axial compressor de 
sign for accommodating a high rate of mass flow. In particular, 
the influence of hub ratio and axial Mach Number on the stage- 
pressure ratio based on assumed limitations of all relative and 
absolute Mach Numbers is given for compressor types designed 
according to various principles such as free vortex flow, rigid 
body prerotation, and arbitrary distribution of rotor work. 

The important relationship between maximum temperature 
rise and mass-fiow ratio is derived for several compressor types so 
that a genera! qualitative picture is given of the influence of com 
pressor type on the relationship between temperature rise and 
mass-flow rati This study represents a qualitative method of 
comparing compressor types and refers to the inlet stage alone, 
which is usually critical for mass flow. 


NOTATION 


= ratio of hub diameter to tip diameter 


I = area 

( = velocity of sound 

( = specific heat at constant pressure 

Cc = absolute velocity 

Cr = blade lift coefficient 

Cr = temperature rise coefficient 

g = gravitational constant 

J = mechanical equivalent of heat 

WV = }imiting intake Mach Number 

M.F.R. = mass-flow ratio 

“i = absolute pressure 

R = radius at tip 

4 = absolute temperature 

Al = stage temperature rise 

l = circumferential speed of blade clement 
Uy = tangential component of whirl induced by rotor 
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V = axial velocity 
WwW = relative velocity 
Subscripts 
aR = conditions at rotor hub 
ax = axial component 
cr = conditions at local sonic velocity 
m = mean velocity of intake and exit velocities 
0) = stagnation conditions before rotor 
r = conditions at any radius 
R = conditions at rotor tip 
u = tangential component 
1 = jntake conditions of cascade 
Y 4 = exit conditions of cascade 


INTRODUCTION 


D' RING RECENT YEARS, the art of axial compressor 
design has become increasingly more complex. 
In particular, there has been a steadily growing require- 
ment for moving larger masses of air at higher pressure 
ratios. 

With the large amount of excellent research and de- 
velopment work being conducted in the field of axial 
compressors by the N.A.C.A., the industry, the Mili- 
tary Services, and the universities, it becomes increas- 
ingly more difficult for research planners to obtain an 
overall picture of the field. With this in mind and 
from the point of view of the research planner, the 
present paper deals with a simplified qualitative com- 
parison of axial compressor types, with special emphasis 
on the relationship between mass flow and the stage- 
pressure ratio. 

It should be emphasized that, although many effects 
are neglected in the interest of obtaining rapid gener- 
alized comparison, continued and accelerated research 
on these effects is strongly urged. 

The present paper contributes to two projects in par- 
ticular. 

(1) Potential 
pulsion systems (project of the U.S.A.F. Scientific Ad- 
visory Board under Dr. Theodore von Karman). By 
increasing the mass flow through the compressor and 


improvements of gas-turbine-pro- 


the remaining internal system of a turbojet for the same 
frontal area the following improvements in the perform- 
ance of an aircraft or missile are made possible: (a) 
increase of maximum flight Mach Number and (b) in- 
crease of range at high Mach Number. 

(2) Potential improvements of flow generating sys- 
tems for test facilities (project of the Air Engineering 
Development Division under Major Gen. F. O. Car- 
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Fic. 1. Maximum possible mass-flow ratio versus flight Mach 
Number for a turbojet engine suspended in the free air stream. 
Isentropic flow in the intake diffusor is assumed. £& equals the 
ratio of maximum nacelle diameter to inner diameter of the com- 
pressor housing. 
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AXIAL MACH NUMBER 
Fic. 2. Influence of axial Mach Number and hub ratio on the 
mass-flow ratio of an axial compressor. 
roll). For test facilities of the magnitude and type of 


those under construction for the Arnold Engineering 
Development Center, the compressor or flow generating 
system is usually the largest single item of expense and 
complexity. Increasing the mass flow for the com- 
pressor system with no increase in external dimensions 
leads to the following advantages: (a) saving of ex- 
pense due to smaller size; (b) greater compactness 
(savings in ducts, building, and auxiliaries); and (c) 
utilization of higher speeds of rotation for electric 
motors. 

To explain these potential design improvements, the 
following comparative analysis of the relationship be- 
tween stage temperature rise and mass-flow ratio is 
developed for different compressor types, with simpli- 
fying assumptions. This is to facilitate the presenting 
of an overa'l physical picture of the problem and to spot- 
light areas worthy of further detailed study and ex- 
ploitation. 


HISTORICAL BACKGROUND 


Since the thrust of a turbojet is a function of the mass 
flow through the unit and the jet velocity, studies of 
potential improvements of a turbojet naturally involve 
study of possibilities of increasing the mass flow through 
the system and increasing the jet velocity—i.e., by 
utilizing higher operating temperatures. Studies of this 
nature were initiated in 1944 and 1945 by the A.A.F. 
Scientific Advisory Group.’ * A further study of this 
series was carried out in 1946 by Morton Alperin, who 
showed the potentialities of the turbojet for supersonic 
flight for increased values of thrust per unit frontal area. 

Over a period of years, the N.A.C.A. has been con- 
ducting extensive studies of potential turbojet im- 
provements. For instance, a report of Hensley, San- 
ders, and Kress* shows the importance of the mass-flow 
factor in improving turbojet performance, and Wood 
shows potentialities of the turbojet for supersonic flight. 

To the best of the author’s knowledge, the first de 
sign of a turbojet with special emphasis on the attain 
ment of high values of mass flow per unit frontal area 
was proposed in 1947 to the Air Materiel Command 
by Northrop Aircraft, Inc. In particular, Eugene 
Hunsaker presented original ideas on axial compressor 
design for unusually high mass flow. Discussions be- 
tween the author and Mr. Hunsaker were essential 


contributions to the present paper. 


DEFINITION OF Mass-FLow RATIO 


The factor used in this paper to represent mass-flow 
characteristics of axial-flow compressors is the ratio 
of the actual mass flow moved by the compressor to the 
maximum mass flow that could be passed at local sonic 
speed through the hollow compressor housing under 
unchanged isentropic stagnation conditions. This ratic 
is identical with the ratio of the critical area of the pass 
ing mass flow to the cross-sectional area of the hollow 
compressor housing. 

The relationship between maximum possible mass 
flow ratio and flight Mach Number is shown in Fig. 1. 
The diagram is based on the assumption that in the 
supersonic speed range the highest possible mass flow 
per second for an engine suspended in the free air stream 
cannot exceed the amount of air which would be passed 
with flight velocity through an area equal to the frontal 
area of the engine. 

The relationship between mass-flow ratio, axial Mach 
Number, and hub ratio is shown in Fig. 2... By defini 
tion, the combination of zero hub diameter and axial 
Mach Number of 1 corresponds to a mass-flow ratio 
of 1. The importance of having a high inlet axial Mach 
Number and a small hub diameter for obtaining high 
mass-flow ratios is shown by the curves. Two facts 
are especially significant--namely, (a) increasing the 
inlet Mach Number pays good dividends up to a value 
of about 0.7 to 0.8, above which the returns are dimin 
ishing, and (b) decreasing the hub ratio is advantageous 
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down to 0.3 and 0.4, below which returns are diminish- 
ing. 

It should be pointed out that the mass-flow ratio 
peaks at Mach Number 1, above which it decreases 
again. This shows that the utilization of supersonic 
axial velocities would have no advantage from a mass- 
flow basis, since it would decrease rather than increase 
the mass flow. Most axial compressors as used for 
many years in stationary gas turbine practices have 
hub ratios above 0.7 and axial flow below 0.517. In 
recent years, aircraft and wind-tunnel applications have 
brought an ever increasing need to increase the mass 
flow to a considerable extent. It, therefore, becomes 
important to study potential limits of increasing the 
mass flow. Referring again to Fig. 2, the question is: 
How near a mass flow ratio of unity is it feasible to ap- 
proach? This naturally leads to the question of the 
limitations on inlet, Mach Number and hub ratio in 
terms of the compressor blading design and the stage- 
pressure ratio or temperature rise. 

A series of blading design methods will be considered 
in this paper based on the limitations of high subsonic 
flow with several simplifying assumptions. 


TEMPERATURE RISE OF A BLADE ELEMENT 


The approach of this paper will be, in general, to 
study the relationships between stage-pressure ratio 
and mass flow for a number of different assumptions in 
order to give an overall physical picture of potential 
compressor performance, especially from the mass-flow 
point of view. 

Before starting this more general study, it is helpful 
to isolate one blade element and examine the influence 
of axial velocity, which is characteristic for mass flow 
and stage temperature rise. Fig. 3 shows a simplified 
velocity diagram for the maximum temperature rise 
for a blade element, assuming subsonic flow.® 

As shown in reference 5, the maximum temperature 
rise of a blade element is obtained when the relative 
velocities of the air entering the rotor and the down- 
stream stator have the highest permissible values under 
the Mach Number limitations. With both of these 
relative velocities limited by the approach to the speed 
of sound, a symmetrical velocity diagram is obtained. 
The temperature rise is proportional to the product of 
rotor speed and flow deflection through the rotor. Since 
the base of the velocity triangle is equal to the sum of 
rotor speed and flow deflection, the maximum value of 
the product is obtained when both legs are equal. 

The following simplifying assumptions are made: (a) 
constant axial velocity through the stage at the location 
of the element; (b) constant value of the velocity of 
sound through the stage, which appears justified since 
the stage temperature rise in a subsonic compressor 1s 
small relative to the inlet temperature; and (c) limiting 
Mach Number equal for rotor and stator. 
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Fic. 3. Theoretical maximum temperature rise for a blade 
element as a function of relative intake Mach Number and axial 
Mach Number. 
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Free vortex with positive prerotation and axial relative 
discharge velocity at hub. 


Fig. 3 shows a typical plot of the relationship between 
axial Mach Number and blade element temperature 
rise for several assumed values of the limiting relative 
Mach Number. It is seen that the temperature rise 
goes in the opposite direction to the inlet axial Mach 
Number and decreases rapidly at high Mach Num- 
bers; therefore, it is important for any specific design 
to determine the best compromise between these two 


factors. 

With this introductory view of the behavior of a 
blade element, the rotor as a whole will now be con- 
sidered, rather than one element. This results in an 


additional variable—namely, the hub ratio a. 
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discharge velocity over blade length. 


INVESTIGATION OF COMPRESSOR TYPES 


(1) Free-Vortex Flow Principle 


The so-called free-vortex method of design has been 
well known for many years. It is called free vortex 
because the tangential component of absolute velocity 
varies inversely as the radius, as it does in a free vortex. 
This is the easiest type to work with theoretically, since 
it involves a constant value of work or temperature rise 
over the whole area, constant axial velocity, balanced 
pressure, and no radial flow. For this reason it has 
been widely used in practice for many years. 
into use in the 1920's when such scientists as von K4r- 
man, Prandtl, Betz, Ackeret, and Troller made popular 
the application of airfoil theory to rotating-blade sys- 
tems. 

(a) Positive Prerotation, Axial Relative Exit Velocity. 

Fig. 4 shows typical velocity diagrams for free-vortex 
The tangential com- 


It came 


flow with positive prerotation. 
ponent developed by the rotor varies inversely as the 
radius. Study of the velocity diagrams shows that at 
the tip the relative intake velocity is critical for limiting 
Mach Number and that at the hub the absolute dis- 
charge velocity is limiting. Now, since in the present 
study increasing the mass flow is an important objec- 
tive, low values of the hub ratio are of particular in- 
terest. As the hub ratio decreases, the relative exit 
velocity approaches the axial direction at the hub and 
can even exceed it. However, for the first investiga- 
tion, the simplifying assumption of axial relative exit 
velocity at the hub will be introduced. This assump- 
tion simplifies the calculations and is in approximate 
agreement with many existing compressors of this 
type. With the above assumption and the Mach Num- 
ber limitation, the following formula can be derived 
(as shown in the Appendix) : 


Free vortex with negative prerotation and axial absolute 


Jg¢,AT = UsUR = 
[1 + (1/a) — (1/a?)]c?(M? — M,,°) 


giving the relationship between temperature rise, hub 
ratio, and relative and axial Mach Numbers. In order 
to isolate further the hub ratio effect, a temperature 
coefficient (Cr) is introduced by dividing the temper- 
ature rise by the difference in the squares of the relative 
and axial Mach Numbers—1.e., 

C Jee, AT l l 

i) es ae 

This formula immediately shows the inherent disad- 
vantage of this type of compressor for high mass flow, 
since the temperature coefficient drops to zero for a hub 
ratio somewhat above 0.6. In practice, vortex com- 
pressors of this type usually have hub ratios of 0.7 to 
0.8 and axial Mach Numbers below 0.5. Referring to 
Fig. 2, it is seen that the mass-flow ratio for this type of 
compressor would amount to less than one-third of the 
ceiling value. With this in mind, other design methods 
will be discussed from the point of view of essential in- 
crease in the mass-flow ratio. 

(b) Free-Vortex Flow—Negative Prerotation, Axial 
Absolute Discharge Velocity.—This example will show 
typical limiting considerations for an inlet compressor 
stage’ with free-vortex flow combined with negative 
i.e., prerotation opposite to the direction 
This method has been widely used 


prerotation 
of rotor rotation. 
in axial ventilating fan practice, since it utilizes lower 
rotational speeds than positive prerotation and is there- 
fore efficacious when tip speeds are limited. Positive 
prerotation is more widely used in gas-turbine practice 
where higher rotational speeds more nearly match those 
with high-speed turbines. Velocity diagrams for this 
type of flow are shown in Fig. 5, which shows the com- 
monly used method whereby the absolute velocity 
leaving the rotor is axial. In other words, the prerota- 
tion vanes introduce against the rotor the same amount 
of spin which the rotor takes out. Using simplifying 
assumptions as in the previous case, the formula for the 
temperature coefficient is Cr = a/(1 + a)*. It is seen 
from this formula that the cutoff of temperature rise 
with decreasing hub ratio has been eliminated. 

(c) Free-Vortex Flow—Negative Prerotation, Axtal 
Relative Discharge Velocity at Hub.—When the limita- 
tion of axial absolute velocity behind the rotor is re- 
moved and replaced by axial relative discharge at the 
rotor hub to correspond with the previously considered 
case of free-vortex flow with positive prerotation, the 
temperature coefficient is Cr = a/(a + 1). It is seen 
that these two types with negative prerotation have 
smaller temperature rise or rotor work for the higher 
range of hub ratios than the previous positive prerota- 
tion type (see Fig. 10). However, this design method 
does open up a lower range of hub ratios making pos- 
sible higher mass-flow capabilities. The disadvantage 
of negative prerotation compressors for turbojet appli- 
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cation is the lower rotor speeds involved. This method 
was used, however, in the compressor design of the 
Jumo 004, where low rotor speeds were also necessi- 
tated by the use of relatively low stress materials and 
simplified construction methods suitable for mass pro- 
duction at that time. 

(d) Free-Vortex Flow Rotor Work. 
In the previous investigations of the free-vortex type, 
there was an arbitrary limitation of axial discharge 
velocity at the hub. Now it is of interest to investigate 
the potentialities of removing this limitation whereby 
the theoretical maximum rotor work is obtained. The 


Maximum 


analysis of this case is contained in the Appendix, and 
the temperature coefficient is Cp = (1 + a)*/4 (see 
Fig. 10). 

The theoretical maximum rotor work for free-vortex 
flow involves negative prerotation as shown in Fig. 6. 
However, for a practical range of hub ratios, the increase 
over the zero prerotation type is small. 

This maximum type is given to show potential ceil- 
ings of the free-vortex principle. In actual practice, 
this maximum rotor work would involve considerably 
higher blade loading than is commonly used. Blade- 
loading limitations will be discussed later in this report. 

A rotor of this type was designed in 1943 by the 
author and Elmer Johnson for wind-tunnel application. 
A scale model was built at Wright-Patterson Air Force 
Base for the Wind Tunnel Branch, Air Materiel Com- 
mand. 

It is noted that the blade element at the hub ap- 
proaches the impulse type. Interesting investigations 
on the impulse type have been made by Erwin and 
Schulze, of the N.A.C.A. 


(2) Positive Rigid Body Prerotation Constant Rotor Work 


In this type of flow, the upstream stator introduces a 
rigid body type of swirl into the air stream-—namely, 
with the tangential component of absolute intake ve- 
locity directly proportional to the radial distance to 
the center of rotation (Fig. 7). For this type there is 
a deviation from the constant axial velocity distribu- 
tion which is characteristic for all constant rotor work 
types not designed according to the free-vortex prin- 
ciple. This influence is not included in the present 
report but is being studied separately. The subject 
has been treated by Rannie, Wu and Wolfenstein,’ and 
The analysis has been carried out in the 
namely, for axial rela- 


others.’ 
Appendix for two assumptions 
tive discharge at the hub for which Cr = a, and for 
maximum constant rotor work for which Cy, = (1 + 
a)?/4. It is seen that for the axial relative discharge 
velocity limitation the rigid body design allows lower 
ranges of hub ratios to be utilized than does the corre- 
sponding free-vortex design with positive prerotation 
(Fig. 10). It is interesting to note, however, that the 
rigid body maximum constant rotor work type has 
exactly the same expression for temperature coefficient 


as the maximum free-vortex type. However, there 
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Fic. 6. Free-vortex maximum constant rotor work. 
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Fic. 7. Rigid body prerotation constant rotor work. 


is a significant difference between the two types in that 
the tangential velocity imparted by the rotor is smaller 
for the rigid body type. Further comparisons of the 
two types will be made later in this report. 


(3) Maximum Constant Rotor Work 


In the cases so far considered, the conditions of maxi- 
mum rotor work were determined (1) for free-vortex 
prerotation and (2) for rigid body prerotation. Now it 
becomes of interest to investigate the more general con- 
dition of maximum constant rotor work by removing 
any specification as to the type of prerotation. The 
analysis given in the Appendix leads to a temperature 
coefficient of C, = 4a/(1 + a)? (see Figs. 8 and 10). 
For this case the prerotation varies from positive to 
negative over the blade, as shown in Fig. 8. _ If the pre- 
rotation is restricted either to positive values over the 
blade length or to negative values over the blade length, 
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Fic. 9. Maximum unconstant rotor work. 


the expression for temperature rise coefficient in both 
cases is Cr = a(2 — a) (see Fig. 10, Curve 6). 


(4) Nonconstant Distribution of Rotor Work over the 

Radius 

In all cases investigated so far, the condition of con- 
stant rotor work was assumed throughout. This con- 
dition obviously limits the total rotor work because the 
work goes to zero as the hub ratio approaches zero. 
Therefore, an increase of total rotor work can be ex- 
pected for compressor types that have the rotor work 
increasing toward the blade tip. Two cases will now 
be considered. 

(a) Constant Rotor Swirl, Maximum Rotor Work for 
a Blade Element at the Tip. 
Appendix shows the temperature-rise coefficient to be 
Cr, = r/R for a blade element at radius ry. For any 
given hub ratio, an average C; value can be obtained by 


The analysis given in the 


integration as follows: 
(Cr)av, = (2/3) [1 — a*)/(1 — a?)] 


see Fig. 10, Curve 7. This case illustrates a possible 
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way of combining high mass flow with reasonable stage 
pressure or temperature rise. 

(b) Maximum Nonconstant Rotor Work.—For this 
case all restrictions as to prerotation and work dis- 
tribution over the blade length are removed and only 
the Mach Number limitation remains. 

The analysis of this case is made in the Appendix and 
the resulting value of temperature coefficient is 

8 (1 — a)? 
(Cr)av. = 9 (1 — a*) (1 — a’) 
see Figs. 9 and 10, Curve 9. 

It can be seen that the value of the hub ratio has little 
influence on the average temperature coefficient. 

(c) Inducer Stages.—An inducer stage is a prelimi- 
nary rotating stage with unconstant rotor work distribu- 
tion over the blade length. Its primary purpose is to 
achieve a desirable axial velocity distribution in a multi- 
stage compressor which deviates from the free vortex 
design principle.’ Since the unconstant rotor work 
types are capable of high temperature rise along with 
high mass flow, it appears attractive to combine this 
characteristic with the primary function of producing 
a desirable axial velocity distribution. Studies in this 
direction are being made, and further work should be 


encouraged. 


(5) Comparison of the Investigated Compressor Types 


(a) The Influence of Hub Ratio upon the Tempera- 
ture Rise.—The previously derived temperature-rise 
coefficients for the investigated compressor types and, 
in addition, the temperature-rise coefficient of a hypo- 
thetical type, which is characterized by the absence of 
the hub ratio effect, are plotted in Fig. 10. In compari- 
son to this hypothetical type, the curves of the investi- 
gated compressor types show a considerable variance 
as to the influence of hub ratio upon the temperature- 
rise coefficient. The rigid body types and the free 
vortex types are greatly influenced by the hub ratio, 
while the maximum rotor work types are little affected 
by the hub ratio within the practical range of hub ratios 
of 0.8 to 0.9. The greatest difference in temperature 
rise among the constant rotor work compressor types is 
shown by the maximum constant rotor work type and 
the free-vortex type with positive prerotation and axial 
While the tem- 


perature rise of the latter type decreases to zero when the 


relative discharge velocity at the hub. 


hub ratio decreases to 0.62, the temperature rise co 
efficient of the first type is still above 0.9. 

Another interesting result is that the temperature- 
rise coefficients for the maximum free-vortex and the 
maximum rigid body constant rotor work types are the 
same. Also, the temperature rise coeflicients for the 
maximum constant rotor work positive prerotation 
type and the maximum constant rotor work negative 
prerotation type are the same. This shows that the 
larger flow deflection of the negative prerotation types 
is compensated by the higher circumferential speed 
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RELATIONSHIP BETWEEN TEMPERATURE RISE COEPFICIENT AND HUB 
RATIO FOR THE COMPRESSOR TYPES LISTED BELOW: 
1, PREE VORTEX-POSITIVE PREROTATION-AXIAL RELATIVE 
DISCHARGE VELOCITY AT HUB 
2. FREE VORTEX-NEGATIVE PREROTATION-AXIAL ABSOLUTE 
DISCHARGE VELOCITY OVER RADIUS 
3. FREE VORTEX-NEGATIVE PREROTATION-AXIAL RELATIVE 
DISCHARGE VELOCITY AT HUB 
4. RIGID BODY PREROTATION-CONSTANT ROTOR WORK-AXIAL 
RELATIVE DISCHARGE VELOCITY AT HUB 
5. FREE VORTEX-MAXIMUM ROTOR WORK 
RIGID BODY PREROTATION-MAXIMUM CONSTANT ROTOR WORK 
6. MAXIMUM CONSTANT ROTOR WORK-POSITIVE PREROTATION 
MAXIMUM CONSTANT ROTOR WORK-NEGATIVE PREROTATION 
7. CONSTANT ROTOR SWIRL-MAXIMUM TEMPERATURE RISE OF 
BLADE ELEMENT AT TIP 
8. MAXIMUM CONSTANT ROTOR WORK 
9. MAXIMUM NON-CONBTANT ROTOR WORK 
10. HYPOTHETICAL (NO HUB RATIO EFFECT) 
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TEMPERATURE RISE COEFFICIENT (Cy) 
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HUB RATIO 
Fic. 10. 


of the positive prerotation types. The problems in- 
volved with the blade loading, which is much higher 
for those investigated types with negative prerotation 
than for those with positive prerotation, will be discussed 
later in this paper. 

(b) Temperature-Rise Ratio Versus Mass-Flow Ra- 
tio.—Having isolated the hub ratio effect, it is now of 
interest to combine this effect with the axial Mach 
Number effect to obtain the relationship between mass- 
flow ratio and temperature-rise ratio. This is shown 
in Fig. 11. These curves are derived by using Figs. 2 
and 10 and are based on a limiting relative Mach Num- 
ber of 0.8. For other limiting Mach Number values 
(.\/’), the curves in Fig. 11 can be used by multiplying 
the temperature-rise ratio by the factor (/’/0.8)? and 
by correcting the mass-flow ratio by using Fig. 2. Fig. 
2 gives for any mass-flow ratio a manifold of hub ratio- 
axial Mach Number combinations. The temperature- 
rise ratio A7/T, for each combination is given by the 
equation 

Af = ——_ (Mt Me) (Cr) 

Ty J8lyTo 
However, for one given mass-flow ratio there is only one 
combination of hub ratio and axial Mach Number which 
gives a maximum value of AT 7, for each compressor 
type. These maximum values make up the curves of 
Fig.11. The curves show that with an increasing mass- 
flow ratio the temperature-rise ratio decreases rapidly 


RELATIONSHIP BASED ON MACH NUMBER LIMITATION BETWEEN MASS 
FLOW RATIO AND MAXIMUM IDEAL TEMPERATURE RISE RATIO FOR 
THE COMPRESSOR TYPES LISTED BELOW: 


1. FREE VORTEX-NEGATIVE PREROTATION-AXIAL ABSO- 
LUTE DISCHARGE VELOCITY OVER RADIUS 
2.5 2. FREE VORTEX-POSITIVE PREROTATION-AXIAL RELA- 
TIVE DISCHARGE VELOCITY AT HUB 


P 3. FREE VORTEX-MAXIMUM ROTOR WORK 
P, RIGID BODY PREROTATION-MAXIMUM CONSTANT ROTOR 
lar WORK 
= 4. MAXIMUM CONSTANT ROTOR WORK 
° 5. MAXIMUM NON-CONSTANT ROTOR WORK 


6. HYPOTHETICAL (NO HUB RATIO EPPECT) 











0.2 0.4 0.6 0.8 1.0 


MASS FLOW RATIO 
Fic. 11 


to zero for all compressor types, including the hypo- 
thetical type where the hub ratio has no influence on 
the temperature rise. Furthermore, it is apparent for 
higher mass-flow ratios, in a range of about 0.5 to 0.7, 
the maximum temperature-rise ratio is greatly influ- 
enced by the compressor type. 


(6) Investigation Based on Blade Loading Limitation and 
Mach Number Limitation 


In the previous investigation of several compressor 
types, the primary limitation was that of the intake 
Mach Number. This corresponds to the limiting case 
of infinite number of blades of infinitesimal thickness, 
where the intake Mach Number is the highest occurring 
Mach Number value and blade-loading effects do not 
exist. However, for a finite number of blades, the 
limiting intake Mach Number cannot be considered 
as independent of the compressor type, because the 
higher the blade loading which is imposed by the type 
of compressor, the higher will be the overvelocity at the 
suction side of the blades. Therefore, the higher the 
blade loading, the smaller will be the limiting Mach 
Number. 

In addition, the blade loading itself is limited by the 
tendency of boundary-layer separation if the pressure 
increase in the direction of flow becomes too great. 

For cascades with conventional airfoil blades, Keller'*® 
derived a limitation for the amount of deflection which 
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can be imparted to the air stream because of mutual 
blade interference. 

If C, is the blade lift coefficient and a is the ratio of 
blade chord to spacing, then the approximate limitation 


Keller is 


2Ur/ Wn 


proposed by 
= o(C,) = 1.2 


see Fig.12). This limiting value of flow deflection will 
g g 


AT 


Curves are plotted in Fig. 14 giving the relationship 
between maximum temperature-rise ratio and mass- 
flow ratio for the above example, for the free vortex 
positive prerotation axial relative exit velocity, for the 
maximum free vortex, for the maximum rigid body con- 
stant rotor work, and for the maximum nonconstant 
rotor work. However, these curves represent only 
relative values, since the limiting value o(C,) = 1.2 
is by no means a universal constant. Its value is in- 
fluenced by Reynolds Number, Mach Number, cas- 
cade configuration, etc., which have not been suffi- 
ciently investigated at this time by experimental 
work." 


7) Comparison of Compressor Types with Blade-Loading 

Limitation 

Fig. 14 shows that the type of compressor has great 
influence upon the temperature-rise ratio. This is true 
especially in the higher range of mass-flow ratios. In 
contrast to the results obtained with the Mach Num- 
ber limitation only, the rigid body type with positive 
prerotation enables higher values of temperature rise 
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be basic for the analysis. In addition, the same limiting 
Mach Number of 0.8 will be applied for comparison 
purpose. 

An analysis is given in the Appendix for the maximum 
constant rotor work type as an example. The velocity 
diagram is shown in Fig. 13. The analysis results in 
the following simplified expression for the temperature- 
rise ratio of the maximum constant rotor work type 


= Gi) ome (SK 


than the free-vortex type, except for the small range of 
extremely high mass-flow ratios where the blade loading 
limitation does not come into effect. The nonconstant 
rotor work type gives reasonable values of stage tem- 
perature rise up to extremely high mass-flow ratio 
values of about 0.8. 

Although it would be desirable to compare efficiencies 
of the various types, this is not done here, since an effli- 
ciency investigation based on two-dimensional flow 
theory is insufficient. The important three-dimensional 
effects, such as boundary layer, centrifugal effects, 
induced drag, gap influences, and influence of surface 
friction, are not sufficiently known* to be 


analytical investigation. 


used in an 


CONCLUSIONS 


Simplified analyses such as presented above help 
research planning and integration and spotlight areas 
worthy of further investigation. 

The following are typical research items brought out 
by such simplified analyses: 


(1) 


(2) Studying limits of loading and deflection possi- 


Increasing relative Mach Number. 


bilities. 

(3) More complete comparative analyses with more 
items considered. 

(4) Radial flow and equilibrium. 


(5) Losses—especially mixing with work unbal- 


ances~ secondary flows. 

(6) Reynolds Number effects (high-altitude prob- 
lems). 

(7) Blade thickness effects on velocity limits. 

(S) Inducer stages. 

(9) Staging and matching. 

Much excellent research is in progress in the above 
fields. It is hoped that simplified analyses like the 
above will give planners greater appreciation of the 
need for expediting and supporting present programs 


and for stimulating new programs 
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APPENDIX 


1) Free Vortex 


(a) Positive Prerotation, Axial Relative Discharge 
Velocity at Hub.—Study of the velocity diagrams (Fig. 
4) shows that the limiting Mach Number (J/) is attain- 
able only for the relative intake velocity at the tip and 
the absolute discharge velocity at the hub. This is be- 
cause the relative intake velocity decreases toward the 
hub, while the absolute discharge velocity increases 
toward the hub. 

As the relative intake velocity, Wj, is limited at the 
tip section of the rotor by the limiting Mach Number 
(\/), the magnitude of the tangential component (W;)z 
is fixed for a given value of axial velocity (V’)—.e., 


(Wise = ¢V M? — Miz? 


lf C\, represents the whirl component introduced by 
the upstream vanes, then the blade element velocity is 
given by the sum of W,, and C), as 
U, = (Wiu)r + (Ciu)> 
At the blade tip, where r = R, 
Un = (Wire + (Cie 
The blade element velocity at the hub is 
Uar = a(UpR) 
For free vortex flow the following relationships hold: 


(Cin)ar = (1 a) (Ciu)r 


and 
(Urdar = (1/a) (Ur) e 
Since 
(Co)ar = (Wi)e = cM 
then 
(Cowar = (Wiwe 
so that 


(Ciwar + (Urdar = (Wide 


For vortex flow, in which case (C2,), varies inversely 
with the radius, this can be written 


[(Ciu)e a| — [(Ur)e a| = (Wiu)e 


From the assumption of axial relative discharge ve- 
locity at the hub, it follows that 


aUp = |(Ur)e al) + (Cwr/a)] = (Wide 
Multiplying by a the equation 
((Ur)e/a] + [(Cw)e/a] = (Wide 
and adding 
Ur — (Cider = (Wide 


results in 
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Fic. 13. Maximum constant rotor work with blade-loading 


limitation applied. 





RELATIONSHIP BASED ON BLADE LOADING AND MACH NUMBER 
LIMITATION BETWEEN MASS FLOW RATIO AND MAXIMUM IDEAL 
TEMPERATURE RISE RATIO FOR THE COMPRESSOR TYPES LIST- 
ED BELOW: 


« FREE VORTEX-MAXIMUM ROTOR WORK 

« RIGID BODY PREROTATION-MAXIMUM CONSTANT 
ROTOR WORK 

MAXIMUM CONSTANT ROTOR WORK 

MAXIMUM NON-CONSTANT ROTOR WORK 

FREE VORTEX-POSITIVE PREROTATION-AXIAL 
2. RELATIVE DISCHARGE VELOCITY AT HUB 











0.2 0.4 0.6 0.8 1.0 
MASS FLOW RATIO 
Fic. 14. 


(Ur)re + Urn = (1 4+ @) (Wiude 
Substituting for LU’, the value (W1,,)e/a gives 
(Ur)r = [1 +a a sen (1 a) (Win)e] 


The stage temperature rise is then given by the expres- 
sion 
Jgcy AT = UpUr = 

{1 + (1/a) — (1/a*)] [ce?(? — M,,")| 


For vortex flow this value is constant over the radius. 
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Since the formula given above for temperature rise is 
a product of a Mach Number effect and a hub ratio 
effect, comparisons may be further simplified by intro- 
ducing the following coefficient that isolates the two 


effects—namely, 


: J gc, AT ' l l 
—s c?(M? — M,,) "9 a @ 

(b) Negative Prerotation, Axial Absolute Discharge 
Velocity over Blade-—A velocity diagram for such a 
blade element is shown in Fig. 5. 

As the absolute discharge velocity is assumed axial 
over the blade length, 


(Ur)r + Ur = (Wiwe 


so that there is no right-hand velocity triangle as in the 
previous case. In order to obtain the highest rotor 
work, (W),), must be the largest possible under the 
Mach Number limitation. It follows from this and 
from the condition of constant rotor work that (ly), 
is smaller than Otherwise, the Mach Number 
limitation would be violated immediately for smaller 
radii. Since (W,,), decreases with the radius until 
(Up), equals U,, at which point (W,,), begins to in- 
crease, the Mach Number limitation will not be violated 
if, at the hub 


Up. 


(Uor ai (Ur)ar < (Wie 


or 
aUp + [(Ur)e/a] S (Wide 
Since 
(Ur)r + Un = (Wiu)e 
aUp + [(Ur)r/a] S Urn + (Ure 
or 


(Ur)r < aur 


For maximum rotor work, (U;)z must have the largest 
permissible value, so that 


aUr, Ur(l + a) = (Wide 


(Ur)r - 


The stage temperature rise can now be expressed as 


Jgc, AT = (Upp (Up) = 


Study of this equation shows that the highest stage 
temperature rise results when k is zero or prerotation 
zero. The equation for highest rotor work is then 

JgcyAT = ac? (M? — M,,;?) 
The temperature rise coefficient is 


Cr = JgcpAT/c?(M* — Maz”) = a 
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c2(M? — M, (a) + 
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J gc, AT = (Up)reUr = la (a + 1)*]c?(M? —_ V/ - 
The expression for temperature-rise coefficient is 
Cr = a(a + 1)? 


(d) Maximum Rotor Work. 
tion of the case of positive prerotation, it becomes eyj- 


First, in the investiga- 


dent that the relative intake velocity decreases toward 
the hub, since U, is decreasing and (Cj,), is increasing, 
Also, the absolute discharge velocity increases toward 
the hub since both (C,,), and (Ul), are increasing, 
This means that, in order to attain the highest possible 
rotor work without violating the Mach Number con- 
dition, only the relative intake velocity at the blade tip 
and the absolute discharge velocity at the hub can 


attain the highest permissible Mach Number. Thus, 


(u wR = (Cowar 
(Cou)ar = (Ur)ar 7. (Ciwar _ (Wiwe 
Since (C\,), and (U;), vary inversely as the radius, 
one can be expressed as a function of the other as fol- 


lows: 


(Ci), = R(U er), 


where k must be positive for positive prerotation. Sub- 
stituting k(Ur)ar for (Ciu)ar in the above equation re- 
sults in 


(U.) (Wie ww A(Wiw)R 
de = or j ») = 
F)aR i+ FJR rr 
because of constant rotor work. 
Since 
Ur _ (Wiwe + (Ciu)r 
then 
Ur = (Wiwe + R(Ur)e 
so that 


a 
1+k 


This gives the stage temperature rise as a function of k 


Ur = (Wwe (1 + 
and hub ratio a as follows: 


" Fea we | 
2— M,,”) a ; 
1+ k(1 — 2k) 


Investigating the case of negative prerotation, we 
find, as in the previous case of positive prerotation, that 
the limiting Mach Number can be achieved only for 
the relative intake velocity at the tip and for the ab- 
solute discharge velocity at the hub. It will be shown 
that, for hub ratios zero up to a definite value (0.41), 
the Mach Number limitation also comes into effect for 
the relative intake velocity at the hub. 
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DESIGN COMPARISONS 


Applying the same method as for the case of positive 
prerotation, the following equation is obtained for the 


temperature rise: 


U; l - = 


a ak 
(M2 — M,,?) =) .) 
(, +e 14h 


where & is negative for negative prerotation. The 


Jgcp AT - 


value of this expression is a maximum for 


k = ~-|[(1 —a)/(1+a)] 
Substitution of this value k& into the above expression 
gives 
wo, ED is ar 
Jgc,AT = ' (M? — her’) (c*) 
The temperature rise coefficient is C7; = (1 + a* 4), see 


Fig. 10, Curve 5. As previously mentioned, this ex- 
pression is limited to values of hub ratio a greater than 
or equal to 0.41. This can be shown as follows: 


The Mach Number limitation will not be violated if 


(Wiwer < (Wiwe 


which means 
(Winder sai al, + (1 a) (Ciu)r < (Wie 


Substituting the functions of (W,)e for Up and (Cider 
results in the equation 


1 <a}l — [k/(1 + R)] [(1/a) — a}} 


and replacing k with —[(1 — a)/(1 + a)] gives a > 
V2 — 1 = 0.41. For hub ratios smaller than 0.41, 
the expression for the temperature-rise coefficient is 


Cr = a(2 + a)/(1 + a)? 


2) Positive Rigid Body Prerotation—Constant Rotor 
Work 


(a) Axial Relative Discharge at Hub.—In this case 
the tangential component of the relative intake velocity 
decreases proportionally with the radius because of 
rigid body prerotation. Therefore, the Mach Number 
limitation will be applied at the tip for the relative in- 
take velocity. For the absolute discharge velocities, 
the Mach Number limitation will not be violated if 
Urn + (Ur)r < 2(Wiw de (1) 
aur S (Wwe (2) 


Eq. (2) follows from the assumption of axial relative 
discharge at the hub for which 
(Uvp)er = (Wiu)er or (Ur) pr a= a(Wiwr 

so that Eq. (1) now becomes 


Ur < (2 — a?)(Wiu)e 


Considering the case 
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Ur 4 (2 - a*)(Wria)r (la) 


which means that the absolute discharge velocity at the 
tip is equal to the limiting Mach Number, and Eq. (2), 
then it follows that the hub ratio a is limited to values 
less than or equal to 0.63, approximately. Next, con- 
sidering Eq. (2) as an equality only, which means the 
absolute discharge velocity at the hub is equal to the 
limiting Mach Number, and the Eq. (la), then it fol- 
lows that the hub ratio a is limited to values greater 
than or equal to 0.63, approximately. 

The expression for stage temperature rise for the 


first case, a < 0.63, is 
Jgc,AT = a?(2 — a*) (M,? — Mez?) (c?) 


so that C; = a?(2 — a’). 
The expression for stage temperature rise for the 


second case, a > 0.63, 1s 
JgcyAT = a(M? — M,,”) 


so that Cp = a (see Fig. 10, Curve 4). 
(b) Maximum Constant Rotor Work. 
rigid body prerotation, the relative intake velocity de- 
creases toward the hub while the absolute discharge ve- 
As in the previous case 


For positive 


locity decreases toward the tip. 
of rigid body prerotation, only the relative intake ve- 
locity at the tip and the absolute discharge velocity at 
the hub can attain the limiting Mach Number. How- 
ever, for hub ratios smaller than 0.41, the Mach Num- 
ber limitation comes into effect for the absolute dis- 
charge velocity at the tip. The analysis is as follows: 


(Wie = (Coular Or (Wide = (Co)ar 


From this it follows that 
Urn — (Cure = (Ciwar + (Ur)ar 
or 
Ur — (Cure = aA(Ci)re + 
(1/a) (Ur)e = CV M? — Maz? 


Now, since (C,,), and U, vary directly as the radius, 
(C,,,), can be expressed as a constant k times L’,. There- 


fore, 
Uri — k — ak)a = (Ur)e 
Up = c[V M? — M,;2/(1 — )| 


This gives the stage temperature rise as a function of k 

and the hub ratio as follows: 

afl— Rk + a))) 
(l1—k)? f 


Url(Ur)e] = yur — M..? 


The value of this expression is a maximum when k = 
(1 — a)/(1 + a). Substitution of this value of & into 


the above expression gives 
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(1 + a)? 
¢ 


/ ge, AT = r2( MM? = M a2”) 


so that Cr = (1 + a)?/4. 


(3) Maximum Constant Rotor Work 


For maximum constant rotor work, (W1,,)2 and (Cox) 
must be the largest possible under the Mach Number 
limitation. From this and from the condition of con- 
stant rotor work, (U,)p must be smaller than Up, 
otherwise the Mach Number limitation would be vio- 
lated for all other radii. Thus, the sum of (IW,,), and 
(C2,), decreases with the radius until (l/y), equals U, 
at which point (IWV,,), begins to increase. However, 
the Mach Number limitation will not be violated at 


any point if at the hub 


Var = (Up)ar < 2 (Wiwe 


or 
aUe + ((Ur)e/a] $2 (Wwe 
Since 
(Ur)e + Un = 2 (Wie 
then 
aun + [((Ur e/a] < (Ure + Ur 
or 


(Ure < aU pe 


For maximum constant rotor work, (l’,), must have 
the largest permissible value, so that (Up)p = aUp. 
Therefore, 


Ur + ale = (Wide 
or 
Up (1 + a) = 2c VM? — M,:? 


The stage temperature rise is then given by the equa- 


tion 


fa 


Jt, 41 = aU,* = 
& —e + a) 


- (M? — M,,*) (c’) 
or C, = 4a/(1 + a)?. 


(4) Nonconstant Rotor Work Distribution over the Radius 


The equation of temperature rise for nonconstant 


work distribution is given as 


9 "RR 
Jyc, AT = — / r(Up)(U,) dr 
: R11 = a”) aR , 


(a) Constant Rotor Swirl—Maximum Rotor Work 
at T11p.—-From the assumption of maximum rotor work 
at the tip, it follows that 


(Ur)r = Ur = (Wie = cv M? — M,.* 
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From the condition of constant rotor swirl, (U,), = 


(Ur)r- 
Since 


(Ur), + Urn = [1 + (7/R)] (Wide 
no Mach Number violation will occur at any radius if 
the prerotation is chosen properly. 
For any blade radius, 


Ce il , r r r » 
J gc, AT = Ur (l FrR= R (U lu)R™ 


R 
The expression for the temperature rise coefficient is 
C J gc, AT r 
"(M2 — Ma?) R 


For any given hub ratio an average value can be ob- 
tained by integration as follows: 


oll 3 2 (1 — a?) 


l } 
(Cr) wv. = —, : f 2ar dr = 
'  @R*(1 — a®) JakRR 


oU — @) 
see Fig. 10, Curve 7. 

(b) Maximum Nonconstant Rotor Work. 
to obtain the maximum total rotor work, the following 
conditions must be considered : 

(1) The Mach Number of relative intake and ab- 
solute discharge velocities over the blade length must 
attain the highest permissible values. 

(2) The theoretical maximum work of a blade ele- 
ment must occur at one point r’ between tip and hub. 
This means that Uy equals U for the radius 7’. It 
should be noted that r’ cannot be expected to coincide 
with blade midpoint, since the mass flow of the annulus 
between blade midpoint and hub is smaller than be- 
tween blade midpoint and tip. 

The ratio r’/R may be designated as k. 
perature rise as a function of k can be obtained as 


In order 


The tem- 


follows: 
U, = (Wie (“) = (Wi)r’ — 
7 k(R) 
and 
(Ur) = 2(Was) — (Wied a ‘ 


Yr 
(Wi) + (2 ~ ) 
R(k) 


Therefore 


Toc AT x 2 Wine? LS. r ( r ) r| 
JS 2Cy = r * ie r 
ms R21 — a?) J or’ R(R) R/k) ‘ 


— Ww)? ke —a‘*) (1 - =] 
~ (1 — a?) 3k 4k? 


The maximum temperature rise is obtained when 


~ 


k = 3(1 — a*)/4(1 — a) 


so that 
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ae ey Bt M2 — Muz?)c? 
Cy - Ma" —~— Mes" )C* 
18% 9 (1 — a’) (1 — a) \ : 


and 


(Cr)av. = (8/9) [A — a*)?/(1 — a*)2(1 4+ a?)] 


6) Maximum Constant Rotor Work with Keller 
Limitation 


The highest possible flow deflection Uy, according to 
the Keller limitation, is Uy = 0.6W,,. This limitation 
applied for rotor and stator at the hub results in 


(1/a) Up = 0.6W,, = 0.6 V (a?/4) U2 + V? 


This results in the following expression for LU’: 
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see Fig. 13. The Mach Number limitation comes into 


effect at the tip where 
cM? = V2 + (0.45U + 0.5U,r)? 


Substituting into this equation the value lL’, results in 


2y = V24 (U +06 fa*U? ") 
= | 0d V { 


(U + 0.6a V (a2U?/4) + V?)? 


cM 
| 


MM ,..*) = 


cV M? — M,,2 — 0.8U = 0.34 V¥ (a2U2/4) + V2 


2c V M? — M,," — 0.6a V V7(1 — 0.09a*) + a*c?(M? — M,,”) 


: l 


U = 2 VM - 


and 


. pa | (1 — a*) (Maz? 
Up = 0.6ac V My, + a*(M? — Maz”) = 0.6a*Me 4/1 + . ) 


This results in the following expression for A775: 


(t= <) (ite 
a” M? 


AT ~ 94 Uiue) 


r, . =(0.24a7M-* \ 1+ 
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The Turbulent Boundary Layer in a 
Compressible Fluid with Positive 
Pressure Gradient’ 


CARLO FERRARIF 


Cornell Aeronautical Laboratory, Inc. 


SUMMARY 


The turbulent boundary layer in compressible flow is investi 
gated for the case where there is a pressure increase in the direc- 
tion of motion in order to determine its effect upon: (a) the ve- 
locity distribution at an arbitrary cross section of the layer; 
(b) the variation of the thickness of the layer as one progresses 
along the boundary and the calculation of the corresponding 
drag coefficient; and (c) the phenomenon of separation of the 
stream from ihe wetted wall. 

A check, resulting in close agreement, is made of the theory 
by comparing the velocity profiles predicted on this theoretical 
basis with those obtained in experiments at low speed under the 
influence of pressure gradients. A further numerical application 
is worked out for the case of a supersonic (Mach Number = 1.5) 
undisturbed stream velocity, but no experimental data at such 
high speed are at present available for confirmation. 


INTRODUCTION 


— INVESTIGATION OF THE turbulent boundary layer 
in supersonic flow was confined in a previous re- 
port! to the case of flow along a flat plate, but now this 
restrictive assumption is relinquished in order to pre- 
sent methods applicable to an airfoil—i.e., to account 
for the presence of a positive pressure gradient. 

A plausible assumption concerning the dependence 
of the turbulent friction stress on the kinematic and on 
the geometric characteristics of the average flow leads 
to the construction of the equation of motion, describ- 
ing the flow through the ‘‘inside’”’ portion of the layer, 
in the same form as that used by von Mises? for the 
case of the laminar régime. A procedure similar to 
that used by von Karman® for incompressible flow is 
employed in order to obtain the solution of this equa- 
tion of motion and for determining the position where 
the separation of the stream from the wetted wall even- 
tually takes place. 

It is further assumed that there is no heat transfer 
between the surface of the constraining wall and the 
fluid and that the flow is wholly turbulent. In addi- 
tion, the equations of motion are considered to be valid 
only at a sufficiently great distance from the place where 
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Torino. 


the increase in pressure in the direction of flow origi- 


nates. 
SYMBOLS 
R, = Ugx/vq = Reynolds Number of the undisturbed stream 
r i = absolute temperature 
p = pressure; r = tangential friction stress 
be = viscosity coefficient; » = kinematic viscosity 
p = density of the fluid; pe = density in the external 
stream at a cross section of the layer having an 
abscissa equal to x 
4 =l]-—- (p/pe); [= (zo — 2)/2 
U = component of the velocity along the x-axis; |’ = 


component of the velocity along the y-axis; Z = 
(0.-2/Ug?) — (U2/U¢2); OU. = velocity of the ex- 
ternal stream (its magnitude is taken to be equal 
to the component in the direction of the x-axis of 
the actual vector) at a location along the boundary 
having an abscissa equal to x 

W = VU? + V2; U* = friction velocity [U* = (rp/p))"); 
Vi = limit velocity of the undisturbed stream when 
expanded toa vacuum 

M = Mach Number; L = reference length; & = x/L; 


n = y/L; x being taken in the direction of the 
wetted surface 
6 = thickness of the boundary layer; § = 6/L 


. p U y 
6* = af l->—H: 
0 P, U. 5 
5* 5* 15 OU U g** 
; = —s sr maf 2 7 (1 - ay 6** = 
L 0 pe U, U. 6 i 


y = )Uy = stream function; ¥; = ¥/pa Val; ¥i* =a function 
that is proportional to y, and defined precisely in 
Section (B-1) 

* = functions of £ defined in Sections (B-1) and (C-1) 


fr 
fr 


Bars and primes will be used to denote, respectively, the mean 
values of the above quantities and their fluctuations about their 
mean values due to turbulence. The subscripts 0, p, e, and a 
correspond, respectively, to the stagnation condition of the un- 
disturbed stream and to other special conditions that exist at the 
wall, in the external stream, and in the undisturbed stream. 


(A) GENERAL EQUATIONS 


(A-1) Equations of Momentum, Energy, and Continuity 


Within the region external to the laminar sublayer 
and through neglect of the momentum transfer due to 
thermal molecular vibration in comparison to the 
momertum due to turbulence, one obtains 
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TURBULENT BOUNDARY 
° So* Gm UV 
Ox p= ox -? Oy \ ; 
: §<— Gh) +— ore we . 
— » = 2 = 
ae? Ox (p oy tel”) ' } 


It is being assumed here that the curvature of the sur- 
face of the obstacle along which the boundary layer is 
built up is sufficiently small to permit the supposition 
that the tangent to the surface, at every point P, is 


parallel to the x-axis for all intents and purposes. Be- 
cause 
pUr=G+teP)U+UYS 


U(pU + 2 p’U’) = pu? (2) 
inasmuch as it is generally true that pl’ > p’U’, 
while, also, 
pUV = pUV + U p'V' + BU'V" (2’) 
then one obtains 
op oO 


= oO =e 
pl?) (pUV) 
Ox = OX - sa ov ‘ T 


PS) - Ps 
— (BU'V’ + Up'V’) (3) 
oy 


Analogously and with the same approximation, one 
derives, from the equation of energy, 


Lost) bold) 
Uli + 5 dy L? i+3m) - 
Let H = 2 + (1/2) W” in the above expression and, 
upon noting that 
pUH = pUH + pU'H’ + op’ U' A+ oH’ U 
and 
epVH = pVA+ Ap’V’ +3, VV’ 


it results that 


0) se oO én js 
(BUH) + (jU'H' + p'U' H+ p'H’ U) + 
Ox Ox 
° GVA <= BV'H'’ + HT p’V’ 0 | 
) j= 
oy p oy \p p (4) 


Finally, from the equation of continuity, it is found 
that 
0 (307 oe - 0 —— a) —_— r 
Ox said oy al ld Ox vO) + oy eT ee 
It should be observed that, in the case of flow past a 
flat plate (i.e., with pressure gradient equal to zero at a 
sufficiently large distance from the transition point 
from laminar to turbulent flow), it has been found! ad- 
missible to write p’V’ = const. Bearing in mind that 
in the problem now under consideration the pressure 
gradient is assumed to be independent of y, it still ap- 
pears allowable to write (0/Oy) (p’V’) = 
to consider this derivative so small as to be negligible. 


0 or at least 
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It follows that p’U’ may be neglected in comparison 
with pl’, and, hence, 


(0/dx) (BU) + (0/dy) (AV) = 0 (6) 
while one may write 
'V’ = op.U, (7) 
as was done previously in the case of flow along a flat 
plate. 
(A-2) First Integral of Energy 


By comparison of Eq. (4) with Eq. (5) and remem- 
bering that if #7 is a constant it is likely to be true that 
U'H’ = p’'H’ = V'H’ = 0, one finds that 


H = i+ (1/2)U? = const. = 


lp = ly = 
i, + (1/2)U,.2 (8) 


is a particular integral of Eq. (4). This integral cor- 
responds to the case of a nonconducting wall, and the 
following investigation is restricted to this condition. 
(A-3) Turbulent Friction Stress 

From Eq. (3) it may be inferred that the turbulent 


friction stress is given by 


—(pU'V' + Up'V’) (9) 


T= 


If one admits that U’V’, at every given point of the 
flow, depends only upon the average stream character- 
istics, one can generalize the formula used in the pre- 


vious note! of this series. Letting 


U'V’ = (00 /dy) f(U, OU /dy, .. . ) 


and expanding f in a power series of the parameters on 
which it depends, and keeping only the first terms, one 
may write 


U'V’ = —(0U/dy) (C.U* + GU) (10) 


where U* is a velocity depending on the velocity U, 
at infinity but independent, at least formally, of the 
local velocity (. The use of Eq. (10) in the case of 
turbulent flow along a flat plate, with pressure gradient 
equal to zero, led previously' to the replacement of 
U* by the friction velocity in contact with the wall, 


which is written here as 


U* = U,* = Vr, D- (11) 


On the other hand, in proximity to the wall, if dp/dx is 
taken equal to zero, Eq. (10) reduces to 


U'V’ = —(0U/dy) C,U* (10’) 


It follows that C,; must depend on dp/dx and vanish 
along with it. Furthermore, because of the fact that 
(again in proximity to the obstacle S) the physical 
quantities that probably affect C2, at least for the sec- 
tions of the boundary layer far from the one at which 
separation occurs, are the thickness, 6, of the boundary 
layer (at the same cross section, viz., for the same value 
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of x, as that defining the point P of S in the neigh- 
borhood of which we are evaluating C)), the friction 
velocity, U’,*, and the kinematic viscosity v,, then one 
obtains, because of the principle of homogeneity of the 
physical equations, 
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At a sufficient distance from the wall S and therefore 
in proximity to the outer contour of the boundary 
layer, C2, as well as Ci, is to be considered as independ. 
ent of those parameters that are more directly bound 
to the viscosity of the fluid, as, for example, the coeffi- 
cient of kinematic viscosity, v. 


: dp | Vy 1 dU, dv, : ; 
ty eh o27 4 = —ke O, & U,* (11°) It seems allowable, therefore, to set simply 
where y must be sufficiently small and x must not be Cs = be6 and C, = hé (12 
too close to the value defining the location at which 
separation occurs. for v large enough. 
Thus it appears that one may write 
uy = (: ee ) f 1 nf 
= — nyU,* — ko = or y small enough 
Oy lel U, & U,* ; . (13) 
U'V’ = —(OU/dy) 6 (4U,* + he ,) for y large innit 


Now, if —dl,/dx is large enough in order that the related values of y require the use of the first of Eqs. (13), 


ky and 


—k(1/U,) (dU,/dx) (6v,/U',*) 


are of the same order of magnitude. 
thus it will be possible to neglect in 


Eqs. (13) the terms containing k; and /,; with respect to those containing k 


It is further worthy of note that (’, in general, is much larger than U,*; 


and Ay. Thus it will be true that 
uy roy i P l dU, 6v, CE 1 ; | 
= b> - or ysmall enough , 
Ov 5 ae 4,* . . (13’) 
U'V’ = —(dl /dy) h.6U, for vy large enough ) 
This assumption enables one to write 
1 aU. & 5 CO. f il values of y | 
— ere — «p,l’, r sim ré sooty 
, P oy By @ U,” ia ee ee (14 
r= —p(0U OVhdU, — opU,0 for large values of st 


On the other hand, the investigation carried out in the previous report! of this series showed that, in the case of 


flow along a flat plate, the term corresponding to ¢,j,U’,’ in the expression for 7 is of the same order of magnitude 


as that corresponding to —k,yl,* (OU /dy). 


Because the latter term has been neglected in the present investi- 


gation, it is permissible, therefore, to neglect the former also, and it is thus deduced that 


00 1 dU, dv, 
T= —p 2 , 

oy l", dx i” 
r= —p (Ol /dy) bd, 


(for small values of y { (14! 


for large values of y 


The assumption being made here is not essential to the aims of the present investigation, however, since a pro- 


cedure like that followed in the succeeding paragraphs would be possible even if the term involving c, were re 


tained. 


It should be noted that the results now attained, 
particularly in the simplified form of Eys. (14’), may 
hold only at sufficiently large distances from the pro- 
file (cross section of the boundary layer) at which dp/dx 
becomes positive, in order that, away from such a pro- 
file, the influence of the pressure gradient on the motion 
of turbulent vibration will be of such a magnitude as to 
offset the effect of the friction stress at the wall upon 
this motion. On the other hand, at least in regard to 
the first of Eqs. (14’), the remark already made holds 


1.e., this equation is only applicable to those sections 
that are not close to the place where separation occurs. 

The derivation of Eqs. (14’) was just accomplished 
by following a wholly formal procedure based on the 
assumption that the component of the double tensor 
of correlation between the fluctuations of 
the velocity, to which the tangential friction stress 1s 


turbulent 


due, depends solely on the average flow characteristics. 
These relations may, however, be given the following 
physical interpretation : 
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TURBULENT BOUNDARY 

Near the cross section ~ of the boundary layer at 
which dU’,/dx < 0, the influential effect of the friction 
in contact with the wall still prevails in determining the 
correlation between L’’ and V’, especially within the 
region closest to the wall itself (region of wall turbu- 
lence). Moving away from this region and also from 
the above-mentioned cross section L, the influence 
that the pressure variation dp/dx has on the correlation 
increases in importance and probably tends to produce 
a dependence of the correlation itself on the mean 
velocity, analogous to that occurring in the region of 
free turbulence,’ as is substantially characterized by 
the second of Eqs. (14’). These considerations lead 
one to foresee that the nearer the section under investi- 
gation is to the place where separation occurs, the larger 
will be the region within which the second of Eqs. (14’) 
holds. Thus it is implied that, downstream of a cer- 
tain section, at least, it will be possible to assume that 
this equation defines 7 within the whole layer. 


(B) REDUCTION OF THE EQUATION OF MOTION 
WITHIN THE BOUNDARY LAYER TO THE EQUATION OF 
EXTERNAL AND INTERNAL SOLUTIONS 

SEPARATION POINT 


VON MISES 


(B-1) Transformation of Eq. (3) 


In accordance with the method suggested by von 
Mises,? which has been improved by von Karman* 
and applied to the study of the laminar boundary layer 
in flow with a pressure gradient for the case of an in- 
compressible fluid, £ and the reduced stream function 
i = ¥/p.U,L are taken as independent variables, in- 


stead of & = x/Land yn = y/L. Itis found that 


re) gs ra) gina ~— oe 
3 (pUV) = pU U(é, Wa) 
de in p p dE V1 


while it is also true that 


Op Ee Bs. pe , 2x. p? oU? 

dt 2° OE 27 

Since 
p T, = 1/2) U, l— (U,2 ] ) 
ae we —— = ——"—-— (36) 
D. 1 tg — (1/2) UV? 1 — (U?/V,?) 


then, consequently, 


0 /p l l 
(2) - = —— X 
Of \D, 1 — (U,?/V,7) V;? 


From Eqs. (17) and (15) it is found that 


a 1 O02 2B U2\ . d /p 
ee == (1 - “) Ve (2) (18) 
oF 2 0& 2p V;? Of De 


On the other hand, by considering now the external 
solution of the equation of motion within the boundary 
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layer and, consequently, by assuming 7 to be equal 
to the expression given as the second of Eqs. (13), one 


obtains 


Or a rey __ 520 (pe UO ol’ 
= ie hy bU, (> ) = hodU, . (? . ) 


p ae 
On On \ On Pn On \p,” U, OW 
_U, p,? 1 U,2/ V2) al’ O2 /p 
ae hs p ( l Ve p . : (?) 
l a Pa 2 Pa l a Oy" Dp. 


Hence, taking Eq. (1S) into account, one infers that 


fe) A U2 l a2 De eo = 
(2) = ind ep (2) (19) 
Of \B, U,? U. B.? pa? Of? \B, 


Letting 


l U,* _* fa) = ra = 
z ao ES (*) = (7) (20) 
OU,” p.* OF \B, Of \Dp, 
and, hence, 
(1/0) (U,?/U.") (p47/p.7) = dt/dé; (21) 
that is, 
*é 
&= (6U',°5,?/ Ua? pa?)dé (21’) 


Now setting 





| 2 = pp/pe = 1 — (U,?/V;*) ; 
one finds that 
O/U, = (% — 2)'*/29'%(1 — 2)’ (22) : 
and, hence, Eq. (19) yields ‘ 
Oz a (* — *) Cs F 
= (1 — 8) (23) : 
Of; Zs oy,** q 
where i 
Wi* = ~i/ Vk (24) 
(B-2) External Solution of the Equation of von Mises 


It is evident that Eq. (23) has the same form as the 
equation of von Mises’ for the case of laminar flow in 
an incompressible fluid. Thus there may be tound an 
external solution—that is, a solution holding in the 
outer portion of the boundary layer and rigorously on 
the very outer contour of the layer. This solution may 
be easily determined as indicated below. 

Following the procedure suggested by von Karman® 
for treating the analogous equation just mentioned 
[and upon letting z take on the value that it attains on 
the outer contour—i.e., from Eq. (20), = 0] one ob 
tains 


Oz, Of; = 072, Oy, *? 25) 


wherein 2, denotes the function z corresponding to the 
external solution. This equation is the same as the 
equation of heat propagation within a medium having 


unit diffusivity (in accordance with the interpretation 
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given by Lord Kelvin). In order to determine z, un- 
equivocally, it is still necessary to assign the boundary 
conditions. Now, one of these conditions, which has 
already been made use of (viz., the one corresponding 
to ¥i* = ©), is given straight off by zs, = O for all 
values of &. 

In order to establish the other boundary condition 
corresponding to y¥,* = 0, it may first of all be said 
that, according to what has been stated in Section 
(A-3), the region within which the second of Eqs. (14’) 
holds is supposed to be large for those cross sections of 
the boundary layer that are sufficiently far away from 
the initial section =. It follows that, by imposing the 
same boundary condition upon z, that conforms with 
the one existing for the internal solution of Eq. (3), one 
may expect to obtain a solution describing the fluid 
motion within practically the whole turbulent layer 
with sufficient accuracy. If attention is now focused 
upon the boundary of the laminar sublayer, it is evi- 
dent that the values of ¥,* are so small here that it is 
permissible to set them equal to zero. If one then 
imposes the condition that the turbulent shearing stress 
at this boundary be a given multiple m of the viscous 
stress, for m suitably large, the first of Eqs. (14’) vields 


rey Oi 1 dU, dv, 0 roy 
— Pp 'a : = Mp» — 
Pay OC, dx U,* M dy 
whence 
1 dU, U m 
r _ = ¢ (26) 
( e dé l . ke 
and, consequently, for ¥,* = 0, 
U/U, = —c [U,*/5(dU,,/dé) | (26’) 


wherein ¢ is a constant whose value can be obtained 
from experiments. Because of the way in which c 
has been introduced, it is likely to be independent of 


the compressibility of the fluid. Noting that 


z= {1 — (0/02))/((V2/0) — (02/02)] (27) 


This solution is formally given by 
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thus Eq. (27), together with Eq. (26’), affords the means 
of obtaining the distribution of the values of =z for 
y,* = 0; these values are denoted here by 2. The func- 
tion z, satisfying 
conditions just defined, is given by 


tq. (25), together with the boundary 


aa vi¥2/4 (1 — x1 
y,* s e 1 “1) 

= 5 21(%1) — a7 da (28 

“VT 0 ie: — yy 


If the function 2;(x,) may be expanded in a power series 
having the form 


n 
s = > B, x (29 
r =0 


(which is now assumed to be possible), the second mem 
ber of Eq. (28) may be evaluated, and one finds that 


n ok 
2 = » B, 2*" r! é,’ | erfc ( v1 )] (30 
r=0 2V & 


where, according to the notation used by Carslaw, 
erfe x = | — erf x and 7” represents the nth integral 
of the complement of the error function. 


(B-3) Complete Solution of von Mises’ Equation 


For those cross sections that are not too close to the 
section at which separation occurs, the complete solu- 
tion of Eq. (23) may be worked out by the method oj 
successive approximations suggested by von Karman 
in his paper* already referred to above. Writing Eq 
(23) in the form 


Oz c 22 
0£; Oy,*? 


= E ya/s (1 =) 1 H\ 3/ 
— 2 — _ = 2) 6t 
oy, *? 20 


and replacing z in the second member of Eq. (23’) by 


te 
w 


2,, one obtains for z a further approximation as a result 
of solving the equation 
Oz 22 


o = H(s, 
eu 


“a ti ; 
s=23.- I J inf dx; [ -(yi* m)?/4(&1 a0 e (Yi* + )?2/4(& — | A(x, ”) (31 
2V rJ0 0 Va m— X% 


The calculations required by Eq. (31), in order to obtain the values of z in the whole field, are undoubtedly labor- 
ious; on the other hand, if one denotes by /7,,,,, the maximum absolute value of 7, the maximum absolute value 


of z — 2, resulting from Eq. (31) is found to be 


a 


9 


where, for sake of simplicity, the symbol x represents y,* 


” & ae 
— Ze'marz. <= i: / nf dx, le v1 n)2/4(E) 
2V 0 0 
l 


4H mar, & {eric x + 2x(i erfe x) + (27? erfe x) — [222 erfe (O)]} (31 


2V é1. 


x1) __ e vik + n)?/4(i — °| a 
/ : 
Vii — % 


> 


Now, for x comprised in the intervals 1.5 to 1.0 


and from 1.0 to 0.5, the coefficient of 47,7, falls in the intervals 0.2479 to 0.2355 and 0.2355 to 0.18, respec- 
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TURBULENT BOUNDARY 


and because, as a result of calculations to be 


tively ; 
described in the following Sections (C-1) and (C-2), it 
appears that 4H nar. is small besides, it is allowable to 
conclude that, for the sections of the boundary layer 
not too close to the section in which separation takes 


place, even the first approximation z = 2, yields good 
results over the whole interval of variation of y,*. 

In fact the error is certainly somewhat smaller than 
that resulting from Eq. (31’) because not only is /7 not 
constant but it reverses its sign within the range of 
variation of y,* [as one may observe from the calcula- 
tions presented in Sections (C-1) and (C-2)]. Further- 
more, /1 approaches zero for large values of ¥,*—that 
is, it approaches zero right where [according to Eq. 
(31’)] the coefficient of the error attains its largest 
values. 

In proximity to the section in which separation occurs 
and particularly for the purpose of determining the 
separation point, the above given approximation cor- 
responding to the external solution is no longer suffi- 
cient. In order to determine such a separation point, 
a procedure is suggested herewith which is similar to 
that used by von Karman. In fact, the condition to 
be satisfied for ¥,* = 0 is now s = 2; this situation is 
identical with that dealt with by von Karman in his 
above-mentioned paper. 


Letting 
¢ = (% — 2)/2 (32) 
be the independent variable, Eq. (23) becomes 
ore ” 02/0E: 
Me 2-241 — a + QE)? ; “3 
E —2 = oe (33) 
which, for ¢ small enough, is approximated by 
o'¢/owi*? = [A(é)/o""] [1 — Vo] (34) 


The latter equation is obtained by setting 


A(&) = —[207?/2°*(1 — 29)**] (020/081) 


and by replacing the coefficient of A/¢'? in the second 
member of Eq. (33) by a function of ¢ only (in accord- 
ance with the procedure used by von Karman). This 
function of ¢ attains, for ¥,* = 0, the same value as 
the coefficient itself does. Furthermore, let a denote 
a constant or, better, a function of £ only, which is 
left undetermined and which will be at free disposal to 
be employed when convenient in order to ensure the 
connection between the internal solution given by Eq. 
(34) and the external solution presented by Eq. (30). 
Setting 
B, = (o¢ OY" )y.*=0 (35) 


one obtains 


d¢/dvi* = VB? + 4Av/¢ — 2Aat —_ (36) 
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and 


| ; ———— 
” = 5 Bi — Vv B,? + 1A Vv ¢ ated 2Aat ie i 
as 


9 | l 
\ es | aresin — 
Aa Vi1+ (B,?/2A)a 
1—av¢ 
coeuial | (37) 
1 + (B,?/2A)a 


aresin 


In order to join the internal to the external solution, 
the following method is proposed: The external solu- 
tion is assumed to hold up to a value of ¥,* which is 
close to the one (¥,*), for which the velocity profile 
(2, Wi" 2v/t1) has a point of inflection and correspond- 
ing to a value of ¢ (denoted here by ¢,) which is small 
enough with respect to 2/2 so as to permit the replace- 
ment of Eq. (33) by Eq. (84). The values of B, and of 
a are then to be determined. This determination re- 
quires that, for ¢ = ¢,, the following conditions be 
satisfied simultaneously: (a) that ¥,*, obtained from 
Eq. (37), be equal to (¥,*),, corresponding to the value 
¢ = ¢, given by the external solution; and (b) that 
(Of /Oyi*)-—;,, calculated from Eq. (36), be equal to the 
analogous value given by the externa/ solution. 

Since it seems likely that the complete profile of z 
(or of ¢) for every & would have only one point of in- 
flection, the point of inflection related to Eq. (34), given 
by Vet = 1/a, will be associated with a value of 
Wi* > (y4*), or, at least, it will be close to (¥:*),. It is 
hardly necessary to remark that the method described 
here cannot be applied to cross sections of the boundary 
layer which are too close to the one at which separation 
occurs, because for these sections ¢ is everywhere too 
large and substitution of Eq. (34) for Eq. (33) is no 
longer allowable. 


(B-4) 

The location of the cross section at which separation 
occurs can be obtained by following a method anal- 
ogous to that suggested earlier by von Karman. 
In fact, such a section is always defined by the condition 
that B, = 0, since it is true that 


Separation Point of the Stream 


pe U,? 0 U 
0, On U, 


o¢ 1 V;? 2 
oy,* 20 U2 (V,?, U,?) = (U? (7,2) p 


at such a point of separation. 

In the case of supersonic flow, however, it should be 
pointed out that the demarcation of the separation 
point, evolved with aid of the same assumptions upon 
which the theory of the boundary layer is based, is 
somewhat less correct than in the case of subsonic flow. 
In the latter case, in fact, separation ameliorates the 
causes that produce it, because of the fact that, when 
the divergence of the stream tube flowing along the wall 
decreases, the component of grad p in the direction of 
motion decreases as well. In the former case, however, 
separation increases the convergence of the tube, thus 
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RECIPROCAL OF BOUNDARY LAYER THICKNESS __ 1/€ 


( NON- DIMENSIONAL ) 


Fic. 1. Experimental verification that the turbulent shearing stress is proportional to the viscous stress and elicitation of the value of 
c in Eq. (26). 


rendering matters worse. This state of affairs gives 
rise to an inconsistency, and in order to eliminate it the 
flow must become subsonic; however, this alteration of 
the flow cannot happen without the generation of a 
shock wave. This wave must start in front of the sep- 
aration point (which would be determined by applying 
the present theory) in order that the stream may still 
have the kinetic energy necessary to overcome the 
rapid increase of pressure corresponding to the forma- 
tion of the shock wave. On the other hand, one must 
take into consideration the fact that the increase in 
pressure is felt upstream of the wave through the sub- 
sonic boundary layer. 

For this reason, the actual separation point will be 
shifted upstream of that which is predicted by the cal- 
culations based on the foregoing assumptions. 


(B-5) Experimental Verification of Eq. (26) and Deter- 

mination of c 

Before showing how the formulas obtained above may 
be applied to the determination of the motion within 
the boundary layer, it is advisable to ascertain whether 
the experimental tests confirm Eq. (26). As a conse- 
quence of such verification, the value of c will be de- 
duced. With this desirable objective in mind, one 
should remember that in establishing Eq. (26) no as- 
sumption regarding the compressibility of the fluid was 
made; it is therefore reasonable to expect that c does 
not depend on this property of the flow. In order to 
check the above equations, one may thus refer to the 
results obtained for streams having a very small Mach 


Number. In view of this likelihood of slight com- 


pressibility effect on c, the experiments of Gruschwitz' 
and Kehl’ are selected for analysis. 

The values of —(1/U,) (dU,/dé) (U U,*), calcu- 
lated on the basis of the above-mentioned experimental 
results, are shown in Fig. | as a function of the corre- 
sponding values of 1/6 derived from these experiments. 
It clearly appears that the experimental points lie close 
to the straight line. 


1 dU, U 4014+ =0 
U. dt U,* we a 


so that Eq. (26) is sufficiently satisfied, and the con- 


stant is determined asc = 0.14.{ 


(B-6) Approximate Evaluation a priori of 5 and of U,* 


In order to apply Eq. (28) to the determination of 
the motion in the boundary layer, it is necessary to 
know a priori the friction velocity U,* for the various 
values of £, as well as the thickness 6 of the layer. It 
would appear at first glance that these two quantities 
are actually only known a postertort. 

An approximate value a priori for U, 
tained by regarding the friction coefficient c, as constant 


” is now ob- 


+t The length of the experimental wing is taken as the reference 
length, Z, for the experiments by Gruschwitz; while for those by 
Kehl, referring to the rectangular channel designated as AK III, 
the value 1 = 172.5cm. has been assumed. 

t The detailed steps in the calculation of these data, as well as 
those presented in subsequent tables, are contained in an ex 
tended version of this investigation which is available from The 
Johns Hopkins University, Applied Physics Laboratory as refer 
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over the whole portion of the wing along which the 
flow is certainly turbulent and where the further re- 
striction is enforced that such a region be not too close 
to that cross section where separation occurs. The 
coefficient of friction is obtained from the formula that 
is applicable for a flat plate (dp/dt = 0). The Rey- 
nolds Number employed in the formula is taken as the 
one corresponding to the final section of this portion 
of the wing, and the value of the velocity U’,,, is selected 
to lie intermediately between those related to the first 
and last sections of the portion of wing now being con- 
sidered. 

As far as the variation of the thickness 6 of the layer 
is concerned (considering first the case of an incompres- 


sible fluid), one obtains, from the equation 


7, 1 aU, 5* 4 08%) 4 dé** - 
ne ( Z ) (38) 
pl,’ I’, dx dx 
wherein 6* = ¢)6 and 6** = ©26, that 
| 1 dU,. i’ dé 2 3 (38" 
i = = 0 (Cc ZC») C2 38 ) 
2 li, dx dx dx 
Since dc./dx = O if p is constant, let 
de./dx = —b(1/U,) (dU,/dx) 
and hence write 
dé ] 1 dU.q + 2a — b 
= ee 6 (38") 
dx 26; U. dx Ce 


Furthermore, assuming that it is permissible to con- 
sider ¢, ‘Co, as well as the coefficient of (1/U,) (dU,./dx), 
as constant (or at least by assuming it is permissible 
to substitute for the variable values a constant mean 
value), it follows that 


; * 6 aU, 
6=ayx — a& r dx = ayx — aod, X 
ii 


. ‘ dx 
i x dU, 
= dx (39) 
Jf 0 U dx 


€ 


From the above-mentioned experiments by Grusch- 
witz and by Kehl it is possible to deduce the values of 
(4, @, and 6 and, hence, of (c, + 2c — b)/c. The 
variation of the values of the last parameter is not small, 
since its value ranges between 1.32 and 2.158; however, 
this variation is of little consequence as long as the 
second term in the second member of Eq. (38”") is suffi- 
ciently small in comparison with the first. 

Thus it is justifiable to take for this coefficient a mean 
value deduced from those presented above [for which 
purpose the value 1.7 is selected to represent (¢, + 
20. — b)/es Furthermore, because of the fact thot, 
for the case of flow along a flat plate, the term dU, dx 
vanishes, one has simply 6 = 0.38 Yc, x, so it is pro- 
posed to take 


a 
ll 


*x dU, 
0.38 Wc (x _ ivf = dr) (39’) 
0 U, dx 


as a first approximation to the evaluation of 6. 

If the values of U,* and of 6 obtained in a first ap- 
proximation by means of the above-suggested methods 
are compared with those obtained from the experiments 
of Gruschwitz and of Kehl, it appears that, in spite of 
the crude processes employed in obtaining the for- 
mulas that yield U,* and 6, the values deduced are in 
satisfactory agreement with the experimental ones; 
thus it is apparent that these methods may be applied 
to advantage despite their inelegance. 

As a typical example one may cite here a comparison 
between the calculated and the experimental values 
that apply in the case of the experiments by Grusch- 
witz for the wing with 12° angle of attack and for a 
cross-section location having an abscissa § = 0.587. 

It is found that (U,*) ac. = 177 cm. per sec., while 
(U,* exp. = 180 em. per sec.; and (4) care. = 0.626 em., 
while (4),x,. = 0.62 cm. If the stream is compressible, 
Eq. (38) is replaced® by 
d6** L. «@ 1 dU, 


4 
= —— BU .*)8** = é* =-—-c, (40) 
ds" 30 de" * Fe 9% ' 


Letting 6* = c6 and 6** = ©6 and noting that 


dey | d “ 


= —b — (p.U,”) 
dx BU 2dx "* ‘ 
one obtains (by following the same procedure as in the 
case for which the Mach Number is zero) 


* f wy, 
6=avVvlo|1—aq = x dx + 
o U, dx 
~ £88; . ; 
ay = xM? dx (40’) 
o U, dx 


wherein @ is a function of the Mach Number, which 
may be given the value derived in the first report! of this 
series. Experimental data are lacking which would 
enable one to determine approximately the values of the 
coefficients a; and a, analogously to what has just been 
done in the case of an incompressible fluid. It is pro- 
posed, therefore, for the time being, to use Eq. (39’) in 
the case of compressible fluid flows. After having per- 
formed the calculations, one can always check the values 
of 6 and U’,* to see whether they are sufficiently good 
approximations. 


(C) NUMERICAL APPLICATION AND VARIATION OF THE 
DRAG COEFFICIENT ALONG THE BOUNDARY LAYER 


(C-1) Numerical Application; Case of an Incompressible 
Fluid 
In order to compare the numerical results that are 
derived upon the basis of theory just expounded with 
those obtained from experiment, a detailed application 
will now be made for the conditions under which the 
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TABLE 1 a 
Computation of the Velocity Parameter z, As a Function of the Integrated Thickness Parameter £*—Case of Incompressible Flow 
t = 0.396 0.484 0. 587 0.676 0.780 0.872 ; 
i = 1.51 1.41 1.31 1.21 1.14 1.06 
Us 
dU, sealed a il i a o es ° 
li = —3,950 cm. per sec. —3,430 —3,340 —2,800 —2,530 —2,100 
a 
6 = 0.80 X 10°? 1.15 X 107? 1.55 X 107? 1.88 X 107? 2.37 X 107? 2.95 X 107? 
a = 1.83 X 107? 2.29 X 107? 2.66 X 107? 2.76 X 107? 2.08 X 1072 3.31 X 10-2 
as 
fi = 0.0096 * 10-2 0.1928 & 10-2 0.444 K 1072 0.667 X 10° 0.991 XK 107? 1.286 K 10°? 
" = 1,182 0.885 0.629 0.570 0.472 0.424 
10%, = &* = 0.0096 0.1928 0.444 0.687 0.991 1.286 
2 = 0.883 1.205 1.321 1.139 1.076 0.964 


experiments conducted by Gruschwitz were carried out. 
In particular, the test speed is low, so that the numeri- 
cal application applies in this instance to the case of an 
incompressible fluid. The velocity distribution as- 
sumed to exist just external to the boundary layer is 
taken to be identical with that which Gruschwitz de- 
termined to exist around the airfoil he tested, when it 
was at a 12° angle of attack. Then, for p, = p,, Eq. 
(21) yields 


dé, /dé = 6 (U2/U,?) 


The values ot U,/U,, &, dU,/dé, and 6, taken as the 
starting point for the calculations, are given in Table 
1. The values of dé,/dt [obtained from Eq. (21)], 
of & (deduced from the preceding quantities by graphi- 
cal integration), of U;/U, = (U0 U.)y,.*-0 [evaluated 
from Eq. (26) by taking ¢c = 0.14 and U,* = 177 cm. 
per sec. as determined in Section (B-6) ], and, lastly, of 


Because of the difficulty encountered in the deter- 
mination of a polynomial in £,* which approximates 
sufficiently well the function z,(¢,*) corresponding to the 
above values at every point of the interval through 
which &, varies, the values of 2; have been represented 
by a polynomial of third degree just in the first part of 
the interval only 
P3(é,*) = 0.90 + 2.31 &* — 4.32 &*2 + 

2.13 &** (forO < &* < 0.69) (41) 
In order to evaluate z, in the interval for which £* 
external to £,* = 0.69, the quantity 2,(&,*) has been ap- 
proximated by the function 
10.90 + 1.558 &* 

(for 0 < &* < 0.28) 
1.442 — 0.379 &* 

(for 0.28 < &* < 1.286) (42 


iS 


R(é*) = 




















a = (U’.?/U,?) — (U*/ U,’) are also given in the same The broken line defined by Eq. (42) and the cubic de- 
table. fined by Eq. (41) are plotted in Fig. 2, on which the 
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points corresponding to the calculated values of 2; are 


marked by crosses. 


Let 2; = P:(&,*). 


|). * 3 * 
Z, ( = ) = p is Biks'Gr ( v1 ) (30’) 
2V/ é* r=0 2V é* 


Wherein the constant 4, contained in the expression for 


One obtains, from Eq. (30), 


v* includes the factor 10 
tution of &,* for & in Eq. (30) itself. 
The functions i*” erfe x 


* which arises upon substi- 
Furthermore, 
n\(i-" erfe x). 
have already been partially tabulated by Hartree. 


o2n 


( I (x) 


10 


Because the values of such functions for values of x 


larger than those considered in Hartree’s Tables are 
needed, the tabulation has now been extended by com- 
pleting the entries up to the value m = 5. In carrying 
out this work, the following recurrence formula is em- 


ployed: 


2(n + 1) (n + 2) t"** erfe x = 


. > eo 
oe . "~~ “ erfc x 
(1 + 2n + 2x*)t" erfe x — - - 
whence 
G 1+ 4 07)G = , 
Fe 4-1 ( % + 2x"), — G 
i 2n +2 2n + 1 
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The functions G, are given in the Appendix. The 
profiles corresponding to the locations for which £,* 
equals 0.444 and 0.687 are computed by means of Eq. 
(30’) and are plotted in Fig. 3. 
Assuming 2;(£,*) to be given by Eq. (42), one finds for 


z, that 
yi* 


Z ( v1" ) = ().90 erfc [ - + 
*Nav/ge) — OM ere la 7s 
* 
1. 55S8E,* eric , M1 | — (1.558 + 0.379) X 
2V ti, 


* 
v1 | (30”) 
OV t,* — 0.28 


The profiles corresponding to the locations for which 
£,* equals 0.991 and 1.286 are computed by means of 
Eq. (30”) and are plotted in Fig. 3. 

The determination of z, has also been worked out for 
1.4, corresponding to & = 0.95, although for this 
value of no experimental graph of U/U, in terms of 
y is given. This extra computation has been per- 
formed in order to determine the point of separation of 


(¢,* — 0.28) [4 erfc 


ee 
a = 


the stream. 
In Fig. 3 are given the profiles of z, in terms of ¥,* 
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appears that the agreement is generally satisfactory 
within the range of variation of y, and for all the values 
of &* considered. The fact that the deviation of 
(Ze)theor. from the experimental values appears to be 
greater for the larger values of y is due to circum- 
stances already pointed out in Section (B-3), while it 
should be noticed that an error in the experimental 
determination of the velocity, even if small, in this re- 
gion of the boundary layer gives rise to a relatively 
large error in the determination of zg. On the other 
hand, the value of z, = 0 corresponds theoretically to 


¥i = ©, while, experimentally, (’ = U, at an actually 


finite distance from the wetted wall. It is worth men- 
tioning that, at a sufficient distance from the wall, the 
theoretical values of 2, are but slightly different from 
one another for equal values of ¥*/2+/£,* but for vary- 
ing &,*. This tendency to uniformity is also in agree- 
ment with the experimental results. 

From the profiles of z, in terms of ¥1/ 2V/t*, the ve- 
locity diagrams of (’/U, in terms of y/L have then been 
computed. 

From the relation (’/U, = 0y,/0n, one gathers that 


; “W/2VE* (yy 2V t,*) 
2v/ é,* J 7 


7 4 . ; = (43) 
2+~/ £,* deduced from the experiments by Gruschwitz. ' Vian — & 
By comparing these results with those obtained from 
the calculations, one finds that V ho = (0.00625. It where a = U,2/U,?; and, consequently, 
1.333 0.444 0.587 
1.658 wiever d (Wy/2VE*) 0.687 0.676 
n = (1.991 7 for £:* = (0.991 or é = (0.780 (43) 
2.268 “° Vm — & |. 286 0.872 
2.366 1.4 0.95 
In the determination of 7 by means of Eq. (43’) the — or 
graph of zs, in terms of y; 2~/i* has been completed — U.* 
up toy; = 0 by extending the velocity diagram through D = | — | = —2R, —— 2v/t,* 
OY; 2N §,*) v¥i=0 Ue 


the laminar sublayer; it has been assumed here that 
within this layer the variation of the velocity is linear 
and the slope is equal to that at the value of U’,* under 
This assumption yields simply a linear 
On the other 


consideration. 
variation of z, in terms of y¥; 2 t*. 
hand, since 

(0s,/O~i)y,-0 = —2(0/On) (U/U,) 


while 


i Mm (5) a l (= SS 7 U,*? 
pU,  pU.? \dv/y-0 =Re \On Up/,-9 Ua? 


one infers that 
|(O On) (U, U,)I, = R, (7, Ue) 


wherein R, = U,L/v is the Reynolds Number of the 
asymptotic (at an infinite distance) stream. One may 


conclude thus that 


(02, Ovr)y, oe —2R, ts Ta U7) 


and, hence, 


3,080 (40) ( 177° ) P 
— 92 2/ ‘* = 
0.15 3,0802 . 


— 10,640 \/¢,* 


D= 


This expression allows the extension, mentioned above, 
of the velocity profile to the “inside” of the laminar sub- 
layer. 

In order to carry out the integration required in the 
second member of Eq. (43’), the interval of integration 
is subdivided into two parts, the first corresponding to 
the laminar sublayer and the second to the turbulent 
For the former it is true that s = 29 — D (wy + 
varies accordingly 


layer. 
2v/5,*), while for the latter z = 3, 
to the law deduced above. The results of the calcula- 
tions are given in the graphs of Fig. 4, upon which the 
points corresponding to the experimental results are 
1.4 it is found, 


also plotted. For the section at &* = 


- = we A 
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by extrapolation of the linear profile depicted in Fig. 2, 
that z,; = 0.9115 there. 

On the other hand, for the value & = 0.95, corre- 
sponding to this value of £,*, the experimental results 
obtained by Gruschwitz show that U,/U, = 0.95; and, 
consequently, it may be assumed that for this section 
% = 2, closely. From the form of the profile (U/ Us, 7) 
for §;* = 1.4, it appears that this cross section is close 
to the one where separation occurs. The determination 
of the internal solution for this section is completed, 
therefore, by following the procedure suggested in Sec- 
tion (B-3). 

Taking into account Eqs. (36) and (43), one obtains 


I - dy, 0.0125 1 
7 = ] x= F ; 
V2I09 Wt “/2 wWZAe 


2A 
| aresin — 
V 4A? + 2AaB;’ 


2A(1 — ave) | 43" 
; 3°) 
V 4A? + 2AaB,’ 


aresin 


0.00625 l 


0.02 = 


Now the function consisting of the second member 
of this equation decreases as B, increases, and since, 
for B, = 0, its value is already equal to 0.01715, it is 
obvious that the junction is impossible. It therefore 
appears evident that separation has already taken place 
by the time this cross section is reached. Then, if it is 
assumed that B, = 0 and it is admitted that U/U, 


(aresin 
V 0.526-6.26B,? V 14+ 22.95B,? — 273B,! 


From the profiles determined by Gruschwitz, the value 
of A is derived as 


1-2 om _ 2 ae: (5°) dg 
~ Gage 9 ‘ dé \U,/ déi* — 


~ a 


(ee \ es) = 
= (0.107 
1.41 3.31 


For the point of connection between the external and 
the iniernal solution (see the profile of Fig. 4), the one 
for which the coordinates are U/U, = 0.4, n = 0.02 is 
selected; this location is somewhat close to the in- 


~ 


flection point. 
From the same profile one perceives that, at this in 
flection point, (d/dn) (U/U,) = 30.2, and, since 


le didn d i 
dg _ dn % _ 900605 


dy,* dn dy, dy,* dn Us 
it results that V Be + 4A Y¢; — 2Aah, = (0.00625) X 
(30.2) = 0.189 or a = 4.92 — 58.5B,";_ while, from Eq. 


(43”), we obtain 


—arcsin 


—0.39 + 16.55B,? ) 
V1 + 22.95B,2 — 273B,4 


vanishes for 7 = 0.02 — 0.0175 = 0.00285, the value of 
a is obtained as 4.92. By use of this value of a and by 


means of Eq. (43”), values in Table 2 are found. 


The graph describing this internal solution is shown 
as a dotted line in Fig. 4, on which the experimental 


profile obtained for £ = 0.98 is also drawn (dash-dot 
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TABLE 2 
Velocity Distribution for Internal Region of Boundary Layer 


U/U, = 0.4 0.3 0.2 0.1 0.05 0 
n = (0.2 0.01675 0.01382 0.01029 0.00802 0 


TABLE 3 


Computation of the Velocity Parameter 2; As a Function of the Integrated Thickness Parameter, £;*—-Case of Compressible Flow 


780 





& = 0.396 0.484 0). 587 0.676 0 
6 = 1.048 X 10 1.505 X 107? 2.03 X 10°? 2.46 X 107? 3.86 X 10-2 
0. 
— = 1.51 1.41 1.31 1.21 1.14 | , 
Ua 
Pe we on ae 
= 0.1175 0.2298 0.379 0.556 0.696 
Pa 
dé, . : : 
at = 0.033 K 107? 0.158 XK 107? 0.498 XK 107? Lia xX wo 2.43 X 107? 
i* = 0 0.0077 0.0384 0.111 0.292 
U 
=) = 0.598 0.481 0.365 0.36 0.254 
0.? acs e ‘ " 
— = ().707 0.617 0.5382 0.455 0.404 
V7? 
zi = 0.606 0.557 0.495 0.42 0.388 
} 
line). The qualitative and the quantitative agree- attains the same values of (’, l, at the same loca 


ment between the two curves appears to be excellent. tions, &, as considered in Section (C-1), for the region | 
comprised between the velocity just external to the 
(C-2-1| Numerical Application: Case of a Compressible layer and the velocity at infinity. 
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In order to calculate the approximate values of 6 
(40’), the following values are em- 


by means of Eq. 
1.7, and a2 = 0. The same 


ployed: G¢ = 0.5, a= 
value as employed in Section (C-1) is assumed for Cy. 
The value of a is the same as that obtained in the first 
report! of this series, at 7, = 1.5, for the case of turbu- 
lent flow along a flat plate. 

The values of 6 in the case now considered are ob- 
tained at once from those calculated in Section (C-1) 
just through multiplication by 0.56/0.38 = 1.31. Be- 


cause 


The profile showing the variation of 2 


2,(£,*) was approximated (used in the determination of 2,) is also given. 


§0.600 — 2.8648:* 
10.428 — 0.166,* 


L(é,*) = 


LAYER 


2, with respect to &, 
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the values of di,/dé may be easily obtained. The re- 
sults of the calculations are listed in Table 3, where 
values of &* = 10°&; U,/U, (assuming U,*/U, = 
177/3,080 = 0.0575, as in the previously considered 


the 
vile 


case of an incompressible fluid) ; 
U? U l 
2/(y — 1)] (AL) 


and 


are also given. 


* is given in Fig. 5, where the broken line with which 


This broken line describes the function 


(forO0 < &* < 0.0636) 
(for 0.0636 < &* < 0.292) 


-— § 


The reduced velocity z, then takes on the following form: 


* 

z. = 0.600 erfc 
9 ¢.* 

“V e1 


* 
s, = 0.600 erfc — 2.864£,* | 42? erfc 
9 f, + » 


aad s 


yi" 
2y =3|+ 


(2.8645 — 0.1657) (&* — 0.0636)| 4°erfe : 
2% 


— 2.864&* | 47? erfc P= 


vi" ‘ id a _ 
x= (for 0 < &* < 0.0636) 


“Vil 


y* = Oya * 90° 
-| (for 0.0636 < & < 0.292) 
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Supersonic flow velocity profiles calculated for three of the same locations on airfoil as previously computed for the 


incompressible case. 








474 JOURNAL OF THE AERONAUTICAL SCIENCES—JULY, 1951 


The calculation of z, for the different values of the parameter y,*/2+/£¢,* has been carried out for £,* = 0.0384, 
0.111, and 0.292. Starting from these values of z,, the corresponding values of l’/ L’, will now be calculated. Let 


Z = (U2/U,?) — U?/U,? 
then one obtains 


Z = (V?/U,?) [1 — (U,2/V?)] [2/1 — 2) (44) 


and, consequently, upon noting that V,?/U,"? = 3.22, one has 


1.51 ] 40.0384 
Z = 1.755 z-/(1 — Ze) for &* =<0.111 (45) 
1.92 j 0.292 
The graphs of Fig. 6 depict the final result of calculating Z in terms of ¥,/2+ €:* = (1.00625 (W.*/2+/,*), where 
the same value of /, is employed as that derived in Section (C-1). 
Letting Z) = U,?/U,?, one has 
“W/2VE* | y ‘W/2VeE* yp | | y 
7 = over | - a ( ‘ ) — over | ; : a( 1 ) 
0 pl 2v/t,* Jo BIl-—-zV¥2—-Z \2QVii* 
pa Uy 
35 . . 0.038 
1.035) W/2V i" | vy ” iS4 
1.2 : =e = = for &* 0.111 (46) 
[55 \ /0 l—zsvV/Z—Z “Vv &1 ho 299 


The evaluation of » from Eq. (46) is accomplished by following the same procedure as described in the preceding 
Section. The slope of the straight-line portion on the graph of Z versus ¥,/2/ &,*, which correspond to the laminar 


sublayer, is now 


(02 /0p1)y,=0 = —2(Pa/ Be) (Be/By)” Rp (U,**/ U2") 
and since 
(B./pp)? = 1/[1 — (U2/V,?)}? 
and 
R, = LU,4/vp = Ra (Be/ pa) (Ma/ He) (Ve/Yp) 
while 
Me/B. = (T,/T.)"= 1 + 0.2M,? E a (7) | 
U2) Js 
and 
v/v, = [1 — (0.2/V,2)]'*" 


it 1s found that 
oZ ] , U,*? | 
= — (2)27/:,*R, — rene : 
vi U? ua . U, : 
le —— | | p=-——3 1ta-eeMels ~ — | 
=-V i ¥ 0 I lv" l 1 
It appears from the form of this expression that the slope of this portion of the curve is obtained from that cal- 


culated in the preceding Section, multiplied by the term standing between the brackets. If this adjusting term 1s 


denoted by the symbol £, then, for &,* = 0.0384: 


+ (0.2)2.25(1 — 1.716)°-* 


for &,* = (0.111: 





T | l 
-— —— (1 — 0.455)'-5 = 1.39 
L(1 — 0.455)2] 1 + (0.2)2.25(1 — 1.465)°-% 
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0384, | | T | | | | | | 
Let 
2.6 { | 
24 
(44) 5.107 
2.2 
2.0 
(45 
09 18 
vhere 
08 16 
07 14 
5 * 107 
q 
06 12 
5**\0? 
; 05 10 + 
(46) 
} 
04 08 | 1 | 
ding | 6 (BOUNDARY LAYER THICKNESS) 
in .* 
a 03 06 6" (DISPLACEMENT THICKNESS) 
ee = 
5°" (MOMENTUM THICKNESS) 5** 
0.2 04 
} 
Ol 0.2 —— ¢ 
call a See | 
0.50 0.55 060 0.65 0.70 075 0.80 0.85 
LOCATION ALONG AIRFOIL - & 
Fic. 8. Theoretical variation along airfoil of boundary-layer thickness parameters for the case of compressible (AZ, 1.5) flow 
} 
for &* 0.292 
| | Sas " 
KE = - - ——— (1 — 0.404)": = 1.27 
(1 — 0.404)2] 1 + (0.2)2.25(1 — 1.3) 
where 7 has been assumed to be equal to 0.75. 
The profiles giving the law of variation of U/U, in terms of n, calculated from Eq. (46), are presented in Fig. 7. 
From these graphs it appears that the value of 5 assumed a priori for evaluating &, as well as for determining 7), 
is somewhat larger than that resulting a posteriori from the above calculation. The magnitude of 6 evinced in 
al- Fig. 8 is practically equal to that obtained in the case of an incompressible fluid; this agreement justifies the state- 
Ms ment made at the end of Section (B-6). 


(C-2-2) Drag Law 


The graphs presenting U/U, versus n and the values of j/ p, calculated in the preceding Section allow us to de 


termine 6* and 6**. By graphical integration it is found that, for &,* = 0.0384: 


= 0.31: = 0.0718: 5 = 1.55 X 10-2; 5* = 0.481 X 10-2: 6** = 0.1114 X 10 
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for &:* = 0.111: 


* 5k 
= 0.285; 


\| 


0.1459 K 10 


’ 


= 0.0775; 5 = 1.88 X 10-2; 6* = 0.536 x 10-7; 6** 


for £:* = 0.292: 


6* oe” “ 

5 = 0.3325; — = 0.091; 5 = 2:37 «K 10-*- o* = 0.780 X 10°: o** = 6.216 X 10 
6 

The values of 6, 6*, and 6** in terms of ~ are presented in the plots of Fig. 8. From these diagrams one can 


estimate that, for § = 0.587: 


Sé 


as** Bt! 4° dix, UF 1 dU, 
= 0.0023; while “—*- = 0.65; (? ) “ase —— = «oon 
dé Pava™ dé Pa l - ff € dé 
for & = 0.676: 
da” B.U,” ad (3, U,* 1 dU 
= 0.00534; while * — = 0,815; (2: “) = = ‘= —0.75 
dé pee? dé \ pq Us? CU, dé 
for — = 0.780: 
dé** 5, U,? d-73; U7 l au, 
—— = (),00835; while ~ = 0.905; (? :) = 0.38; = = —().722 
dé pe US dé \ pq U?? U, dé 
Appendix 
Values of the Functions G,,(y,*/2 VN &*) Tt 
* 
v1 Ge Gy Gs Ge Gi. G; 
bv &,* 
0 l 1 1 1 1 l 
0.05 0.943628 0.8921 0.8591 0.8337 0.8124 0.7940 
0.1 0.887537 0.7936 0.7363 0.6933 0.6585 0.6290 
0.15 0.832004 0.7040 0.6292 0.5749 0.5321 0.4967 
0.2 0.777297 0.6227 0. 5862 0.4753 0.4285 0.3907 
0.25 0.723674 0.5491 0.4556 0.3922 0.3449 0.3077 
0.3 0.671373 0.4828 0.3861 0.3226 0.2764 0.2409 
0.35 0.620618 0.4233 0.3263 0.2647 0.2212 0. 1884 
0.4 0.571608 0.3699 0.2748 0.2165 0.1764 0.1471 
0.45 0.524518 0.3223 0.2308 0.1765 0.1402 0.1141 
0.5 0.479500 0.2799 0.1934 0.1435 Oo. 283% 0. 0884 
0.55 0.436677 0.24238 0.1616 0.1164 0.0878 0.0683 
0.6 0.396144 0.2090 0.1344 0.0941 0.0692 0.0526 
0.65 0.357971 0.1798 0.1116 0.0758 0.0544 0.0404 
0.7 0.322199 0.1541 0.0992 0.0609 0.0426 0.0309 
0.75 0. 288844 0.1316 0.0761 0.0487 0.0332 0.0235 
0.8 0.257899 0.1120 0.0625 0.0389 0.0257 0.0176 
0.85 0. 229332 0.0950 0.0512 0.0309 0.0198 0.0132 
0.9 0. 203092 0.0803 0.0419 0.0246 0.0154 0.0100 
0.95 0.179109 0.0677 0.0341 0.0196 0.0122 0.0081 
1.0 0.157299 0.0568 0.0277 0.0154 0.0093 0.0060 
4 0.119795 0.0396 0.0180 0.0094 0.0052 0.0029 
1.2 0. 089686 0.0272 0.0116 0.0058 0.0032 0.0020 
3 0.065992 0.0184 0.0074 0.0035 0.0020 0.0014 
1.4 0.047715 0.0122 0.0046 0.0022 0.0012 0.0009 
1.5 0.033895 0.0080 0.0028 0.0012 0. 0006 0.0004 
1.6 0.023652 0.0052 0.0017 0.0008 0.0004 
es 0.016210 0.0033 0.0011 0.0005 
1.8 0.010909 0.0021 0.0006 
1.9 0.007210 0.00138 0.0004 
2.0 0.004678 0.0008 0.0002 
= % 0.002979 0.0005 0.0001 
22 0.001863 0.0003 
+ Note: The values of Go and G; are given here for sake of convenience. They have been reproduced from Carslaw and Jaeger, 


Conduction of Heat in Solids, page 373, Appendix III. The values of G2, G3, G4, and G; are given to four decimal places, although the 
fourth place is not accurate because of the inaccuracy inherent in the fifth decimal place in the expansion for G;, from which the G’s of 
higher ordered subscripts are derived by means of a recurrence formula. This inaccuracy is aggravated at the larger value of 
y1 4 2 =. 








can 


Be 
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From Eq. (40) it follows that, for & = 0.587: 


2.4 ] 
0.0023 + — — (0.001114) — 0.898(0.00481) = 0.0024 = 5 


-00 yA 


whence c, = 0.0048 and U,*/U, = (U,/Us) » ‘(0.0024 = 0.0641. 


For & = 0.676: 


1.16 
0.00534 + (0.001459) — 0.75(0.00536) = 0.00339 


0.815 


so c, = 0.00678 and U,*/U, = 0.0705. 
For ¢ = 0.78: 


0.00835 + om (0.00216) — 0.722(0.00789) = 0.00392 
£U0 
so c, = U.00784 and U,*/ U, = 0.0713. 

The values of U,*/ U, just found are slightly larger than those assumed in carrying out the first approximation, 
but, on the other hand, it also turns out that the assumed values of 6 are larger than those calculated. Because 
these trends yield values of 6* and 6** larger than the correct values, it may be presumed t hat the deviation be- 
tween the assumed values of L’,*/ U, and the actual ones is smaller than that now in evidence. 


Turbulente Reibungsschichten, Ingenieur Archiv, Volume xiii, 
p. 293, 1943. 
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Tables of the Theodorsen Circulation 
Function for Generalized Motion’ 


YUDELL L. LUKEt ann MAX A. DENGLER? 
Midwest Research Institute 


(1) SumMMaARY 


Important in all stability, aeroelastic, and flutter problems is 
the Theodorsen C(k) function. Tables of this function are well 
known if the motion is harmonic. Heretofore, tables of the 
Theodorsen function for generalized motion have not been avail 
able. In this report, the latter tables are given for a range that 
should be sufficient to meet the needs of aerodynamicists. W. P. 
Jones has given a discussion of the generalized Theodorsen func- 
tion, but his results appear incorrect, for, in his analysis, the 
Theodorsen function as derived by approaching harmonic motion 
from the side of stability and from the side of instability is dis 
tinct. If the Jones analysis were true, the nonuniqueness of the 
Theodorsen function would imply that the Bessel functions are 
not continuous, which is indeed false. The difficulty is clarified 
in this paper. 


(2) INTRODUCTION 
tm OF THE MOST BASIC FUNCTIONS in subsonic 
unsteady aerodynamics is the Theodorsen func- 


tion.' 
bility, aeroelastic, and flutter problems. 


This function is important in all types of sta- 
In the case 
of incompressible flow, it is needed to compute the 
values of the flutter coefficients. 
ative method, Dietze? uses the function to calculate 
An alternative 


Employing an iter- 


flutter coefficients for subsonic flow. 
iterative procedure for calculating subsonic flutter 
coefficients has been suggested by Miles,* and here, 
too, the starting point is a knowledge of the Theodorsen 
function. Another application is found in the aero- 
dynamic theory of oscillating wings of finite span.‘~® 
The curt discussion rendered is by no means complete 
but is presented only to show the significance of the 
Theodorsen function. 


(3) DEVELOPMENT OF THE THEODORSEN 
CIRCULATION FUNCTION FOR GENERALIZED MOTION 


In the treatment of aerodynamic oscillatory prob- 
lems, it is assumed that the motion follows the law 
e*'’ where ¢ is the time variable and 


ik = uw + Ww (1) 


Thus, w gives the frequency of motion, while » defines 
the damping. If wu < 0, the motion is convergent or 
stable; if « > 0, the motion is divergent or unstable. 

Received July 5, 1950. Revised December 22, 1950. 

* The authors desire to thank Elizabeth Kahn, Gloria Flake, 
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When uv = 0, the motion is of constant amplitude and 
is said to be neutrally stable. 


Heretofore, for the most part, the nature of the 
Theodorsen function has only been considered when 
ua = 0. Ina paper by Jones,’ the generalization of this 
function is studied. His conclusions are rather dis- 
turbing, and it seems appropriate at this time that the 
difficulties be clarified. To accomplish this effectively, 
a derivation of the Theodorsen function C(k) for gen- 
eral motions will be presented. For a starting point, 


consider the equality*® 


where it is supposed that 
U - Une b x) + o| (3 


and that w > 0. It is well known in the theory of 


Bessel functions? that 


r zcosh y 
K,(z) = / e °°" cosh ny dy (4 
0 


which is valid if |arg 2}< 2/2. This requires that 
Re(z) > 0. Let 


o 


s= ik is 


Thus, if Re(z) > 0, it follows that u > O, so that the 
motion is assumed to be divergent. Combining Eqs. 


(2) through (5) with the well-known relation 
K,(2) = (9/2) te?" HT, (iz) (6) 
it is readily deduced that 
C(k) = My (k)/ (Eh (k) + ify (k) | (7) 


which is the generalized Theodorsen C(k) function. It 
is true that, in the development of Eq. (7), it was neces- 
sary to require that » > 0. However, Eq. (7) has no 
need for such a restriction, and therefore, by the 
method of analytic continuation, one can argue that 
Eq. (7) is valid for all y's. 

In the preceding analysis the frequency of the motion 
was taken as positive. If, in Eq. (3), 7is replaced by —7 
so that Eq. (1) is replaced by its complex conjugate, 
then it is easy to show that the Theodorsen function is 
given by the complex conjugate of Eq. (7), which is 


C(k) = Th? (k) [2 (Rk) — ie? (R)] (8) 
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93993 
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TABLE 1 
Tables of the Theodorsen Circulation Function for Generalized Motion 
0° ‘ 0 = 5° 6 = 10° 
— G(pe?) F( pe’) — G( pet) F( pet? — G( pet?) 

0 1.0 0 1.0 0 

0.04565 21 0.98549 60 0.04717 45 0.98871 13 0.04851 44 
0.07520 79 0.96833 68 0.07839 16 0.97323 65 0.08136 36 
0.07971 35 0.95048 58 0.10269 36 0.95637 13 0.10729 86 
0.11600 13 0.93260 31 0.12225 04 0.93900 62 0.12839 71 
0.13064 44 0.91503 33 0.13821 25 0.92162 O1 0.14576 94 
0.14259 44 0.89797 i9 0.15132 68 0.90468 45 0.15909 51 
0.15238 17 0.88153 O00 0.16212 89 0.88782 24 0.17206 30 
0.16040 21 0.86576 68 0.17102 38 0.87169 73 0.18192 70 
0.16696 15 0.85070 83 0.17832 92 0.85618 78 0.19006 56 
0.17230 22 0.83635 89 0.18429 96 0.84132 73 0.19674 41 
0.17661 95 0.82270 87 0.18914 23 0.82712 77 0.20218 13 
0.18007 27 0.80973 89 0.19302 83 0.81358 69 0.20655 96 
0.18279 33 0.79742 48 0.19609 96 0.80069 26 0.21008 19 
0.18489 04 0.78573 82 0.19847 53 0.78842 69 ().21272 77 
0.18645 56 0.77464 93 0.20025 57 0.77676 64 0.21475 66 
0.18756 59 0.76412 76 0.20152 57 0.76568 61 0.21621 22 
0.18828 65 0.75414 26 ().20235 76 0.75515 91 0.21717 45 
0.18867 27 0.74466 47 0.20281 28 0.74515 89 0.21771 24 
0.18877 18 0.73566 50 0.20294 44 0.73565 74 0.21788 42 
0.18862 42 0.72711 59 0.20279 72 0.72662 89 0.21774 11 
0.18826 47 0.71899 12 0.20241 10 0.71804 69 0.21732 70 
0.18772 32 0.71126 59 0.20181 90 0.70988 75 0.21668 02 
0.18702 55 0.70391 66 0.20105 08 0.70212 62 0.21583 31 
0.18619 40 0.69692 11 0.20013 03 0.69474 138 0.21481 47 
0.18524 80 0.69025 86 0.19908 O08 0.68771 04 0.21364 96 
0.18420 43 0.68390 96 0.19792 10 0.68101 42 0.21235 97 
0.18307 76 0.67785 59 0.19666 75 0.67463 30 0.21096 37 
0.18188 07 0.67208 05 0.19533 48 0.66854 96 0.20947 84 
0.18062 46 0.66656 71 0.19393 53 0.66274 63 0.20791 75 
0.17931 91 0.66130 10 0.19248 04 0.65720 81 0.20629 4: 
0.17659 29 0.65145 54 0.18944 11 0.64686 76 0.20290 27 
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0.51295 
0.50874 
0.50628 
0.50471 

0.50367 
0.50239 
0.50168 
0.50124 
0.50095 
0.50076 
0.50061 


86 
60 
00 
79 
63 
92 
75 
79 
21 
61 
22 
16 
71 
95 
99 
66 
66 
86 
83 
18 
32 
48 
40 
00 
92 
09 
72 
46 
70 
96 
O09 
78 


~ 


frequency of the motion is negative. 
involving the Theodorsen function, current 
practice is to assume that the frequency part of the 
In this manner, many difficulties 


motion is positive. 


can be avoided. 


0.17375 
0.17085 
0.16792 
0.16498 
0.16205 
0.15915 
0.15629 
0.15347 
0.15070 
0.14405 
0.13778 
0.13191 
0.12641 
0.12129 
0.11650 
0.112038 
0.10784 
0.10393 


0.10027 2 
0.07356 - 


0.05769 
0.04729 
0.04000 


0.03462 2 


0.08049 
0.02459 


0.02059. 5: 
0.01770 


0.01552 
0.01381 
0.01244 


In all 


If the motion is nonoscillatory but divergent so that 
® = Oand yu > O, then, 


80 0.64244 24 0.18628 06 0.63741 98 0.19937 58 
75 0.63417 33 0.18304 78 0.62877 OO 0.19576 92 
44 0.62657 03 0.17978 00 0.62083 48 0.19212 55 
40 0.61956 55 0.17650 60 0.61354 07 0.18847 72 
57 0.61309 89 0.17324 76 0.60682 32 0.18484 91 
43 0.60711 77 0.17002 15 0.60062 52 0.18125 98 
10 0.60157 53 0.16684 03 0.59489 64 0.17772 35 
40 0.59643 05 0.16371 30 0.58959 22 0.17425 06 
95 0.59164 67 0.16064 67 0.58467 30 0.17084 84 
33 0.58106 19 0.15327 49 0.57383 76 0.16268 29 
52 0.57211 97 0.14634 82 0.56474 49 0.15502 89 
06 0.56450 10 0.13987 06 0.55705 02 0.14788 87 
89 0.55796 OO 0.13382 82 0.55048 91 0.14124 39 
05 0.55230 52 0.12819 74 0.54485 59 0.13506 60 
24 0.54738 56 0.12295 10 0.53998 90 0.12932 26 
Ol 0.54308 11 0.11806 02 0.53576 OO 0.12398 02 
96 0.53929 48 0.11349 72 0.53206 60 0.11900 61 
77 0.538594 82 0.10923 53 0.52882 39 0.11436 94 
69 0.10524 96 0.52596 57 0.11004 17 
65 0.07642 70 0.50980 34 0.07900 88 
13 0.50843 14 0.05951 18 0.50364 62 0.06104 38 
69 0.50490 39 0.04853 74 0.50087 53 0.04950 59 
39 0.50296 19 0.04089 26 0.49950 39 0.04153 138 
23 0.50180 75 0.03528 34 0.49878 89 0.03571 63 














61 0.50108 18 0.03100 27 0.49840 94 0.03130 O07 
86 0.50028 47 0.02491 64 0.49811 86 0.02505 83 
55 0.49990 47 0.02080 81 0.49808 84 0.02086 98 
50) 0.49971 11 0.01785 40 0.49814 96 0.01787 12 
20 0.49960 99 0.01562 99 0.49824 16 0.01562 10 
61 0.49955 78 0.01389 59 0.49834 O08 0.01387 15 
67 0.49953 29 0.01250 69 0.49843 75 0.01247 29 

C(u) = Ki(u)/[Ki(u) + Ko(u)] (9) 


If the motion is nonoscillatory and convergent, then 
Eq. (9) must be replaced by 
Ki(u) + atli(u 


= - (10) 
Ki(u) nial Ko(u) 5 mi[1;(u) + To(u) | 


C(u) 


where u is now treated as a positive number. 











480 JOURNAL OF THE AERONAUTICAL SCIENCES—JULY, 1951 
TABLE 2 
Tables of the Theodorsen Circulation Function for Generalized Motion 
—— = 15°———~ ————# = 20°-—.. ———P = 25°-—___-—. —§ = 30° a 

p F( pe?) — G(pe*®) F( pe?) — G( pe?) F(pe*®) — G(pe*#) F( pe*®) — G(pet#) 
0 1.0 0 1.0 0 L.O 0 1.0 ) 

0.01 0.99206 O1 0.04965 82 0.99553 36 0.05059 20 0.99912 17 0.05130 20 1.00281 29 0.05177 39 
0.02 0.97842 52 0.08409 88 0.98390 25 0.08657 02 0.98966 56 0.08874 85 0.99570 96 0 09060 23 
0.03 0.96268 22 0.11169 82 0.96943 31 0.11585 838 0.97663 69 0.11974 08 0.98430 49 0. 12330 07 
0.04 0.94594 13 0.138441 22 0.95348 98 0.14026 15 0.96153 35 0.14590 42 0.97025 50 0.15129 23 
0.05 0.92881 33 0.15329 25 0.93666 02 0.16075 28 0.94521 21 0.16811 41 0.95452 27 0.17533 21 
0.06 0.91167 80 0.16904 21 0.91956 81 0.17799 25 0.92823 74 0.18697 28 0.93775 96 0. 19594 84 
0.07 0.89478 31 0.18218 25 0.90248 35 0.19248 14 0.91000 37 0.20294 77 0.92048 42 0.21356 23 
0.08 0.87829 12 0.19312 30 0.88562 69 0.20462 04 0.89379 32 0).21642 38 0.90289 20 0.22853 18 
0.09 0.86230 80 0.20219 53 0.86914 99 0.21474 18 0.87680 76 0.22772 70 0.88538 98 0 24116 99 
0.10 0.84689 91 0.20967 30 0.85315 52 0.22312 47 0.86019 O9 ().23713 84 0.86811 75 0 25175 38 
0.11 0.83210 18 0.21578 56 0.83770 94 0.23000 71 0.84404 34 0.24490 23 0.85121 36 0.26053 04 
0.12 0.81793 37 0.22072 61 0.82285 30 0.23559 28 0.82843 28 0.25123 11 0.838477 75 0.26771 92 
0.138 0.80439 79 0.22465 88 0.80860 73 0.24005 68 0.81340 17 0.25631 O8 0.81887 77 0.27351 58 
0.14 0.79148 73 0.22772 40 0.79497 94 0.24355 O04 0.79897 44 0.26030 37 0.80355 90 0.27809 34 
0.15 0.77918 84 0.23004 12 0.78196 61 0.24620 41 0.78516 10 0.26335 16 0.78884 79 0.28160 58 
0.16 0.76748 18 0.23171 38 0.76955 67 0.24813 10 0.77196 11 0.26557 87 0.77475 68 0.28418 88 
0.17 0.75634 62 0.238282 99 0.75773 55 0.24942 92 0.75936 65 0.26709 34 0.76128 74 9 28596 23 
0.18 0.74575 73 0.23346 61 0.74648 31 0.25018 36 0.74736 37 0.26799 04 0.74843 35 0). 28703 21 
0.19 0.73569 O07 0.23368 82 0.73577 78 0.25046 82 0.73593 49 0.26835 27 0.73618 27 0 28749 10 
0.20 0.72612 11 0.23355 34 0.72559 64 0.25034 72 0.72506 OO 0.26825 25 0.72451 86 0). 28742 10 
0.21 0.71702 36 0.23311 13 0.71591 54 0.24987 64 0.71471 68 0.26775 33 0.71342 18 0 28689 38 
0.22 0.70837 38 0.23240 46 0.70671 O08 0.24910 46 0.70488 26 0.26691 02 0.70287 07 0.28597 24 
0.23 -0.70014 78 0.23147 08 0.69795 89 0.24807 42 0.69553 40 0.26577 16 0.69284 27 0.28471 22 
0.24 0.69232 28 0.23034 238 0.68963 66 0.24682 20 0.68664 77 0.26437 98 0.68331 46 0.28316 16 
0.25 0.68487 68 0.22904 75 0.68172 14 0.24538 06 0.67820 O07 0.26277 17 0.67426 30 0.28136 32 
0.26 0.67778 73 0.22761 02 0.67419 16 0.24377 79 0.67017 O05 0.26097 97 0.66566 46 0.27935 438 
0.27 0.67103 96 0.22605 42 0.66702 66 0.24203 86 0.66253 55 0.25903 21 0.65749 67 0.27716 74 
0.28 0.66460 96 0.22439 60 0.66020 63 0.24018 43 0.65527 47 0.25695 36 0.64973 71 0.27483 11 
0.29 0.65848 13 0.22265 30 0.65371 22 0.23823 38 0.64836 82 0.25476 58 0.64236 46 0.27237 04 
0.30 0.65263 79 0.22083 93 0.64752 62 0.23620 35 0.64179 67 0.25248 77 0.63535 85 0.26980 72 
0.32 0.64174 38 0.21704 88 0.63601 22 0.23195 93 0.62958 74 0.24772 47 =: 0.62236 79 ~=—s- 0.26 444 75 
0.34 0.63181 18 0.21310 74 0.62554 O1 0.22754 68 0.61851 24 0.24277 37 0.61061 91 ).25887 79 
0.36 0.62273 97 0.20907 85 0.61599 95 0.22303 82 0:60845 16 0.23771 79 0.59998 06 ) 25319 438 
0.38 0.61443 76 0.20501 05 0.60729 25 0.21848 92 0.59929 81 0.23262 O08 0.59033 46 0) 24746 97 
0.40 0.60682 60 0.20094 04 0.59933 26 0.21394 16 0.59095 68 0.22753 03 0.58157 62 0.24175 91 
0.42 0.59983 47 0.19689 62 0.59204 34 0.20942 74 0.58334 38 0.22248 23 0.57361 26 0.23610 37 
0.44 0.59340 19 0.19289 92 0.58535 70 0.20497 O05 0.57638 46 0.21750 51 0.56636 12 0 23053 41 
0.46 0.58747 28 0.18896 49 0.57921 36 0.20058 84 0.57001 31 0.21261 71 0.55974 88 0). 22507 25 
0.48 0.58199 88 0.18510 51 0.57356 00 0.19629 41 0.56417 10 0.20783 32 0.55371 06 0.21973 48 
0.50 0.57693 69 0.18132 79 0.56834 92 0.19209 66 0.55880 64 0.20316 32 0.54818 94 0.21453 19 
0.55 0.56584 42 0.17227 91 0.55699 58 0.18206 19 0.54719 42 0.19202 44 0.53632 72 0.20215 34 
0.60 0.55660 58 0.16382 00 0.54762 14 0.17270 89 0.53770 03 0.18167 66 0.52673 87 0.19069 58 
0.65 0.54884 81 0.15594 97 0.53981 92 0.16403 24 0.52987 92 0.17210 81 0.51893 34 0.18013 89 
0.70 0.54228 50 0.14864 44 0.53327 82 0.15600 18 0.52339 18 0.16327 98 0.51254 02 0 17043 22 
0.75 0.53669 48 0.14186 98 0.52775 87 0.14857 51 0.51797 76 0.15514 O1 0.50727 51 O 16151 21 
0.80 0.538190 41 0.13558 68 0.52307 35 0.14170 59 0.513848 45 0.14763 32 0.50291 88 0.15331 19 
0.85 0.52777 53 0.12975 63 0.51907 53 0.13534 76 0.50960 37 0.14070 42 0.49930 03 0. 14576 60 
0.90 0.52419 89 0.124384 O1 0.51564 68 0.12945 57 0.50635 98 0.134380 07 0.49628 50 0 13881 27 
0.95 0.52108 64 0.11980 23 0.51269 40 0.12398 85 0.50360 26 0.12837 42 0.49376 60 0.13239 56 
1.00 0.518386 60 0.11461 00 0.51014 O08 0.11890 78 0.50125 14 0.12288 04 0.49165 77 0.12646 29 
1.50 0.50362 58 0.08126 98 0.49705 73 0.08316 52 0.49009 40 0.08464 49 0.48273 53 0.08565 25 
2.00 0.49860 38 0.06225 44 0.49331 16 0.06310 75 0.48778 00 0.06356 40 0.48202 38 0.06358 14 
2.50 0.49666 87 0.05017 60 0.49229 71 0.05051 97 0.48777 71 0.05050 72 0.48312 90 0.05010 74 
3.00 0.49591 88 0.04189 87 0.49222 20 0.04197 32 0.48843 22 0.04173 21 0.48457 11 0.04115 26 
3.50 0.49567 73 0.03590 37 0.49248 87 0.03582 86 0.48924 21 0.03547 38 0.48595 91 0 03482 23 
4.00 0.49566 74 0.03137 62 0.49287 21 0.03121 56 0.49004 15 0.03080 57 0.48719 63 0 03013 39 
5.00 0.49591 21 0.02501 48 0.49367 99 0.02477 70 0.49143 8&3 0.02433 68 0.48920 54 0 02368 73 
6.00 0.49624 77 0.02077 38 0.49439 59 0.02051 41 0.49254 72 0.02008 56 0.49071 73 0.01948 44 
7.00 0.49657 33 0.01775 17 0.49499 38 0.01749 14 0.49342 39 0.01708 72 0.49187 70 0 01653 67 
8.00 0.49686 44 0.01549 20 0.49548 90 0.01523 97 0.49412 63 0.01486 23 0.49278 82 0 01435 85 
9.00 0.49711 88 0.01373 98 0.49590 14 0.013849 88 (0.49469 83 0.01814 73 0.49352 02 0 01268 47 
10.00 0.49733 97 0.01234 22 0.49624 81 0.01211 33 0.49517 17 0.01178 56 0.49411 98 0.01135 90 

(4) JONES’s DISCUSSION OF THE THEODORSEN IE (se"™* sin(l — m)nr 1" 
‘ " ‘ n ze ) = » w= 
CIRCULATION FUNCTION FOR GENERALIZED MOTION sin ur 
sin mn , 
nwt 


The essential starting point for Jones’s discussion of 
He replaces uw by 
uw + tw and employs the form of Eq. (6). These equa- 
tions are then combined with the following formula 


the Theodorsen function is Eq. (9). 


given by Watson: 


”"(s) «=(11) 


sil 7 
1 and derives the form of 
—1 in Eq. (11) is used, 
By letting 


If un > O, Jones takes m = 
Eq. (7); but, for n.:< 0, m = 
and the result is different from Eq. (7) 


0 
0.002 
0.004 
0.006 
0.008 
0.016 
0.01: 
0.01 
0.0l¢ 
0.01 
0.026 
0.02: 
0.02 
0.02! 
0.02! 
0.03! 
0.03: 
0.03: 
0.03) 
0.03: 
0.04 
0.04 
0.04 
0.04 
0.04 
0.05 
0.05 
0.05 
0.05 
0.05 
0. 0€ 
0. 06 
0. Of 
0.06 
0. 0¢ 
0.03 


dif 
ars 


the 
fac 
for 
tic 
sh 


an 
do 


(5 


99 


92 





aw 
48 
19 


a 


ee 


OU 











TABLES OF THE 


THEODORSEN CIRCULATION F 


UNLTION +81 


TABLE 3 


Table of the Theodorsen Circulation Function for Harmonic Motion 


— G(k) k F(k) — G(k) 

0.072 0.87269 76 0.15411 63 
0.074 0.86958 41 0.15578 31 
0.076 0.86650 19 0.15738 48 
0.078 0.86345 12 0.15892 36 


k F(k) 


0 1.0 0 

002 0.99670 96 0.01257 97 
0.004 0.99325 79 0.02225 66 
006 ©60.98970 82 0.03077 52 






0 
0.008 0.98608 97 0.03851 30 0.080 0.86043 18 0.16040 21 
0.010 0.98242 15 0.04565 21 0.082 0.85744 35 0.16182 23 
0.012 0.97871 74 0.05230 15 0.084 0.85448 62 0.16318 64 
0.014 0.97498 80 0.05853 45 0.086 0.85155 99 0.16449 63 
0.016 0.97124 18 0.06440 43 0.088 0.84866 43 0.16575 41 
0.018 0.96748 56 0.06995 15 0.090 0.84579 92 0.16696 15 
0.020 0.96372 52 0.07520 79 0.092 0.84296 44 0.16812 03 
0.022 0.95996 57 0.08019 99 0.094 0.84015 97 0.16923 22 
0.024 0.95621 11 0.08494 91 0.096 0.83738 49 0.17029 88 
0.026 0.95246 51 0.08947 42 0.098 0.83463 98 0.17132 16 
0.028 0.94873 08 0.09379 10 0.100 0.83192 41 0.17230 22 
0.030 0.94501 10 0.09791 35 0.110 0.81877 79 0.17661 95 
0.032 0.94130 79 0.10185 41 0.120 0.80632 73 0.18007 27 
0.034 0.93762 38 0.10562 38 0.130 0.79453 72 0.18279 33 
0.036 0.93396 03 0.10923 23 0.140 0.78337 15 0.18489 04 
0.038 0.93031 92 0.11268 88 0.150 0.77279 46 0.18645 56 
0.040 0.92670 18 0.11600 13 0.160 0.76277 19 0.18756 59 
0.042 0.92310 94 0.11917 72 0.170 0.75326 99 0.18828 65 
0.044 0.91954 30 0.12222 33 0.180 0.74425 70 0.18867 27 
0.046 0.91600 36 0.12514 61 0.190 0.73570 31 0.18877 18 
0.048 0.91249 21 0.12795 13 0.200 0.72757 99 0.18862 42 
0.050 0.90900 90 0.13064 44 0.210 0.71986 08 0.18826 47 
0.052 0.90555 50 0.13323 04 0.220 0.71252 11 0.18772 32 
0.054 0.90213 07 0.13571 41 0.230 0.70553 73 0.18702 55 
0.056 0.89873 65 0.13810 00 0.240 0.69888 79 0.18619 40 
0.058 0.89537 27 0.14039 21 0.250 0.69255 26 0.18524 80 
0.060 0.89203 97 0.14259 44 0.260 0.68651 25 0.18420 43 
0.062 0.88873 77 0.14471 06 0.270 0.68075 02 0.18307 76 
0.064 0.88546 68 0.14674 43 0.280 0.67524 91 0.18188 07 
0.066 0.88222 74 0.14869 86 0.290 0.66999 42 0.18062 46 
0.068 0.87901 93 0.15057 68 0.300 0.66497 11 0.17931 91 
0.070 0.87584 27 0.15238 17 0.320 0.65556 86 0.17659 29 


u — 0 from both the positive and negative sides, two 
different limits of the Theodorsen function for real 
argument are obtained. Actually, no limiting process 
is required, since the Hankel functions are regular in 


the complex plane except at zero and infinity. A prima 


facie examination shows that his results are incorrect, 


for, if not, one then concludes that the Bessel func- 
tions are not continuous, which is false. It can be 
shown that the correct branch of the Hankel function 
is taken when m = 1. By utilizing the principle of 
analytic continuation, the resulting form of the Theo- 
dorsen function is valid for all y's. 


(5) DISCUSSION OF THE TABLES OF THE THEODORSEN 
CIRCULATION FUNCTION FOR GENERALIZED MOTION 


In the sequel, the Theodorsen function is defined as 


H(z) 


9 2 >) = F(p, 0) + iG(p, 0) 12) 
H(z) + 1y‘’(z) ‘ p ( 


C(z) = 


where 


k F(k) —G(k) A F(k) G(k 
0.340 0.64694 60 0.17375 80 0.940 0.54315 33 0. 10469 96 
0.350 0.64290 06 0.17231 37 0.950 0.54249 83 0 10393 77 
0.360 0.63902 00 0.17085 75 0.960 0.54185 81 0.10318 56 
0.386 0.63171 79 0.16792 44 0.980 0.54061 97 O 10171 05 


000 0.53943 48 0.10027 29 
100 0.53421 47 0.09360 62 
200 0.52995 60 0.08770 89 
300 0.52643 67 0.08246 43 


0.400 0.62497 63 0.16498 40 
0.420 0.61873 92 0.16205 57 
0.440 0.61295 75 0.15915 43 
0.450 0.61022 37 0.15771 73 


0.460 0.60758 79 0.15629 10 1.400 0.52349 57 0.07777 59 
0.480 0.60259 21 0.15347 40 1.500 0.52101 32 0.07356 41 
0.500 0.59793 61 0.15070 95 1.600 0.51889 92 0.06976 29 
0.520 0.59358 96 0.14800 19 1.700 0.51708 45 0.06631 74 
0.540 0.58952 58 0.14535 41 1.800 0.51551 56 0.06318 17 
0.550 0.58759 22 0.14405 33 1.900 0.51415 02 0.06031 72 
0.560 0.58572 05 0.14276 82 2.000 0.51295 48 0.05769 13 
0.580 0.58215 22 0.14024 50 2.500 0.50874 40 0.04729 69 
0.600 0.57880 16 0.13778 52 3.000 0.50628 00 0 04000 39 
0.620 0.57565 12 0.13538 85 3.500 0.50471 92 0 03462 23 
0.640 0.57268 53 0.13305 44 4.000 0.50367 09 0.03049 61 
0.650 0.57126 71 9.13191 06 1.500 0.50293 41 0.02723 63 


0.660 0.56988 98 0.13078 22 5.000 0.50239 72 0 02459 86 
0.680 0.56725 18 0.12857 08 5.500 0.50199 44 0 02242 16 


0.700 0.56475 95 0.12641 89 6.000 0.50168 46 0.02059 55 
0.720 0 56240 26 0.12432 52 6.500 0.50144 14 0.01904 21 
0.740 0.56017 12 0.12228 82 7.000 0.50124 70 0.01770 50 
0.750 0.55909 99 0.12129 05 7.500 0.50108 $4 0.01654 23 
0.760 0.55805 67 0.12030 65 8.000 0.50095 96 0.01552 20 
0.780 0.55605 09 0.11837 84 8.500 0.50085 17 0.01461 96 


0.800 0.55414 66 0.11650 24 9.000 0.50076 09 0.01381 61 


0.820 0.55233 69 0.11467 70 9.500 0.50068 38 0.01309 58 
0.840 0.55061 58 0.11290 05 10.000 0.50061 78 0.01244 67 
0.850 0.54978 66 0.11203 O01 20.000 0.50015 58 0.00624 32 
0.860 0.54897 74 0.11117 14 30.000 0.50006 94 0.00416 46 
0.880 0.54741 66 0.10948 83 40.000 0.50003 90 0.00312 42 
0.900 0.54592 86 0.10784 96 50.000 0.50002 50 0 249 96 
0.920 0.54450 89 0.10625 38 100.000 0.50000 62 0.00124 99 


@ 0.50000 00 0 00000 00 


the motion is divergent, —7/2 < 6 < 0. However, if 
the motion is harmonic, 6 = O and, therefore, p = k. 
In this case, values of C(k) are well known.'! However, 
if it is desired to study an aerodynamic system for non- 
harmonic motion, it is then necessary to know the 
value of C(z) for general z. Heretofore, tables of C(z) 
have not been available. Many studies in stability 
have been hypothesized on a _ pseudo-expansion of 
C(k) about k = 0. That this can lead to serious dis- 
crepancies is the subject of a recent report by Goland." 
In other studies, the necessity of tables of C(z) has 
been obviated by using curve fits of C() for k real and 
then substituting s = pe” for k in the curve fit expres- 
sions. An example of this approach is to be found in a 
paper by Goland and Luke."* 

Since the tendency in aerodynamics is always toward 
greater accuracy and refinement, it has become apparent 
that a table of C(z) would prove most useful for future 
research. In Tables 1 and 2, C(z) is tabulated for the 
following range of p and 6: 


p: 0(0.01)0.3(0.02)0.5(0.05) 1.0(0.5) 10.0 

6: —5°(5°)30° 
For 6 = 0, the range of values is a bit more extensive 
than that given above. All the tables were computed 
using values of the Bessel functions tabulated in refer- 
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FIG. 2 - PLOT OF THE THEODORSEN FUNCTION FOR GENERALIZED 
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TABLES OF THE THEODORSEN CIRCULATION FUNCTION 183 


ences 14, 15, 16, and 17. Seven decimals were carried, 
and it is estimated that the error is < 3 X 1077. It is 
believed that the present range is sufficiently large to 
cover current needs. The range of 6 can easily be ex- 
tended using the last two references. Plots of F and G 
versus p for different 6 are shown in Figs. 1 and 2. 

For a number of applications, linear single or double 
interpolation, as the case may be, will afford sufficiently 
reliable results. To achieve greater accuracy, one can 
employ the series expansions of the Hankel functions 
in the neighborhood of the origin. Elsewhere, Taylor 
series expansions of C(z) are valid, for, if z is small, 
neglecting terms of second order and higher, 


C(z) ~ 1 + is{log (2/2) + y + u(mr/2)] (14) 
dC/dz ~ t{log (2/2) + y + i(m/2)] (15) 
where y is Euler’s constant. Thus, dC/dz does not exist 


atthe origin. It follows from Eq. (14) that 


F(p, 0) as l nasi 


p (10g ‘ +r 7) sin 6 + C + 0) cos 6 (16) 








G(p, 6) ~ 
i p T a - 
p (tog - + v) cos @ — E + 0) sin @ (17) 
Thus, lim C(z) = 1. Note further that lim |OG/dp| > 
z—>0 po 
© for all @ and that lim |OF/0p| — © for@ #0. But, 


po 

if 6 = O, then the latter limit becomes —7/2. It can 
easily be shown that the higher derivatives of F with 
respect to p do not exist at the origin. Differentiating 
Eq. (12), one easily derives that 

2(dC/dz) = C(C — 1) + 22(2C — 1) (18) 
and, therefore, calculation of all higher derivatives is 
simple. 
If zis real, 

k(dG/dk) = (2F — 1) (G+ k) (19) 

and 


k(dF/dk) = F(F — 1) — G*? — 2Gk (20) 


It can be shown that 


gor ) (21) 
ak YoN"'"1 -f- JoJi 


and, if Eq. (21) is placed in Eq. (20), then dF/dk has 
the form as given by Rauscher and Heilenday.'* This 
reference contains a typographical error in that the 
coefficient of J)J; should be positive. 


ce : 2G l 
F(F — 1) + G = ( 


Values of C(z) for z large can be calculated rapidly 
by employing the asymptotic expansions of the Hankel 


functions. By employing these expansions, it is readily 
deduced that lim C(z) = 1/2. 


3? « 
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Three-Dimensional Shock-Wave 


Reflections: 


E. MIGOTSKY?7 ano M. V. MORKOVIN? 
Unwersity of Michigan 


SUMMARY 


Reflections of steady shock waves from three-dimensional 
bodies are analyzed by applying the sweepback principle locally 
along the intersection between shock and body. It is found that 
the mathematical analyses of two important special cases 

namely, a plane shock intersecting a circular cylinder and a coni- 
cal shock impinging on a plane surface—are identical because of 
the underlying similarity of local geometries. In both cases it is 
concluded that some type of Mach reflection must occur over 
part of the body even if regular reflection is possible at the point 
where the shock first reaches the body. Furthermore, the exist- 
ence of an expansion wave originating on the lee side of the circu- 
lar cylinder is indicated. The practical significance of these 
phenomena and their possible modifications in presence of a 


boundary layer are briefly discussed. 


INTRODUCTION 


nce REFLECTION in steady two-dimensional 
supersonic flow has been studied at length both 
theoretically'~* and experimentally.*~* * No work, 
theoretical or experimental, on reflection of steady 
shock waves from three-dimensional bodies has yet 
been reported in accessible literature. The extension 
of the relatively meager results to three dimensions is 
likely to be slow; the considerable difficulties en- 
countered in the two-dimensional problem are here 
further complicated for both the experimental and 
analytical approach. Yet, there are some configu- 
rations of practical interest where the knowledge of the 
three-dimensional reflection would be useful: wing- 
body interactions, static pressure probes in proximity 
of shock waves, three-dimensional models in closed 
supersonic wind tunnels, etc. In the present paper, 
the principle of local sweepback is developed for the 
study of these reflections in absence of boundary layer. 
This principle is applied to the case of a plane shock im- 
pinging on an infinite cylinder aligned with the free 
stream and to the case of a conical shock, with its axis 
parallel to the free stream, reflected from a plane wall 
(which is equivalent to the case of intersection be- 
tween two equal, symmetrically placed conical shocks). 
In both cases it is found that, at all Mach Numbers 
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and at all skewed shock orientations, there exists a re- 
gion where regular shock reflection cannot occur; if 
the shock remained plane for this irregular region, there 
would also be a region where no reflection at all is pos- 
sible. It is planned to follow this theoretical study 
with experimental investigations in the University of 
Michigan Supersonic Wind Tunnel in the near future. 


List oF SYMBOLS 


curve along which the shock-wave surface inter 


i = 

sects the body surface 

t = tangent toi 

n = normal to the body 

c = a constant (equal to component of Mach vector 
normal to shock wave) 

q = velocity 

M = Mach Number 

0 = shock-wave angle 

6 = angle of velocity deflection through a shock wave 

y = angle of yaw of t (Fig. 2) 

w = meridian angle (Fig. 4 and 6) 

5 = plane tangent to the shock surface 

T = plane tangent to the body 

N = plane normal tot 

V = plane containing the local Mach vector and a gen- 
erator of the body surface 

x,y, 5 = right-handed rectangular coordinates (Fig. 4 and 6) 

Subscripts 
0 = region upstream of oncoming shock 
1 = downstream of oncoming shock, upstream of reflected 
shock 

2 = downstream of reflected shock 

01 = across oncoming shock 

12 = across reflected shock 

e = “effective” 

m = maximum 

n = normal 

t = tangential 


Prime denotes parallelism to unprimed element 


(1) Two-DIMENSIONAL SHOCK REFLECTIONS 


In absence of the boundary layer, the supersonic 
flow is tangent to the body upstream of the line, i, 
where the shock impinges on the wall. Regular (also 
called simple) reflection of the shock is determined by 
the condition that the flow remain tangential to the 
boundary downstream from i. Thus, the deflection 
angle 69; across the oncoming shock must be equal and 
opposite in sense to the deflection angle 6). across the 
reflected shock (Fig. 1). 

It is well known that, for a range of upstream Mach 
Numbers M, and of oncoming shock angles 6, such de- 
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Fic. 2. Initial and final flows viewed in the plane of the wall 


flections 6y, and ‘‘middle’’ Mach Numbers J/; (upstream 
M of a regularly reflected shock) occur that the maxi- 
mun deflection 6:2, across the regularly reflected shock 
is less than 69, and the condition of tangency is vio- 
lated. The boundary of this “‘forbidden’’ range is 
represented by the curve A of Fig. 7; it is also given in 
terms of the critical 69; and /) in Fig. 4.14 of reference 
2. A second, somewhat fictitious boundary also ex- 
ists; this is the curve that represents the condition for 
which the Mach Number /, behind the initial shock is 
equal to unity. If the shock could continue straight 
all the way to the wall, then this curve, denoted by B in 
Fig. 7, would represent the limit for any reflection, in- 
asmuch as the flow beyond the first shock is no longer 
supersonic. Presumably, a Mach-type reflection oc- 
curs when the parameters \/) and 6; are in the region 
A and B. It is often forgotten that steady- 
state reflection in the form of a Mach JV with straight 
shocks (which is equivalent to the “stationary’’ Mach 
Either 


between 


reflection of reference |) is then also impossible. 
the stem of the Y curves or the boundary-layer inter- 
action transforms the reflection more radically, skew- 
ing the V pattern or producing A- and X-type patterns, 


oiten accompanied by separation. At present, such 
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cases cannot be handled by analytical methods. In 
the analogous three-dimensional situation, we shall 
therefore restrict ourselves to the location of the for- 


bidden zones on the body. 


(2) Basic PRINCIPLE OF THREE-DIMENSIONAL 
REFLECTION 


The transition from two-dimensional to three- 
dimensional treatment is accomplished by first applying 
the sweepback principle to shock reflection. In Fig. 2, 
a plane skewed shock S impinges on a flat wall 7’ (paral- 
lel to the paper) along the line i yawed by an angle ¥ 
from the perpendicular to the free-stream Mp direc- 
tion. Since the component q), (parallel to i) of the 
velocity vector is tangent to both the oncoming and 
reflected shock, it remains unchanged and therefore 
tangent to the body. The actual condition of tan- 
gency needs to be imposed only upon the velocity 
component q»z, in the plane N normalto i. Hence, in 
this normal plane, one can apply directly the two- 
dimensional analysis corresponding to Fig. 1, using the 
‘effective’? Mach Numbers AJ), for the upstream Mach 
Number 1/) and measuring the ‘‘effective’’ angles 41, 


bore, ete.,in NV. Clearly, one has 


Mo. = My cos ¥ (1) 


Furthermore, dropping a plane perpendicular to the 
shock plane S through My, one forms (Fig. 3) the fol- 
lowing three triangles: one in the vertical plane V 
containing 9);; one in the horizontal (i.e., tangential) 
plane 7 containing y; and one in the normal plane NV 


containing 1.. From these triangles, it follows that 


Sill Iie COS Y sin 4) (2) 
Eliminating y between Eqs. (1) and (2), the important 


relation 


Moe Sin O1¢ = Mo sin O01 (3) 
is obtained. The two sides of Eq. (3) are merely the 
expression for the component of the Mach vector per- 
pendicular to the shock wave, which should govern the 
whole phenomenon. Therefore, Eq. (3) expresses the 
fact that, by considering the reflection in the normal 
plane N, nothing of the essential character of the phe- 
nomenon has been omitted. It also indicates that 
the basic parameter governing the reflection is the 
product My sin 0, rather than Jo and 4 separately. 

When the oncoming shock forms a continuously 
curved surface, the line ¢ is curved and has a continu- 
ously changing tangent ¢. It can be shown rigorously 
that the previous analysis still applies Jocally at any 
point on the curve of impingement if, in Fig. 2, the 
straight line ¢ is interpreted as the tangent line ¢. 
This is analogous to the application of the sweepback 
principle to the local direction of the sharp leading edge 
of a supersonic wing with variable sweepback. The 
proof rests on the fact that finite changes across the 


+n ni, 
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shocks at the reflection curve ¢ occur within an “‘infini- 
tesimal’’ distance, whereas the changes in the same 
quantites associated with a comparable distance along 
i are only “‘infinitesimal.”’ 

A similar argument applies if the curve i is a space 
curve in three dimensions. The plane of Fig. 2 now 
becomes the plane 7’ tangent to the body, defined by 
the local upstream Mach vector A/) and the tangent 
line t, since the condition of tangency behind the re- 
flected shock now demands that the resulting velocity 
vector q2 remain in 7. In Fig. 3, S represents the 
plane tangent to the shock surface at the point and N 
the plane normal to the local space direction tf. 

Eqs. (1), (2), and (8) remain valid. From them, 
the basic reason for the anomaly of the reflection (say, 
in the case of a conical shock impinging on a flat wall 
in an otherwise uniform supersonic stream) can be 
understood. The curve iis a hyperbola with a vertex 
at the point nearest the apex of the cone. At the ver- 
tex there is no yaw. Let us assume that 1/9; and 4); 
are such that at the vertex regular reflection occurs. 
As one follows the hyperbola on the wall, the tangent 
t turns continuously and the angle of yaw y increases. 
Then, according to Eqs. (1) and (2), the effective 
Mach Number J/, in the normal plane decreases, while 
the effective shock angle increases. It is clear that the 
effective deflection angles 6o;, and 6, also change and 
that the ‘“‘forbidden”’ 
41. might be entered at some point on the hyperbola. 
That this is actually always the case will be shown later. 


range of parameters MM), and 


Just what would happen beyond the critical point 
on the hyperbola an analysis of the present type can- 
not show. It should be added that, during testing of 
cones in rectangular supersonic wind tunnels, the coni- 
cal shock is not usually of constant strength, since it 
is softened by the expansion wave originating at the 
shoulder of the finite cone. Also, there is usually a 
well-developed boundary layer on the wall of the tunnel 
which may “‘accommodate”’ the needed flow deflection 
well beyond the critical point predicted by the above 
theory. 


(3) SPECIAL CASES 


Consider a uniform supersonic flow with an infinitely 
long circular cylinder oriented along the horizont -1 
axis directed upstream (Fig. 4). A plane shock wave 
S, generated, say, by an unyawed two-dimensional 
wedge, impinges upon the cylinder at the angle 0); which 
is measured, as customary, in a vertical plane V parallel 
to the xy-plane. For the regular reflection, the curve 
of impingement i, is then the elliptical intersection of 
the cylindrical surface with the shock plane S. At 
any point P on i, a plane 7 tangent to the cylinder can 
be constructed; this plane contains the free-stream 
Mach vector J/) attached to P and the line t tangent 
to i. As explained in the previous section, it is now 
possible to apply the sweepback principle locally and 
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Nominal and effective Mach Numbers and shock angles 
according to the sweepback principle. 
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a circular cylinder, 


Fic. 4. 


Plane shock wave impinging on 


to obtain Eqs. (1), (2), and (8) with the aid of a con- 
struction completely analogous to that of Fig. 3. 

It is still necessary to express the two fundamental 
“effective’’ variables, 1p, and 6;,, as functions of the 
location of the point P on the intersection curve i. In 
the present case of a uniform flow and a plane shock 
wave, the component of the Mach vector normal to 
the shock, J/o sin 01, has a constant magnitude larger 
than unity, say c. Eq. (3) then yields the following 
relation between the effective variables: 


Mo. sin O1- = € (4) 


If the meridional angle w between the normal to the 
cylinder m and the y-axis (Fig. 4) is used as the par- 
ameter defining the location of point P on the inter- 
section i, it is shown in the Appendix that the “effective” 
wave angle is given by the simple relation 


COS Joi. = COS w COS Dy (5) 


Eqs. (4) and (5) completely determine the ‘‘effective” 
variables Mo, and j¢. 
able to derive from Eq. (5) the expression for csc Box 


For later purposes, it is desir- 


csc O16 = csc Oo,/V 1 + sin? w cot? 0; (6) 


It is shown in the Appendix that Eqs. (4), (5), and (6) 
remain valid for the case of a conical shock, with axis 
parallel to the uniform stream, impinging upon a plane 
wall, also parallel to the free stream (Fig. 6). The 
reasons for this duality between the reflection of a 
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THREE-DIMENSIONAL 


plane shock from a cylinder and the reflection of a 
conical shock from a plane are also discussed in the 
Appendix. 

(4) DISCUSSION OF RESULTS 


In Section (2) it was shown that the principles of 
two-dimensional shock reflection can be extended to 
three-dimensional reflection if the ‘‘effective’’ variables 
My, and 4. in the plane normal to the intersection 
line i are used rather than the true free-stream Mach 
Number ./o and true shock angle 4). In particular, the 
boundary of the ‘‘forbidden”’ range, where regular shock 
reflection becomes impossible, is again the curve A 
of Fig. 7 if AZo, and cse 4, are the coordinates. Simi- 
larly, the curve B of Fig. 7 bounds the region of the 
My, and ese 61. variables for which the normal com- 
ponent of the velocity behind the oncoming shock wave 
is subsonic. In Section (3) and in the Appendix, the 
particular dependence of these basic variables on the 
parameter w, specifying the position on the intersection 
curve i, is developed in the two special cases considered. 
For these cases, it is now possible to examine what 
type of reflection is to be expected at any point on the 
body. 

The criterion for the existence of the regular reflection 
is that the geometry of the oncoming plane shock wave 
and of the body produces a combination of the effective 
variables Mo, and csc 61. which would fall above the 
curve A in Fig. 7. For values of MJ, and esc 6), fall- 
ing between curves A and B, a Mach type reflection is 
indicated; for the values of 1/9, and csc 6), below curve 
B, no reflection would be possible if oncoming shock 
were to remain plane. In the case of shocks with con- 
stant strength, it was seen that the possible combina- 
tions of 7p, and 4, are governed by Eq. (4), which can 
be put into the more useful form of 


csc 09,. = (l/c) Mo (7) 


The advantage of using ese 4, and AJ), for the co- 
ordinates in Fig. 7 is seen from the fact that, for a 
given shock strength characterized by c, the possible 
combinations of these “‘effective’’ variables lie along a 
straight line. The intercept of this line on the ordi- 
nate csc 1. = 1 is equal to c, the component of the 
Mach vector normal to the shock wave, 
If the oncoming 


free-stream 
which is therefore larger than unity. 
shock strength is increased monotonically from that of 
a Mach wave (c = 1), the straight-line locus of the 
possible combination of the ‘‘effective’’ variables shifts 
to the right and becomes increasingly flatter.* = In 
all cases of shocks with finite strength, therefore, these 
loci pass through both the zones of ‘‘no reflection”’ and 
“Mach-type reflection.” 

In order to verify that the points on the loci in these 
two zones correspond to actual points on the reflecting 


*The loci form a one-parameter family of straight lines 
through (0, 0). The significance of the single parameter was 


pointed out in comments under Eq. (3) 
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Fic. 5. Geometry leading to Eq. (5). 














Fic. 6. Conical shock impinging on a plane surface 


body, consider Eq. (5), which relates the effective wave 
angle 4:, with the position of the point P on the curve 
of intersection i. In both of the special geometries con- 
sidered, let the meridional angle w vary continuously 
from 0° to 90°. Then, by Eq. (6), ese 41, will decrease 
continuously and monotonically from cse 4, to unity, 
the value corresponding to the intercept c. From Fig. 
7 it is therefore seen that, even if regular reflection is 
possible for w = 0, there always exists a set of points 
on the intersection curve t (corresponding tow, <w S 
90°) for which regular reflection is not possible. Fur- 
thermore, there will also exist a region (w, S w < 90°) 
in which no reflection would be possible if the shock 
remained plane. 

As an example, consider the following configuration : 
M, = 1.90 and @; = 40° (c = 1.222). For this case, 
the variation of Mp, and 4%, with w is plotted in Fig. 8 
here. If 


(solid curves). The critical value w, is 23.8 
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the free-stream Mach Number is increased to 2.5 and 
6, is taken to be 29.2° (in order to keep the strength of 
the shock and c¢ constant), the critical angle w, is in- 
creased to 36.2°; however, the region of nonregular 
reflection still exists. 

Whereas the local shock reflection phenomena for 
0° <w < 90° are identical for the two special geometries 
considered, there does exist a region in the case of the 
plane shock reflecting from the cylinder which has no 
counterpart in the conical shock reflecting from the 
plane wall. This is the region on the lee side of the 
cylinder (90° < w < 270°). In this region, the bound- 
ary condition at the (here fictitious) line of intersec- 
tion ¢, which requires the flow to be tangent to the 
body, is of opposite sense to that of the regions 0° < 
w < 90° and 270° < w < 360°. Thus, instead of 
a reflected shock, an expansion wave is required to 
turn the flow back toward the cylinder.* If 
visualized the expansion wave as originating immedi- 
ately downstream from the intersection line i, the Mach 
lines of this wave, corresponding to a reduced Mach 
Number, would intersect the original shock from be- 
hind and thus reduce its strength in the vicinity of the 
In fact, the strength of the expansion wave at 
180° is sufficient to re- 

As a result, there is, 


one 


body. 
the point corresponding to w = 
duce the shock to zero strength. 
on the lee side of the body, a ‘“‘shadow”’ that weakens 
the shock. 

In the region 0 S w < w, where regular reflection 
occurs, the analysis of the shock reflection in the nor- 
mal plane N leads to the initial direction of the reflected 
wave at the line of intersection. However, it is im- 
possible to predict the shape of the reflected shock 
away from the body without a rather difficult analysis 
of the complete three-dimensional flow. Still, an inter- 
esting application of the information on the initial di- 
rection of the reflected shock can be made in supersonic 
wind-tunnel calibration. Plane shock waves origi- 
nating at the sidewalls of a tunnel are not visible with a 
schlieren beam passing through the side windows, but 
the introduction of a cylindrical body into the channel 
renders the reflections off the cylinder visible and there- 
by makes detection of sidewall shocks possible. 

These considerations have assumed the fluid to be 
nonviscous. Actually, a boundary layer will always 
cover the reflecting surface. One effect of the presence 
of the boundary layer has been observed on schlieren 
pictures of bodies of revolution supported by a vertical 
wedge from the floor of supersonic tunnels. As the 
shock wave generated by the support cuts across the 
cylindrical body, the expansion wave on the lee side 
does not seem to build up to its theoretical strength, 


* Here, it is useful to visualize an existing plane shock pierced 
the complete field, which 


” 


by the cylinder and call ‘‘reflection 
for the distortion of the shock from its plane 
The point of view of refraction of a 


would account 
form and constant strength. 
shock wave around the cylinder would probably be more fruitful 


if linearized methods of analysis were used. 


AERONAUTICAL 


SCIENCES—JU 






































Tooe| 
att a4 | 
2.0/2 | aZexe }— 
1.8}—— 
= 
——; 
2 | | 
1.4 [— + : —— +—— 
| Li 20 
i aes 4 aii 
| alee 
| | LAT [— 8 
' Lad | 








) 
1.0 12 1.4 1.6 16 6=6200C 22 
EFFECTIVE MACH NUMBER, M,, 



























































Fic. 7. Representation of the possible shock reflections 
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Variation of reflection conditions along the body for two 
special cases. 


Fic. 8 


and a sharp thickening of the boundary layer occurs. 
In some cases a deflection of 2° to 3° 
dated by the boundary layer. Transition or separa- 
tion of the boundary layer on bodies of revolution 
tested in tunnels has sometimes been ascribed to the 


is thus accommo- 


impingement upon the body of the nose shock re- 
flected from the walls. The preceding discussion indi- 
cates that the region on the lee side of the body is eriti- 
cal in this respect, perhaps more so than the region 
where the shock strikes the body. The implications of 
this tendency toward separation for wing-body con- 
figurations are clear. 

In the region where Mach-type reflection is indicated, 
the boundary layer may thicken and even separate. 


This is the region where strongest anomalies are ex- 


pected because the reflected shock should be relatively 
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strong. It is difficult to visualize a discontinuous 
change in the vicinity of w,. 

Clearly, the location of the critical regions on the 
body, #, <0 < 90°, makes it difficult to observe the ac- 
tual phenomenon by optical methods. This un- 
doubtedly explains why no similar anomalies have been 
seen thus far. Experiments designed to clarify the 
behavior of the shock and boundary layer in some of 
these regions are scheduled in the 8- by 13-in. supersonic 
tunnel at the University of Michigan. 

It is paradoxical that investigations of any shock, 
and in particular of the shock reflected from a body, by 
probing with a thin static pressure ‘“‘needle’’ from be- 
hind are likely to be inaccurate because of the same 
phenomenon; the shock is impinging on the cylindrical 
body of the probe. It would seem desirable to place 
the orifices at meridional positions 90° from that at 
which the shock first intersects the probe; at these 
points the pressure is least likely to be disturbed. 


(5) CONCLUDING REMARKS 


The sweepback principle applied locally at points 
along the line of intersection of a shock-wave surface 
and a three-dimensional body makes it possible to 
analyze what type of reflection should take place in 
absence of the boundary layer. For the important spe- 
cial cases of a plane shock impinging on a cylindrical 
body and a conical shock impinging on a plane wall, 
the mathematical analyses are identical because of the 
underlying similarity of local geometries. In both 
cases it is concluded that, even if regular reflection 
is possible at the points on the zero meridian w = 0 
where the shock first impinges upon the body, there 
will always be a region w, < w < w, where regular re- 
flection must be replaced by some type of a Mach re- 
flection and a region w, < w < 90° where no reflec- 
tion would be possible if the shock remained plane. 
Furthermore, the existence of an expansion wave 
originating along the line of shock intersection on the 
lee side of the circular cylinder is indicated. These 
phenomena, whether appreciably modified by the 
presence of the boundary layer or not, may have prac- 
tical implications for wing-body interactions, testing of 
three-dimensional models in rectangular supersonic 


tunnels, pressure probing near shock waves, etc. 
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APPENDIX 


In general, relations of the type of Eq. (5) can be 
derived most easily with the aid of vector analysis; 
in the present case it is obtained from the following 
geometrical construction shown in Fig. 5: The verti- 
cal plane V and the normal plane N through the point 
P as seen in Fig. 3 are extended beyond P into the 
body of the cylinder if necessary. The yz plane is also 
made to pass through P. Finally, a plane 7” parallel 
to the tangent plane 7 is constructed so that it forms 
a tetrahedron PABC with the three planes V, N, and 
yz. Theedge AC between 7” and V is parallel to the 
free-stream Mach vector and the generators of the 
cylinder; AC makes, therefore, a right angle with all 
the lines in the ys plane, in particular with CB and CP 
(which is parallel to the y axis). The edge PB is identi- 
fied with n, the normal to the cylinder at P—as a line 
in the normal plane JN, it is perpendicular to the tan- 
gent line t; as a line in yz plane, it is perpendicular to 
the generators of the cylinder. Hence, nm (PB) makes 
a right angle with the edges BA and BC located in the 
7” plane. Furthermore, the angle between nm and 
PC in the yz plane is equal to w by definition. Asan 
intersection between the vertical plane V and the nor- 
mal plane N, the edge AP is normal to the shock plane 
S and to the lines DP and EP of Fig. 3 used in defining 
the shock angles* @ and @,.. Hence, by plane geometry, 
the angle between AP and CP in the vertical plane 
V is equal to 0, and the angle between AP and BP in 
the normal plane N is equal to @,. 

Having identified the three angles at P of the tetra- 
hedron as w, 0, and @,, one obtains the required relation 
(5) from the identity of ratios 


BP/AP = (BP/CP) (CP/AP) 


cos #, = cos w cos 6 


The case of a conical shock, with axis parallel to the 


uniform supersonic stream, impinging upon a plane 


* Since the Appendix is concerned with the oncoming shock 
wave only, the subscript 01, which differentiates between it and 
the reflected shock, is dropped everywhere 


(Continued on page 504) 
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The Laminar-Turbulent Transition in a 
Boundary Layer’—Part I[' 


H. W. EMMONS} 


Harvard Unwersity 


SUMMARY 


A theory of transition from laminar to turbulent flow is de 
veloped from some new observations of the transition phenomena. 
Every point of a body is turbulent a fraction of the time, and a 
formula is derived for this fraction in a phenomenological man 
ner. Application of the theory to available flat plate data gives 
some insight into observed effects and shows the way to addi- 


tional experimental results required. 


INTRODUCTION 


Shes BOUNDARY LAYER that grows on the surface 
of any body is at first laminar. Any fluctuations 
with time which are carried in from the free-stream 
turbulence, which are carried in by the passage of 
sound waves, or which are produced by roughness on 
the body or by vibrations of the body are generally 
damped rapidly near the leading edge. 

As one proceeds downstream from the leading edge, 
confused time-dependent velocities arise which con- 
sume the laminar boundary layer and replace it by 
a more rapidly growing and thicker turbulent one. 
Data on the overall drag of bodies clearly show the 
characteristic values of laminar drag at low Reynolds 
values of turbulent 

In between there 


Numbers and the characteristic 
drag at high Revnolds Numbers. 
is a more or less extensive region of transition from 
one to the other. 

The laminar-flow drag and heat-transfer coefficients 
agree well with theoretical results where the calcula- 
tions have been carried out. The turbulent flow, while 
too complex to permit a basic theory as yet, gives ex- 
perimental drag and heat-transfer coefficients that 
correlate well against the appropriate dimensionless 
variables, Revnolds Number and Prandtl Number. 
Transition data have not been treatable in either way, 
so that at present no precision of prediction is possible 
if transition conditions are encountered. 

It has appeared empirically that there is a critical 
value of Reynolds Number at which transition starts. 
Fig. | shows the data of Schubauer and Skramstad! 
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* The author is indebted to A. E 
cussions and assistance with the numerical work, to Dr 
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integral equation, and to the Hughes Aircraft Company for pro 
viding the stimulating atmosphere in which this fundamental 
investigation was carried out 

+ Part II will be published in the Proceedings of the First U.S 
National Congress of Theoretical and Applied Mechanics 
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and of Hall and Hislop’ for a flat plate lined up care- 
fully with the air in a zero pressure gradient wind 
tunnel. The apparent agreement is belied by the 
authors’ statements that the Reynolds Numbers given 
were the ‘‘approximate beginning and end of trangi- 
tion’ and ‘‘the mean transition position’’ in the two 
references, respectively. The unsteady nature of tran- 
sition prevents more precise work on critical Reynolds 
Numbers. 

It has long been known that the theoretically derived 
formulas are correct solutions for all 
Thus, turbulent flow appears as 


laminar-flow 
Reynolds Numbers. 
a more stable, though nonsteady, type of flow which re- 
places the laminar flow. Considerations of stability of 
the laminar boundary layer by many authors (see, for 
example, Lees*) show that the damping of infinitesimal 
disturbances in the laminar boundary layer depends 
upon the Reynolds Number, the disturbance frequency, 
and the rate of heat transfer. A typical curve is given 


as Fig. 2. The amplified frequencies appear as meas- 


ee 


f 
' 


Spe nme sees mR 


urable Tollmein-Schlichting waves in the boundary | 


layer as shown by Schubauer and Skramstad.! 
agreement between computed and experimentally deter- 
mined neutral curves is fair. The differences may well 
be caused by the fact that the experimental disturb 
ances were finite. 

No theory has yet been attempted on the effects oi 
finite disturbances on a boundary layer. The non- 
linear nature of this problem has, and will for some 
time, continue to block further progress in this direc- 
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Schlichting’s laminar boundary-layer curve (solid) 
Experimental 


Fic. 2. 
Frequency of disturbances vs. Reynolds Number. 
points (dashed) by Schubauer and Skramstad. 


tion. On the other hand, correlation of transition data 
has not been successful, as Fig. 3 (taken from Schoen- 
herr?) shows. von Karman® has stated in this con- 
nection: ‘‘However, in view of the Schoenherr curves, 
it must be remarked that the transition phenomenon 
is probably too unstable to be represented by a simple 


equation.”’ 


A NEw OBSERVATION 


During the past year, a simple water-table analogy 
to supersonic flow was completed at Harvard as a stu- 
dent demonstration. On this table, water flows in 
a '/s- to '/4-in. deep layer over a nearly horizontal 
piece of plate glass about 3 by 6 ft. The velocity and 
volume flow of the water are controlled by the adjust- 
ment of the tilt of the glass plate and the adjustment 
of a by-pass to the pump. It was observed that, in 
addition to the anticipated ‘‘supersonic’’ phenomena, 
laminar to turbulent transition was obtained. 

The mechanics of transition could be clearly ob- 
served, since the turbulent motion of the thin layer of 
water disturbed the water surface and, hence, altered 
the appearance of the surface in both reflected and 
Fig. 4 shows a plane view of the 
At some in- 


transmitted light. 
glass plate with water flowing over it. 
stant and at some point (both random), a tiny spot 


of turbulence appeared as at S. The spot was small 


and of irregular shape. This spot grew as it was 
Washed along with the boundary-layer fluid. The 


approximate shape and size of the spot after a period of 
growth is shown at A in Fig. 4, as is also the locus of 
turbulent 

At random times, 


successive spot positions. Each spot was 


washed off the end of the plate. 


TRANSITION 
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spots appeared at random locations over the plate sur- 
face. These tiny turbulent origins will be called ‘‘tur- 
bulent sources’’ (more briefly ‘‘sources’’). Sources 
could be increased in density and frequency at a given 
location on the plate either by increasing the water 
velocity or adding disturbances by placing an object 
(finger) in the deep (low-velocity) water well ahead of 
the plate. A disturbance (as a vertical wire) placed 
directly into the water on the plate acted as a con- 
tinuous source, as at B of Fig. 4, and thus filled a sector 
of the plate surface with turbulent water. 

The growing spot approximately preserved its shape 
as it grew, and the locus of successive spot positions 
appeared straight and at the same angle as the con- 
tinuous turbulence from a wake. These facts suggest 
that the oriented vorticity of the laminar boundary 
layer is confused at the edge of the spot by the velocities 
induced (in the oriented flow) by the already confused 
vorticity of the turbulent region. Thus, the spot 
grows at each point normal to its surface by consuming 
the laminar region. The velocity of turbulence prop- 
agation would be expected to be a function of the angle 
between the spot normal and the local direction and 
strength of the oriented vorticity. (The noncircular 
spot shape shows the direction of the normal to be im- 
portant.) The initial vorticity is important because a 
turbulent region spreads through a boundary layer 
much faster than into the potential flow region outside 
of the boundary layer. 

The turbulence in transition is thus seen to be es- 
sentially intermittent. Each point of the plate is 
covered by a laminar boundary layer, except during 
those periods of time during which a turbulent spot 
moves over it. Since spots grow as they move, the 
downstream points of the plate are covered for greater 
fractions of the time by a turbulent layer. In fact, the 
spots soon overlap and cover downstream portions con- 
tinuously. 

No quantitative measurements have yet been made 
on the water table, so that other evidence must be 
used at present for evaluating the spot characteristics. 
Although the above-detailed process has not been de- 
scribed before, there is abundant evidence that the 
same process occurs generally. 

Charters® discovered that turbulence in a boundary 
layer on a flat plate in a wind tunnel grew over the 
plate from the edges where the plate was supported 
at the tunnel The between the 
laminar and turbulent region was mapped out by hot 


walls. boundary 
wire exploration and was found to be a straight line 
making an angle with the free stream of about 9.6 
Such straight-sided turbulent areas could also be pro 
duced by allowing the (eddying) wake from a body to 
strike the plate boundary layer. These effects de 
scribe precisely the continuous sources mentioned above. 
Schubauer and Skramstad,' while studying the growth 
of the Tollmien-Schlicting waves in the laminar bound 


ary layer, observed that, at a sufficient distance from the 
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plate leading edge, the hot-wire record showed more 
or less frequent bursts of turbulence interrupting the 
otherwise relatively low-frequency smooth wave pat- 
tern. These bursts would, by present theory, be inter- 
preted as the passage of turbulent spots over the wire. 
Charters (personal communication) has recently ob- 
tained some spark photographs of a projectile in flight 
on a range at Aberdeen Proving Ground which show a 
boundary layer partly turbulent and partly laminar 
(see Fig. 5). The turbulent spots were different at 
each spark station along the flight path. 
to be good evidence that the above-described transition 


This appears 


process occurs in a boundary layer in air at super- 
sonic speed. The fact that the boundary layer on 
some projectiles is wholly turbulent can probably be 
ascribed to fixed disturbances in the projectile geom- 
etry. 
The 
therefore to be visualized as follows: 
a laminar boundary layer completely covering the given 
This 


boundary layer is disturbed (in space, time, frequency, 


transition from laminar to turbulent flow is 


Viscosity builds 
body (this may include regions of separation). 


and amplitude) by random motions carried in by the 
fluid from the free stream, carried in through the 
fluid as sound waves, produced in the boundary layer 
by surface irregularities, or produced by the vibra- 
Each of these disturbances ampli- 
The sum 


tions of the plate. 
fies or damps as it moves along the surface. 
total of these disturbances is to be visualized as dis- 
turbing the calm of the laminar boundary layer in the 


same way as random waves disturb the calm of the 
sea. 

At each point of the laminar boundary layer there 
is assumed to be a certain critical (finite) amplitude 
(and perhaps critical frequency range) which, when 
exceeded by a disturbance, is capable of causing the 
oriented vorticity of the laminar flow to be replaced 
the confused motion of turbulence, just as _ the 
As the 


ocean waves approach the shore, their amplitudes (and 


by 
waves of the ocean may break into white caps. 


frequency) change, so that at random points the wave 
tips break. 
the amplitude and frequency of disturbances and the 
critical conditions change with the distance from the 
leading edge. Every now and then at various points, 
the boundary-layer disturbances ‘‘break’’ (whatever 
this may mean), and the flow becomes locally turbulent 
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Similarly, in the laminar boundary layer | 


————— 





the 
flu 
poi 
col 
tht 
by 

tii 


I 
tral 
tiol 
spe 
Eac 
que 
scri 
tha 
we 
(ori 
tior 
of S 

1 
zs 
tim 
port 
that 
mo\ 
spac 
ent 
grol 
(doe 
the 
be 
rela 
fron 
“pre 
surf: 
rays 

\ 
that 
turb 
x) a 
fron 
the ; 
of P 
and 
upst 
volu 
P is 
lami 


toa 


ther 


of the 


there 
litude 
when 
g the 
laced 
s the 
is the 


(and | 


wave 


layer | 


d the 
n the 
ints, 
tever 
ulent 


ae | 








LAMINARYFURBULENT TRANS 


a turbulent spot has been born. From this spot, 
the confusion grows in all directions (relative to the 
fluid), carrying all before it. Thus, the flow at any 
point on the body will be laminar part of the time and 
completely turbulent for the remainder. Transition is 
thus a random phenomenon, which can be 
by a probability function specifying the fraction of the 
time that the flow at eath point is turbulent. 


described 


PROBABILITY TRANSIFION THEORY 


In order to developa quantite tava formulation of the 
transition process, we agsumag the existence of a func- 
tion g(x, y, ¢) of position (on the surface) and time which 
specifies the rate of spot ‘production per unit area. 
Each spot so produced must be followed in its subse- 
quent history. These spots groweanel overlap as de- 
scribed before. Thus, if the fraction of the time f(x, y) 
that a given point P(x, y) is turbulent is to bécomputed, 
we must add up the times that it tsovered. by spots 
(originating upstream somewhere). In this summa- 
tion, however, we must not count twice er 
of spots which overlap. 

To do this, consider an x, y, ¢ space, Figi' 6; ~ we 
x, y are coordinates on the cae surface and! #\is ws 
time—the coordinate x to be taken along spot (trans 
port stream lines. A source at Py produces a ‘Spot 
that, on the body surface (projection in the x, y plate} 
moves and grows as described above. In the x, y, ? 
space, the spot sweeps out a volume of a shape depend- 
ent upon the manner and rate of spot movement and 
Since the spot is swept away by the stream 
did, 


growth. 
(does not propagate ahead of its source—if it 
the spot would be self-sustaining and the body would 
covered by turbulence) and 
the volume will always diverge 
We shall call this the 
its 


be permanently grows 


relative to the fluid, 
from Pp» in a conelike manner. 
“propagation cone’ for a source at Po(xo, Vo, to) 
surface elements are propagation rays of the spot (the 
rays may or may not be straight). 

We are interested in finding the fraction of the time 
that a given point P(x, y, t) is turbulent. P will be 
turbulent if source has occurred upstream (x) < 
x) at such a time that P lies in the propagation cone 
from Py. The locus of all points Py which can influence 
the state of turbulence at point P is a volume upstream 
of P as shown in Fig. 7. This volume is also conelike 
and has generators that are propagation rays drawn 
“the dependence 
then 


This will be called 
If there is only one source, 
otherwise, P is 


upstream from P 

RK for P. 
P is turbulent if that source is in R; 
Thus, the probability of turbulence at P due 


volume”’ 


laminar. 
toa source at Po is 


flit Poisin R 


AP. Po) = % ; ; (1) 
. (0 if Py is not in R 


In a volume element dV» (= dxo dyp dty) at point Po, 


there will be g(Py) dy sources, each of which will cause 
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showing turbulent 


Spark photograph of projectile 
Aberdeen Proving 


(by A. C. Charters, 


Fic. 5. 
boundary-layer spot 
Ground). 
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Fic. 6. Propagation of a turbulent spot in x, y, ¢ space 
P to be turbulent with probability y(P, Po). We can- 


not find the fraction f(P) of the time that P is turbulent 
by integrating ¢(P, Po)g(Po) dV however, since there 
may be more than one source in R and the resultant 
turbulence at P would thus be counted twice. To 
avoid these overlaps, we shall assume that the turbu- 
lence produced nearest the leading edge of the body will 
be counted as causing the turbulence in the spot over 


lap regions. There is some physical justification for 
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this assumption in that the turbulent boundary layer 
grows faster than the laminar one and, hence, the spot 
that has traveled furthest will generally be thicker: 
Until many of the physical parameters are better 
known, this assumption makes little (sometimes no) 
difference as compared with other plausible assumptions. 

Let ¥(P, Po) dV» be the probability that the flow at 
P be turbulent because of sources in dV) at Py but not 
because of any source at Po!(%o!, yo', to!) with xo! < 
xo—i.e., the source at Py is to be counted only if it is 
the point nearest the leading edge that is affecting P. 
The probability ¥(P, Po) dV» is connected with g(Po) by 


VW(P, Po) dVo = o(P, Po) g(Po) dVo X 


E — f V(P, Put) aVa'| (2) 
Ri 


where R' is the portion of the influence volume R for P 
which lies between the leading-edge plane x = O and 
the plane through x = x». The factor 


E 7 Y(P, P,') ave'] 
Ri 


is the fraction of the time that P is laminar so far as 
sources ahead of xp are concerned and, hence, is the 
fraction of the time during which a source in dV) might 
be effective in producing turbulence at P. Eg. (2) 
is to be regarded as an integral equation for ¥(P, Po). 

After Y(P, Po) is found, the probability that P be tur. 
bulent—the fraction f(P) of the time during whicly P 
is turbulent—is 


S(P) = I V(P, Po) dVo (3) 
R 


The general solution of the integral equation is car- 
ried out in Appendix (I), and the fracéion of the time 
during which a point P is turbulent is shown to be 


f(P) = 1 we Zager (4) 


This can now be workéd One Mumefically as soon as 
the function g(P») is determinéd and the influence vol- 
ume R is defined#~Before considering these ques- 
tions, one other general_-qaestion should be examined. 
The record of the turbulenee at a given point on the 
body would show turbulent’ and laminar periods fol- 
The frequency and size of the 
Spots from various 


lowing one another. 
turbulent bursts can be computed. 
sources pass over P, without necessarily being sepa- 
rated by laminar regions. We shall call the number 
of sources that contribute to the 
“the spot frequency N(P),” whether or not the spots 
are separated by laminar flow. The number of dis- 
tinct spots would be somewhat less than this. 

To find the spot frequency V(P), observe that during 
a time 7, point P will be turbulent for 7rW(P, Po) dV 
Each spot from a 


turbulence at P 


sec. because of sources in dV». 
source at Py, will cause P to be turbulent for At(P, Py) 
sec., where Af is the yp = constant chord of the S(P, xo) 


cross section of R as shown in Fig. 7. Hence, the num 
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ber of spots passing P from sources in d Vy as Py in 7 sec. 
is 
TW(P, Po) dVo/ At (5 
The spots per second passing P from dV at Po is 
dN = ¥(P, Pe) dVo/ At (6 
and the total number of spots per second passing P is 
“YW(P, Po) dVo _ 
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The average duration of a turbulent spot is then 


T = f/N (8) 


TRANSITION THEORY FOR A FLAT PLATE 


To make application of the preceding theory to transi- 
tion on a flat plate, it is necessary to know the source 
rate density g(x, y, t) and the propagation constants 
for the resultant turbulent spot. Until further ex- 
perimental work has been done (a theoretical approach 
to source production appears difficult), no definite in- 
formation on g is available. As a first and simplest 
assumption, g will be taken as independent of position 
on the plate. For a steady mean flow, g will be inde- 
pendent of the time. 

Again, for the spot propagation constants more ex- 
perimental work is needed before precise values can be 
known. (A theoretical approach is probably feasible, 
although it has not yet been attempted.) We shall 
assume that each point of the spot boundary propagates 
at a constant rate away from the source. Thus, the 
propagation cone of Fig. 6 is a true cone with straight 
generators. It is not a circular cone, however.) The 
dependence volume R is also a true (noncircular) cone. 

Therefore, the integral of importance is 


ficav=ev (9) 
R 


where |" 
x. The volume of a cone is one-third the area of the 


base A, times the altitude 


V = Ayxe/3 = Ayr'/3 (10) 


cone cross section at unit distance from 
Therefore, 


where A; = 
the apex = area of S(x, x;) atx — x = 1. 


7 


! o 
dy, dty _ : 
\% — Xy)° ./ * constant section of cone l 


(11) 


The x = constant cross section of the dependence 
cone may be considered the ‘‘time-width shape” of a 


spot. It could be measured by noting the time history 
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is the volume of the dependence region of 
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of a spot as it passes over a line at a fixed distance x from 
the leading edge of the plate. Thus, the parameter 
o is a dimensionless propagation parameter of the spot. 

The area of the x = constant cross section of the cone 
can be obtained from the ¢ = constant cross section 
(spot) as shown in Appendix (Il). Thus, the dimen- 
sionless x = constant cross-section area of the propaga- 
tion cone is given in terms of the turbulent spot dimen- 


sions by 
o = a*d/B (12) 
where 
a = tan (one-half propagation angle) = 0.17 [see 


Fig. 4 and Eq. (19) ] 
velocity of center of spot 
‘free-stream velocity U 

area of spot 


square of half width 
The fraction of time that x is turbulent thus becomes 


f(x) =1l1-—e (13) 
which is plotted as Fig. 8. Although g was assumed 
independent of x in the above derivation, it does not 
have to be considered as starting from the leading 
edge of the plate—1i.e., x in the above formulas could, 
if desired, be measured from a critical value x, given 
by Re, = Ux,/v = constant. However, in view of the 
slow start of f(x), there is insufficient data now available 
to show this to be necessary. Furthermore, the trans- 
verse contamination of Charters® shows that the bound- 
ary layer right from the leading edge is unstable rela- 
tive to a sufficiently violent disturbance of which the 
turbulent boundary layer along another body (the 
tunnel wall) is a sufficient sample. 

The remaining turbulence characteristics for a flat 
plate with a constant spot rate density now follow from 
Eqs. (A10), (A12), (A13), (7), and (8). 


Pie, &) = 1 — ¢ SAO - Ge (14) 
og . — .) 7 
¢(x, Xo) = v (x — x9)2e7 (#/8U) bs oe ooo 
- - _— » - (og/3L’) [x (x x9) *] . 
V(x, Xo) = v(x, Xo) ge (16) 

ra 
N(x) = 2age~7***/%! | te ee /3U de (17) 

0 

ers su _ 1 

T(x) = (18) 


"x 
‘ te x8 /3L te 
2ag ie dk 
we 


Of these, the number of turbulent spots passing a 
poiut x is plotted as in Fig. 9. 

The values of the ‘“‘constants’’ o and g can only be 
crudely estimated from currently available data. The 
most important constant g can best be estimated 
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by fitting the present theory to the small amount of 
transition data available. 

The spot characteristic « can, at present, be esti- 
mated as follows: The spread of the permanent turbu- 
lence of transverse contamination was measured by 
Charters. °® Thus, for 
present purposes, 


The angle was about 9.6°. 


a = tan! 9.6 = 0.17 (19) 


The value of 6 can be estimated from the recent 
spark photograph, Fig. 5. A small shock wave is 
visible at the edge of the turbulent spot and is at an 
angle different from the Mach waves in this region. 
Since a spot is a (supersonically) small disturbance, 
the shock angle is evidence of the spot velocity relative 


to the projectile. From this photo, 


B = 0.43 (20) 


that is, the spot travels at about one-half of the free- 
stream velocity. 

The value of \ which depends approximately upon 
the spot area can be estimated at present only by the 
spot appearance as visually observed. The spot looks 
about as shown in Fig. 4 and, hence, 


A = 2 (21) 


Thus, 

¢ = a’\/B = 0.134 = 0.1 (22) 
In view of the crudeness of the approximations, o will 
be taken as 0.1. 

We note that o above was evaluated from some low 
subsonic boundary-layer tests, some supersonic tests, 
and water flow in a sheet (not a real boundary layer). 
It is an assumption that must be eliminated by further 
tests or theory that the spot propagation constants are 
the same for such diverse conditions. 


TRANSITION DRAG ON A FLAT PLATE 


The drag on a flat plate arising from the laminar 
boundary layer near the leading edge is given by the 
theoretical formula 


Cy, = 7/(1/2)pU? = 0.664/Re (23) 
The corresponding mean friction coefficient is 


Cp, = D/(1/2)ApU? = 1.328/Re” (24) 


If the boundary layer is turbulent, the local friction 
coefficient is given by the semiempirical formula 


1/V Cy, = 1.7 + 4.15 logy (ReC,,) (25) 
The corresponding mean friction coefficient is given by 
] \ Cr — 1.15 log io (ReCp,) (26) 


Using a C,, based upon the distance from the leading 
edge, the mean local coefficient is 
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Fig. 10 shows these coefficients for several values of 
the transition parameter ogy’, U*. Measurements of 
the surface velocity gradient at various points along a 
glass plate by Burgers and van der Hegge Zijnen’ 
have been converted to local friction coefficients and 
are compared with the theoretical curves. The rea 
son for the poor agreement in both laminar and turbu 
lent regions has not been investigated. More local 
coefficient data are badly needed. 

In Fig. 10 we note that, although no critical Rey 
nolds Number has been assumed, the different values oi 
ogv®/U* lead to different curves, which could well 
have been interpreted as implying a critical Reynolds 
Number. We further note that, as the spot rate den 
sity g is increased, the transition region moves forward 
on the plate and becomes narrower. 

If we define “the transition point” x, as the point on 
the plate which is turbulent just half of the time, then 


ogx,?/U = 2.077 (28 
or 
x, = 1.278 (U/ag)"” (29 
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This relation, while a basic one, will not become use- 
ful until further study (much of it experimental) per- 
mits the prediction of ¢ and g. 

The transition curves of Fig. 10 are to be considered 
as those appearing at various points along a plate 
not those at a fixed point of a plate as the velocity or 
viscosity varies the Reynolds Number. The source rate 
density g has not been assumed independent of U and 
y, and, in fact, it seems unlikely that this is so. There- 
fore, it is not correct to attempt to deduce a critical 
Reynolds Number from the present theory. In fact, 
such a Reynolds Number probably does not exist with 
any great precision. 

The mean values of the coefficients for the whole 
plate can be derived by integrating the local values 


] “2 | wi’ 3 
Cp = ai C,ds = L I (1 — f) Cy, dx + 
l “L 
f Cydx (30) 
L z re, 


This may be expressed in the form (by assuming C, = 


A Re”) 
Cp = Cp,(1 = f) + Co, f (31) 


where Cp, is the mean coefficient for the plate if the 
whole boundary layer were laminar, where Cp, is the 
mean coefficient for the plate if the whole boundary 


layer were turbulent, and where 


I od Seed 
| a Py _ 25" a I 3" A dz 


“hy 
| 
a 
“iz 
“> 
ty [ey 
< 
~~ 
a 


2° = ogx*/U 


The mean coefficient functions f and 1 — f are plotted 
in Fig. 11. The resultant mean transition coefficients 
are compared with the data of Gebers on Fig. 3. Ge- 
bers’ data were obtained on six different length plates 
and, hence, might be expected to agree better with 
the present theory (in which Reynolds Number is 
varied by varying the plate length) than Schoenherr’s 
data obtained by varying the velocity. 

A curve (not shown) can be fitted fairly well to 
Schoenherr’s 6-ft. plate results, but the 3-ft. plate 
cannot be fitted unless g increases sufficiently rapidly 
with the velocity U. 

The measurements of transition on a flat plate made 
by Schubauer and Skramstad! included hot-wire records 

A few '/s5-sec. strips of 
Fig. 12 of their N.A.C.A. 


showing the turbulent bursts. 
these records were published. 


TABLE 1 


i Nox ( Sources ) 

. ° “voxpe g , 

Ft. per Sec. y 4?/U (Bursts per Sec.) ° \Ft.? Sec. Neale 
50 )_ 0033 ~10 3.1 19.2 
105 0.0034 ~20 63.3 86.3 
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report showed bursts, as in Table 1. The correspond- 
ing maximum burst frequency computed from 


Ninax. 6 */ag"U” = 1.228 


from Fig. 9 are also shown for comparison. 

Considering that the range of g values computed 
from various data sources is from 10~* to 104 and the 
crudeness of the values of several of the constants, the 
agreement of Table 1 is fair. 

A second paper, already submitted for publication in 
the Proceedings of the First U.S. National Congress of 
Theoretical and Applied Mechanics, discusses further de- 
velopments of this transition theory and makes applica- 


tion to plates and cones. 


APPENDIX (1) 


In order to solve Eq. (2), it is convenient to introduce 
several other probabilities. Let ¢(P, x9) dxo be the 
probability that P be turbulent because of a source 
between the planes x = x) and x = x) + dxo but not 
because of a source closer to the leading edge. 


o(P, Xo) = f V(P, P») dVo dto (Al) 
S(P, xo) 


where S(P, xo) is the cross section of R produced by the 
x = x) plane. Further, let F(P, x») be the probability 
that P be turbulent because of a source between the 


leading edge and Xp. 


F(P, Xo) = I o(P, Xo') dx! “a 
0 


The probability of turbulence at P is then 


f(P) = F(P, x) (A3) 


. 


V(P, Po) dV, 
Ri 
(A2) 


Using the definition of Eq. (A2) in Eq. (2) 
V(P, Po) = g(P, Po) g (Po) {1 _ F(P, Xo) | (A4) 


By integration over S(P, xo), the x = x» cross section of 
R, we find 


o(P, x0) = / g(Po) [1 — F(P, x0)| dyo dtp) (AS) 
S(P, xo) 


where ¢(P, Po) is dropped, since it is equal to 1 in the 
area of integration. 
By differentiating Eq. (A2) partially with respect to 
Xo, we find that 
o(P, x) = OF(P, Xo) /OX0 (A6) 


Therefore, 


OF(P, xo , 
Sn =~ Pal 
OX S(P, x0) 


g( Po) dyo dto 
(A7) 


Now, Eq. (A7) may be integrated in x» directly, using 
the boundary condition [from the definition Eq. (A2)], 


F(P, 0) = 0 (AS) 
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In[1 —F(P, xo)| = — i g(Po) dVo (AQ) 
R’ 


Thus, the probability that P be turbulent because of 
a source between 0 and x» is 


F(P, Xo) =]— esr’ g(Po) dVo (A10) 


and, by Eq. (A3), the probability that P be turbulent 
iS 


f(P) = 1 — e~Sr g(Po) dVo (All) 


The other probability functions of interest now follow 
from various of the preceding equations. 

The probability density that P be turbulent because 
of sources between the planes x = x») and x = x) + 
dx) but not from some more distant source x9! < x» is 


— 0 [e~S n’8(Po) a¥ey e 
Oxo 


eS pr’ &(Po) wf g(Po) dx» dto (A12) 
S(P, xo) 


and the probability density that P be turbulent because 
of sources in dV, at Py but not any more distant point 
is [by Eqs. (A4) and (A10)] 


¥(P, Po) = oP, Py)g(Poe 2 ®’ a(Po) dVo 


oP, Xo) = 


(A13) 


APPENDIX (II)—-RELATION BETWEEN PERPENDICULAR 
Cross SECTIONS OF A CONE 


The cross-section area formed by a horizontal plane 
from a cone with an oblique axis is known, and it is 
desired to find the dimensionless cross section formed 
by a vertical plane normal to the vertical plane con- 
taining the axis. Consider a cone in the &, 7, y space 
symmetrical with respect to the y = 0 plane. The 
relation between elements of area on a horizontal cross 
section at 


7 = Ne (dA, = dy» dé») 


and a vertical cross section at 


g= (dA, = dy, dm) 
is 
& & Ne 
dy, = — dy», dn, = ( its) Al4) 
E> . fo \ ko ( 


Thus, the desired dimensionless area is given by 


1 dno dé, 
af” dm = yi 


Since 7 is arbitrary, we can choose it as we please. 
There does not seem to be any good reason for selection 


A, = (A15) 
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Arbitrarily then, the 
Thus, set 


between several possibilities. 
center of area of the x, y spot will be used. 


nm = §/BU 


where ¢ is the distance from the cone apex to the center 
of area of the spot and #U is the rate of propagation of 
the center of area of the spot—i.e., the cotangent of the 
angle between the & axis and cone axis defined by the 
center of area of the n = 7 spot. 

Finally, introduce the lateral dimension of a spot by 


(A16) 


using af in place of —&, where tan! a is the half angle 
between the two loci of successive spot positions as in 
Fig.4. Then, 


a? : dy dé a’h 
A, = : = — 7) 
1 «a {Bl BU (Alj 


A= at fdy dé/(ak)® 


Thus, o the dimensionless vertical propagation cone 


where 


cross section is given in terms of the turbulent spot 
dimensions by 


og = a’hd/B (AIS) 
The integral \ can be simplified by an approximation 
by supposing the spot to be small compared to its 
distance from its source. In this case, & ~ £ and \ 
reduces to 


A = A/(aét)? (A19) 
Thus, A is approximately a dimensionless area of a 
spot—the area A of the spot divided by the square of 


the half width. 
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RIEF REPORTS of investigations in the aeronautical sciences 

and discussions of papers published in the JOURNAL are 

presented in this special department. Publication will be com- 

pleted approximately 8 to 10 weeks after receipt of the material. 

The Editorial Committee does not hold itself responsible for the 

opinions expressed by the correspondents. Contributions should 
not exceed 800 words in length. 


A Possible Relation Between Cellular Flame 
Structure and Quenching 


Raymond Friedman 


Thermodynamics Section, 
East Pittsburgh, Pa. 


April 9, 1951 


Westinghouse Research Laboratories, 


N THE PAPER ‘Experimental and Theoretical Studies of Flame 

Front Stability’? (JOURNAL OF THE AERONAUTICAL SCIENCES, 
Vol. 18, No. 3, pp. 199-209, March, 1951), G. H. Markstein de- 
scribes observations of cellular structure in rich laminar hydro- 
carbon flames and presents a theory that predicts the following 
relationship 


Amar. (Cpp,S,/k) ~ 10 


where Amar. is the wave length for maximum instability of the 
flame front and is identified with the mean flame cell size and 
where ¢,, pu, &, and S, are the specific heat, density, thermal 
conductivity, and normal burning velocity, respectively. 

It is interesting to note that Amar. is of the same order of magni- 
tude as the quenching distance of these flames. Unfortunately, 
it is difficult to make a direct comparison, since Markstein’s 
experiments were made with air containing 30 per cent added ni- 
trogen and since quenching data are not available for such mix- 
tures. However, the present writer! has shown that the dimen- 
sionless group 


X( CpPuSu k) 


does not change greatly from one set of conditions to another. 
X is the distance between plane parallel plates which will just 
quench a laminar flame. Data are available for all the quanti- 
ties in this dimensionless quenching group for various mixtures 
of propane-air, of methane-oxygen, and of hydrogen-air; the 
numerical value of the group is found to range from 14 to 33 
(when & is calculated at 1,800°R.). This range of values is in 
good agreement with 107. 


In addition to this order-of-magnitude agreement, there are 
several other aspects of cellular flames which closely parallel 


quenching phenomena. These are: 


(A) The present writer? has found that the quenching dis- 
tance of moderately rich propane-air flames varies with the in- 
verse 0.9 power of pressure, whereas Markstein’s Fig. 5 shows 
that cell size varies as an inverse power of pressure near 0.75. 
In view of the nature of the observations, great precision is not 
possible, and this agreement appears satisfactory. 


(B) Markstein finds that replacement of nitrogen by carbon 
dioxide slightly decreases cell size, while replacement by helium 
increases cell size considerably. The effects of these diluents 
upon quenching distance are found! to be similar. 


Markstein reports that cell size is constant for fuel con 
centrations between 1.2 and 1.4 times stoichiometric and in 
creases gradually for richer compositions. Quenching distance 
beiuaves similarly, the minimum distance occurring on the rich 


(C) 


side. 


(D) Cell size was found to decrease markedly with increasing 
molecular weight at fuel concentration near 1.3 times stoichio- 
Since no quenching data are available for heavier hydro 
is necessary to consider minimum spark-ignition 
this energy is roughly proportional to 
Fig. 4 of reference 3 shows 


metric. 
carbons, it 
energy measurements ;° 
the square of quenching distance.‘ 
that the minimum value of the minimum energy vs. composition 
plot shifts progressively toward richer mixtures for six hydrocar 
bons ranging from methane to heptane, even though the mini- 
mum value of the minimum energy is nearly the same for all of 
the hydrocarbons except methane. If such a shift were also to 
take place in cell size vs. composition plots, it could conceivably 
“explain” the effect observed by Markstein, although the reason 
for this shift is still obscure. Experiments are in order to de- 
termine how quenching distance vs. composition curves depend 
upon molecular weight of fuel. 


The foregoing arguments [with the possible exception of (D)| 
show rather convincingly that both flame quenching and cellular 
instability are related in similar ways to the same basic param- 
eters, particularly to the reaction zone thickness. 
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The Combined Supersonic-Hypersonic 
Similarity Rule 


Milton D. Van Dyke 
Ames Aeronautical Laboratory, N.A.C.A., Moffett Field, Calif. 


April 18, 1951 


HE HYPERSONIC SIMILARITY RULE of Tsien' and Hayes? can 

be combined with the similarity rule for linearized supersonic 
flow. There results a single rule valid throughout the range of 
supersonic speeds from just above the transonic range up to in- 
finite Mach Number and for general three-dimensional slender 
bodies. Thus, the value of the rule for correlating experimental 
results is considerably increased. 

Similarity rules of the sort considered here correlate com- 
pressible flows past slender bodies that differ from one another 
only by a uniform expansion or contraction of all dimensions 
normal to the free-stream direction. For example, the second 
body in Fig. 1 is derived from the first by reducing its thickness, 
angle of attack, and aspect ratio by the factor r2/7,, where r is 
some representative thickness parameter. (In practice, it seems 
best to take 7 to be the maximum streamwise slope of the body 
surface. ) 

At hypersonic speeds, as the thickness is reduced, flow similarity 
is maintained if the free-stream Mach Number M is simultane 
ously increased so that the product 7M remains constant. Then, 
at corresponding points Q (Fig. 1), the ratio of local pressure co 
efficient C, to 7? remains constant. Hence, the hypersonic 
similarity rule may be written 


Cp = r*P,(rM) (1) 
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where / is some function of the hypersonic similarity parameter 
7M. The range of validity of this rule has recently been in- 
vestigated both analytically*® and experimentally.' 

At moderate supersonic speeds, linearized theory yields a dif- 

/ 

ferent rule. The similarity parameter is 7 M? — 1, and cor- 
responding to Eq. (1) the supersonic similarity rule may be 
written 


Cy = P(r M? — 1) (2) 


These two separate rules can be combined by adopting the 
latter form, because the difference between MM and V i? — lis 
(This can be justified rigorously 
Therefore, Eq. (2) 


negligible in hypersonic flow. 
by re-examining the derivation of the rule.) 
expresses a combined supersonic-hypersonic similarity rule valid 
at all speeds above the transonic range. 


To illustrate the superiority of the combined rule, suppose that 
the pressures on a cone of 5° semivertex angle are known (from 
the M.I.T. Tables*) and that approximate values for a 10° cone 
are derived therefrom using the similarity rule. Fig. 2 shows 
that the results of applying the combined rule almost coincide 
with the exact values down to the Mach Number at which sonic 
velocity is attained at the surface, while the hypersonic rule be- 
comes inaccurate at considerably higher speeds 


The combined similarity rule is obtained at the expense of some 
loss of generality in the supersonic rule. Supersonic similarity is 
independent of the adiabatic exponent y, while for the combined 


rule y must remain fixed. Another penalty attached to the com 
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bined rule is that an alternative form of the hypersonic rule 
which relates local to free-stream pressure, 


p/po = F(rM (3 
loses its simplicity, becoming 
p/p. = 1 + r2M°G(4rV/ M? — 1) (4 


Recently, Smelt® has proposed an alternative similarity pa 
rameter valid at both supersonic and hypersonic speeds. His 
rule, analogous to Eq. (4), is 


p/pbo = F(rM2/V/ M2 — 1) (5 


However, it can be shown that this rule is restricted to two 
dimensional flow or flow past planar systems such as flat wings 
(for which the similarity parameter is not uniquely determined 
whereas the rule expressed by Eq. (2) applies to general thre 
dimensional flow. 


The combined similarity rule may be expressed in other forms. 
such as 

Cp = [1/(M? — 1)| P2(r7V M — 1) (6 

but the first form |Eq. (2)] seems preferable because only in this 

case does the function P remain finite at infinite Mach Number 


The corresponding rules for drag coeflicient Cy and lift coef 
cient Cr are 


Cp = r2t*D(7rV M2? — 1) 7a 
Cr = rit*L(irY M2 — 1) (7b 
where PD and ZL are functions of the similarity parameter 


r\/ M? — 1,k = 
and k = 


pressure position is a function only of the similarity parameter 


1 if some plan-form area is taken for referenc 


0 if some cross-sectional area is used. The center of 


At high Mach Numbers, the hypersonic (or combined) simi 
larity rule applies to rather thick bodies in some cases if 7 is 
taken to be the sine of the maximum angle of inclination of the 
body surface. The reason for this is that, for bodies of any thick 
ness at J = o (and y = 1), Newtonian theory predicts pressure 
coefficients proportional to the square of the sine of the angle of 


inclination. 


Just as the supersonic and hypersonic similarity rules are com 
bined here, Spreiter? has shown that for two-dimensional or 
planar bodies the subsonic, transonic, and supersonic similarity 
rules can be combined. However, it does not seem possible, even 


for planar bodies, to combine all four régimes. 
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On the Velocity of Sound in Fluids 


C. Truesdell 
Professor of Mathematics, Graduate Institute for Applied Mathe- 
matics, Indiana University, Bloomington, Ind. 


April 23, 1951 


OT LONG AGO there was a protracted discussion of the 
N velocity of sound, begun by Professor Morduchow’s! 
objection to the “lack of rigor” in the customary derivations 
and carried on by Professors Morkovin,? Morduchow,* and 
Millikan,’ and by Dr. Munk.’ None of these authors saw fit 
to define the velocity of sound before deriving its value; most of 
the discussion assumed the existence of a barotropic relation 
pb = fle). The present no original 
contribution; its only purpose is to call attention to the basic 
known results on the subject, which were not mentioned by any 


communication makes 


of the four disputants. 


For an inviscid fluid, the local velocity of sound c may be 
defined as the speed of advance, relative to the particles instan 
taneously situate upon it and in a direction normal to itself, 
of a surface (not necessarily plane) across which the velocity, 
pressure, and density are continuous, but at least one of the nine 
components of the velocity gradient suffers a discontinuity. 
In the terminology of Hadamard,® such a surface is called a ‘‘wave 
Appell’ calls such a wave a ‘“‘sound wave.” 
f(p), it was shown by Hugoniot* in full 
Of the customary infinitesimal 


of second order’’; 
In the case when p = 
generality that c? = dp/dp. 
arguments he says: ‘‘|Treatments of] the theory of propagation 
of motion in a fluid , have introduced into the equations of 
hydrodynamics hypotheses, disguised, it is true, under the name 
of approximations, but singularly altering the worth of the re- 
sults which can be deduced from them."’ His analysis is repro- 
duced (in two different forms) in the texts of Webster,’ Appell,’ 
and Lichtenstein," and the matter is treated in detail by Hada- 
mard in his great classical treatise on waves."! 


A part of Professor Morkovin's note? and all of Dr. Munk’s 
note® are devoted to the more interesting case when a barotropic 
relation p = f(p) does not exist. According to Dr. Munk, 
“the necessity to give a Mach angle uw leading to physically 
Op/ Op, the derivative 


significant Mach lines’’ indicates that c? = 
being taken ‘‘in the direction of the stream line of a steady flow. 
It is hard to take this proposal seriously, since, quite apart 
from the purely physical objection that no mention or use what- 
ever is made of any property of sound waves, the classical 
result, c2? = dp/dp when p = P(p), is not even included in all 
cases when it holds. The statement fails to have a meaning 
when the flow is three-dimensional or unsteady, and the adduced 
arguments do not hold if the flow is subsonic, there being then 
no real characteristics, physically significant or otherwise. Pro- 


fessor Morkovin’s proposal is 
c? = (Op/ Op)s 


s being the entropy, but he presents no analysis and makes no 
attempt to discuss the propagation of motion. 


Now in fact the second derivation of Hugoniot!? remains 
valid even if the relation p = p(p) has a form varying from one 
fluid particle to another. Such a relation would be of the form 


p = f(o, X, Y, Z) 


where Y, Y, Z are “Lagrangian” coordinates. In a flow of an 
inviscid fluid when heat conduction is neglected, the entropy 
remains constant for each particle, and hence by Hugoniot’s 
result there is indeed a definite velocity of sound, and it is given 
by Professor Morkovin’s formula. However, in a flow where 
each fluid particle conserved, for example, its own temperature 


6, the speed of sound would be given by 


SO 


FORUM 


c? = (Op/Op)e 


If we really wish the term ‘“‘speed of sound”’ to represent the 
speed of advance of a physical sound wave (not necessarily 
plane), we should recognize the fact that there is no physical 
reason a privri to expect that speed to be a function of the ther- 
modynamic state only. Hugoniot’s result yields a broad class of 
flows in which it /s a function of state, but probably one could 
show that if there is no barotropic relation p = f(p, X, VY, Z), 
then more information than simply the values of p, 8, X, }, Z is 
required in order to determine the value of « 

A closer representation of the physical situation is the quasi 
wave of Duhem,'* which is a layer of small thickness, ¢, across 
which the change in pressure, density, and velocity is of the order 
of «, but at least one of the nine components of the velocity 
gradient suffers a large (though possibly continuous) change 
In a remarkable analysis, which has apparently escaped the notice 
of all current writers on this subject, Duhem shows that in the 
case of a slightly conducting perfect fluid any approximately 
transversal quasi-wave is stationary, but there are two types of 
approximately longitudinal quasi-waves. One, of thickness 
small in comparison to the coefficient of conductivity, propagates 
approximately according the Newtonian 
(O0p/Op)e, while the other, of thickness of the same order as the 
coefficient of conductivity, propagates approximately according 
to the Laplacian formula c? = (O0p/0p)s.. He concludes further 
that when the effect of viscosity is taken into account the fore 


to formula c? = 


going conclusions are not modified, provided that the thickness 
of the quasi-wave be of the same order as the viscosity coef 
ficients. Since these statements cannot be expressed in dimen 
sionless terms, it is clear that they can be given a precise meaning 
only in terms of a limiting process « > 0. 

If, finally, we have no interest in sound waves but merely wish 
a local reference quantity, endowed with the dimensions of 
velocity, to be used in scaling laws (e.g., in definition of the 
Mach Number), then any sufficiently general definition is admis 
sible, providing it be used consistently Thus Bateman'* 
regarded the term /ocal velocity of sound as ‘‘merely a convenient 
name.’ Two different such ‘velocities of sound’ occur in the 


work of Hicks.'® 
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The Analogy Between Irrotational Gas Flow 
and Minimal Surfaces 


C. Truesdell 
Professor of Mathematics, Graduate Institute for Applied 
Mathematics, Indiana University, Bloomington, Ind. 


April 23, 1951 


_ A RECENT NOTE, Professor Szebehely! points out an analogy 
between an open soap bubble (minimal surface) and an ir- 
rotational flow of a perfect fluid. 
if and only if the pressure-density relation is of the Earnshaw- 
Chaplygin type p = A — (B/p) was discovered by Chaplygin? 
and Bateman,’ is mentioned in Sauer’s textbook,‘ and has been 
made the subject of extensive study.® 


This analogy [exactly valid 
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On Broszko's Theory of Turbulent Flow 


William Squire 
Aerodynamic Research Department, Cornell Aeronautical 
Laboratory, Inc., Buffalo, N.Y. 


April 19, 1951 


IT A LITTLE KNOWN PAPER, Broszko! attempted to formulate a 
rational theory of turbulent flow by deriving additional rela- 
tions between the mean motion and the turbulent fluctuations. 
No discussion of the theory seems to have been published; 
therefore, it is believed that some critical comments on the deri- 
vation and consequences of the additional equations may be of in- 
terest 
Broszko employs Reynolds’ equations 


Ow: , dit, 
or Ox? 


1 dp o*, 


Ou;'u;’ 
=~¢ am 
p Ox; Ox Ox? ox? 


(1) 


with ¢ independent of position but much larger than the ordinary 
kinematic viscosity. In order for the results of the theory to 
agree with experiment, « must vary with Reynolds Number. 
The averaging process is defined in Reynolds’ manner as 


g= SedV/fdv (2) 


where V is a small fixed volume around the point considered. 
Aside from such obvious properties as linearity and commutation 
with partial differentiation, the characteristic property of aver- 


aging is 
(U; + C;) uj’ = u,'u;’ (3) 


where U; = 4; + u;' and C; isa constant velocity. This implies 


that both 


ui;’ =0 u;u;’ = 0 (4) 


j ; 

Broszko attempts to derive the additional equations from Eq. 
(3) rather than from a physical principle. Differentiating both 
sides of Eq. (3) with respect to time gives 
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1 " ae ’ ee) 1 _- 
V aut Ci) u;'dV => i u;'u;’ dl (5 
Using the continuity equation, 
dg Og e.. P 
= ~ (tte + Ux" dg (6 


dt ot ax 


and setting the average of terms containing only one fluctuating 
quantity equal to zero, leads to the set of nine equations. 
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There are two principal difficulties in the derivation. First, 
the differentiation of the integrals in Eq. (5) is derived? for mate- 
rial volumes, not the fixed volumes of Eq. (3). Second, it seems 


equally plausible to start with 
Mf’ = 0 (8 

where M is any mean quantity and /’ is any fluctuating quantity 
and to derive 

OMf' ux! /dx* = () (9 
For example, if M/ is a constant and f’ = u;’, the result is 

Ou j'UK’ /ox* = 0 (10 
which, in conjunction with Eq. (1), leads to Boussinesq’s early 
theory of turbulent flow. 

By considering flow through a long pipe of constant cross sec- 
tion, the equations are greatly simplified, since there is only one 
mean velocity component and no variation in the direction of 
The simplified set corresponding to Eq. (7) is satisfied by 


yt a wera yy’ « f1/(U+C)] (11 


flow. 


However, x’? is arbitrary, and, therefore, 


K = u'0'/V uty"? 


is also undetermined. Following Broszko, —C is equal to the 
apparent wall velocity, but infinite values do not appear because 
of transition to the laminar sublayer. The theoretical »’? and 
w’? distribution agrees qualitatively with experiment, but u'v’ 
does not vanish at the center as required by current theory. To 
clarify the problem, consider the equation of motion for a chan- 
nel: 
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7 27° > 
5 a. (12) 
oy oy? Ox d 


Integrating from 0 to y gives 


, 230 ou Toy ; 
(ou'v’ — pu'v'y=0) — (v« oe = =r (13) 





€ 
oy OVy =0 d 
At the center 0u/Oy = 0, and it has been assumed that u’y’ is 
zero also so that 
- os oU - . 
ee ee ee = pu'v (14) 
oy 


since the 0u/Oy term is negligible if « is the usual kinematic vis- 
This was not confirmed by experimental results until the 
se, 3 The earlier results of Wattendorf‘ 
While Laufer’s results show 


cosity. 
recent work of Laufer. 
and Reichardt® assumed Eq. (14). 
discrepancies up to 20 per cent from Eq. (13), the deviation from 
Eq. (11) is even greater, especially at the center. 

In spite of its deficiencies as a complete rational theory of turbu- 
lent flow, Broszko’s theory gives satisfactory results for the mean 
velocity distribution. The writer has obtained the following ap- 
proximate solutions: 


(Om — 0)/U* = 122(1 — V1 — 9) (15a) 
(Om 0)/U* = 3.2n? — 2.51n (1 — 7?) (15b) 
i - [i—¥* 
(Um — U)/U* = 4.9n? — 2.5 [1m ( ) + 2n (15c) 
1+ 


where 9 is r/radius for a circular pipe and y/half width for a 
channel. Eq. (15a) is derived for either case, Eq. (15b) is de- 
rived for a circular pipe, and Eq. (15c) is derived for a channel. 
However, the three expressions give almost identical distribu- 
tions. Fig. 1 compares Eq. (15a) with averaged experimental 
results 
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Equivalence of the Spanwise Lift Distribution 
to the Spanwise Lift-Influence Function for 
Slender Wings and Wing Bodies 


G. K. Morikawa and A. E. Puckett 
Research and Development Laboratories, Hughes Aircraft Company 
April 23, 1951 


i ines USEFULNESS of the concept of the spanwise lift-influence 
function* is apparent in the calculation of lift and rolling 
moment of lifting surfaces with twist (real or apparent). A 
particularly simple result has been obtained by Lagerstrom and 
Van Dyke! (also independently by Ward?), as expressed by their 





* Defined as the total lift due to a chordwise strip of unit width deflected 
at unit angle of attack (see Fig. 1). 
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theorem on preservation of lift for wings with supersonic leading 
edges and straight supersonic trailing edges. The only com- 
parable result available for plan forms with subsonic leading 
edges is given by the modified strip method of Graham,’ appli- 
cable to cases where there is no interaction between edges 

Within the realm of slender-body theory, expressions for the 
spanwise lift-influence function are readily available for both 
wings and wing bodies. The analysis for this limiting case is 
reduced to a two-dimensional potential problem in the cross-flow 
plane,‘ and the problem of finding the lift-influence function is 
analogous to that of obtaining the Green’s function. The solu 
tion for the velocity potential is simply obtained by considering 
symmetrical wings with respect to the y-axis (Fig. 1) and trans 
forming the wing-body cross section into the unit circle 

The complex velocity potential is 
w=¢e+wWe= Q In 


2r 


J V1 — (4/A?) (s + a?/z)? + V1 — (4/A?) (x + a? xo)? 


/ 
(V1 — (4/A® (2 + a/s)? — V1 — (4/ A®) (x + a/xo)?f 
(1) 
where 
Q = VaAx, source strength 
a = body radius 
x9 = position of deflected strip of width Ax 
A = (1/2) (s + a?/s) 
s = semispan of wing 
The portion of the lift acting on the wing is 
4Ayk\™ 8 ‘ 
-—— ¢ dx (2a) 
q V Ja 














a) PLANFORM 





b)CROSSFLOW PLANE = 





FIGURE 1. 
WING -BODY CONFIGURATION WITH SYMMETRICAL 
CHORDWISE LIFTING STRIPS 
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AL (Ww so AL (B 
qa Ax 


where ¢ is obtained from Eq. (1) by setting zs = x, and the por- AL 
tion of the lift acting on the body is 


qa Ax 
AL“ 8 a = 4V (s + a*/s)? — (xo + a?/xo)? (4) 

- gy dx (2b) 
q J 0 which is exactly the spanwise lift distribution. This is a verifj. 
cation of a result that is directly obtainable from the general re. 


verse-flow theorem given by Ward? (also independently by W. D, 


6 


where ¢ is obtained from Eq. (1) by setting = ae” and a cos 


6=Xx. 
Integrating Eqs. (2a) and (2b) gives 


4 | l 
= AB — 2AB |.cos~'4/1 — 
T " (co V A? 


/ 
/A?—~1+AB 
In » = AY LF x 
me eee 
B+ V(A? — 1)(1 — B) 


B — V(A? — 1) (1 — B?) 


Al (W) 
qa Ax 


(3a) 


2AB cos 


/A*—1+AB 
In| — +AVI1 
V A?—1-—AB 


B+ VV (A? — 1) (1 — B?) 


— B?x 


In (3b) 


B — V(A? — 1)(1 — B?) 


ii 1 (Xo + Chae 
\ (s + a?/s) 


Then the spanwise lift-influence function is 


where 


Hayes* whose proof follows that of Munk® and Brown’ for less 
general reverse-flow theorems) or an equivalent statement given 
by Jones’ and some intuitive arguments (on slender-body con. 
those presented by Lagerstrom and 


figurations) such as 


Graham.? 
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Three-Dimensional Shock-Wave Reflections 


(Continued from page 489) 


wall, also parallel to the free stream (Fig. 6), turns out 
to be a complete dual of the preceding case of the plane 
shock impinging on a cylinder. One reason for this 
duality stems from the already mentioned fact that 
the reflection conditions, and therefore the geometry of 
various planes as in Figs. 3 and 5, have a purely local 
character; they are associated with the orientation of 
the tangent line t rather than with the shape of the 
impingement curve i. The lack of similarity between 
the impingement hyperbola located in the wall plane 
here and the impingement ellipse located in the shock 
plane in the preceding case is immaterial. The other 
reason for the duality lies in the fact that the conical 
shock, like the plane one, has constant strength every- 
where. Thus, the normal component of the Mach 
vector Mo sin 6 again has a constant magnitude c and 
Eq. (4) remains valid. 


Referring to Fig. 6, one identifies the wall plane 
which contains Mp, and t as the tangent T plane, and 
the plane through the cone generator OP and t as the 
shock plane S. The shock angle is measured in the 
plane through M/) and OP, which is therefore equivalent 
to the plane V. This plane is no longer vertical; on 
the other hand, the tangent plane T has become hori- 
zontal. As a matter of fact, one obtains the complete 
correspondence between the geometry of Fig."5 and the 
present geometry by rotating the former configuration 
around the M]po vector as axis until the plane 7 becomes 
Inspection of Fig. 4 shows that the needed 
angle of rotation is exactly w and that n, the line nor- 


horizontal. 


mal to T, becomes coincident with the local vertical 
axis. Thus, the relations between the edges of the 
tetrahedron PABC of Fig. 5 remain valid, and the Eq. 
(5) relating 6, to the position parameter w also holds. 
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